Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,121,547 G→A intergenic (‑60/+14) dinI ← / ← pyrC DNA damage‑inducible protein I/dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,121,5470GA94.4% 44.5 / ‑3.8 18intergenic (‑60/+14)dinI/pyrCDNA damage‑inducible protein I/dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (10/7);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTGCCCAGCGCGGTTG  >  NC_000913/1121315‑1121790
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
tcttTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAATTCTCCCGCCACGGCGTCAATAACCCCCGCCGGCCATTGGGAAGTTTTTGCTATTGTTTCTTCACTTTTCATAATATCCCCCCTTTTTATATACTTGTTTTTTTTTCAGGTTAAATAAACTCTTGTCAACAAGAAGCTActttttttttt                                                                                                                                                                                                                                    >  1:285711/1‑241 (MQ=255)
                                                                    aTCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAt                                                                                                                                                                 <  1:43448/250‑1 (MQ=255)
                                                                    aTCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATg                                                                                                                                                                <  1:6421/251‑1 (MQ=255)
                                                                    aTCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATg                                                                                                                                                                <  1:297063/251‑1 (MQ=255)
                                                                                             cttgAAAGTTCGCCAGTCAGGGCGTCAATAGCCCCAGCTGACAATGGCGACGTTTTCGCTATGGTGCCTTCAATTCGCATCATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAACTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGCACCAGCGTACAGTTTCCCCGTCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGGTGCTCTTCAc                                                                                                                                       >  2:108522/4‑251 (MQ=255)
                                                                                                    gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGACGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTAcc                                                                                                                                  >  1:112384/1‑249 (MQ=255)
                                                                                                     ttCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGt                                                                                                                               >  1:262929/1‑251 (MQ=255)
                                                                                                         taagaaagtgcattaaTAGCCCAAGCTGGCAATGGACAATATTTCGCTATGGGGACTTCAATTCGAATAATAGCCCCCTGATGAATACAATGGGTATTTATACCGGTAAAATAACCTAAAGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCGGCGTGTCATCAGTCATTGCGATGCTTTCAGCAACCTCTTGGTCTTCACGTACCAGtt                                                                                                                              >  2:280109/15‑248 (MQ=255)
                                                                                                                          aGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGAc                                                                                                          <  1:10824/250‑1 (MQ=255)
                                                                                                                                              aagatttCGATATGGTGACTTCAAATCGCATAATAGCCACCTATTGAATCTACTGGTTATTTATACAGGTAAAATAACCTACTGACAGCAAGAAGCTACGATTTTTATTATTTAACGGACCAGCGTACCTTTTCCCCGGCGAGGGTAGGCACAAACGTTTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACTGACCAGTTCGTTGAATGGTACGTAGACCGGCAACCCATAGaaatgagg                                                                                     >  2:182302/5‑245 (MQ=255)
                                                                                                                                                           ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTacaga                                                                         >  1:18214/1‑250 (MQ=255)
                                                                                                                                                                          catattcacccacTGTTGAATAGACTGTTTATTTATACAGGTATAATAACCTAATGCCAAGTAGAAGCTACTATTTTTATTGTGTAACGGTTCAGCGTACCTTTTCCCCGTCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAGGCTTTCAGCAACCTGTTGCTCTTCACGTTCCAACTCGATGAATGTGTCGTTGGCCGGCACCCCATAAAACTGCTGGCCGTTTACAGAACAGAATGTTTCAAAg                                                         <  2:4198/239‑1 (MQ=255)
                                                                                                                                                                                    cccTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCACAGTGCTGCaa                                                 >  1:250310/1‑249 (MQ=255)
                                                                                                                                                                                           gAATATACTGGTTATTTATACAGGTATAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGcattca                                         <  2:262929/250‑1 (MQ=255)
                                                                                                                                                                                                               cAGGTAAAATAACCTAATGACAACAAGAAGCTACTATTGTTATTTGTTAACGGACCATCTTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGTTGCTTTCACCTACCTGTTACTCTTCACGTTCCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCTGGCCGTTTACAGAACCGAATTCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGt                    <  2:72954/251‑1 (MQ=255)
                                                                                                                                                                                                                   tAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACt                >  1:214918/1‑251 (MQ=255)
                                                                                                                                                                                                                                tGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTGCCCAGCGCGGtt   >  1:278524/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTg  >  1:94840/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
TCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTGCCCAGCGCGGTTG  >  NC_000913/1121315‑1121790

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: