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,5470GA95.5% 60.9 / ‑5.4 22intergenic (‑60/+14)dinI/pyrCDNA damage‑inducible protein I/dihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/11);  minor base T (0/1);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTT  >  NC_000913/1121307‑1121757
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                  
gTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGTATAACGTTCCGATACGTGTCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTGTCGCGATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTATAATAACCTAATGACAACATGAATCTACg                                                                                                                                                                                                           <  2:413332/250‑1 (MQ=255)
                                                            aTACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                                                                                                                              <  2:156313/251‑1 (MQ=255)
                                                               cGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCATCGTACCGTTTCCCCGGCGAGGAATGGCTCCAGCGTGtcatc                                                                                                                                           <  2:532274/251‑1 (MQ=255)
                                                                gTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGGTCGCCAGCCAGGGCGGCAAAAGCCCCAGCTGGCAAAGGAGAAGTTTTTGCTTTGGGGACTTCAATTCGCATAAAAACCCCCCGTTGGATATACTGGTTATTTATACAGGGAAAATAAACTAATGACAACAAGAAACTACGATTTTTTTTTTTTTACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatc                                                                                                                                           >  1:192137/1‑250 (MQ=255)
                                                                       tCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATcagtcagt                                                                                                                                    <  2:99105/250‑1 (MQ=255)
                                                                                      gCATACTGAATACGGCGGGAAAGTTCGCCATCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCATGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGGTTCCCCGGCGAGGAATGGAACCAGCGTTGCATCAGTCATTGCGATACTTTCAGCaa                                                                                                                    <  2:516435/251‑1 (MQ=255)
                                                                                                       ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGt                                                                                                   <  1:485812/251‑1 (MQ=255)
                                                                                                         aaaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTAc                                                                                                 <  2:167107/251‑1 (MQ=255)
                                                                                                                      cagGGCGTCAATAGCCCCAGATGGCAATGGTGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGGTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGCGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTGGCTCTTCACGTACCAGTTCGATGAAt                                                                                    >  2:437106/1‑251 (MQ=255)
                                                                                                                            gTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCATTTCGATGAATGTGTCg                                                                              >  2:544682/1‑251 (MQ=255)
                                                                                                                                             gCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGATAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATa                                                             <  1:194615/251‑1 (MQ=255)
                                                                                                                                             gCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATa                                                             <  1:178593/251‑1 (MQ=255)
                                                                                                                                                  ggAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACt                                                        >  2:152730/1‑251 (MQ=255)
                                                                                                                                                          tttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACCGGTAAACTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGAt                                                                                        >  2:290282/1‑211 (MQ=255)
                                                                                                                                                          tttCGCTATGGTGACTTCAATGCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGAt                                                                                        <  1:290282/211‑1 (MQ=255)
                                                                                                                                                             cGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGttt                                             >  2:351221/1‑251 (MQ=255)
                                                                                                                                                                  tGGTGACTTCACTTCGCATCATAGCCCCCTGTTGAATATCCTGGTTATTTATACAGGCAAAATAACCTAATGCCAACAAGACGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTACCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGCATGTGTCGTTGACCGGCAACCCATAGCACTGCGGGCCGTTTac                                           >  2:317766/1‑248 (MQ=255)
                                                                                                                                                                     tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACAAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTacagaac                                      >  2:349778/1‑250 (MQ=255)
                                                                                                                                                                            aaTTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACAAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCt                              >  1:40118/1‑251 (MQ=255)
                                                                                                                                                                             aTTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATTCtt                             >  2:343208/1‑251 (MQ=255)
                                                                                                                                                                                  cATAATAGCCCCCTGTTGAATATACTTGTTATTTATACAGGTAAAATATCCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAGCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATATAACTGCGGGCCGTTTACAGAACATAATGCTTCAAAg                        <  2:394562/251‑1 (MQ=255)
                                                                                                                                                                                                         aCTGGTTATTTAGACAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCATCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTTCTCGTCACGTACCAGGTCGATGAATGTGTCGTTGACCGGCAACACATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATttctt  <  1:374666/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                  
GTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTT  >  NC_000913/1121307‑1121757

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: