Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,578 C→T N30N (AAC→AAT nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,5780CT85.7% 24.2 / ‑0.9 14N30N (AAC→AATnhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (6/6);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

CATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCA  >  NC_000913/17351‑17815
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
cATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATTGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCTCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTATCAGTTATTCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTTGCGATTATTATGGCCAATAGCGGCGTAACCATTGGATGGTa                                                                                                                                                                                                                         <  2:220286/251‑1 (MQ=255)
               ggCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACCCTATAATCTGATTTTAACGATGATTTGTTCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGCTGATTATGGCCAATAGCGACGCAACCAGTGGATGGTATCACGTCTTTCTgg                                                                                                                                                                                                           >  1:208271/1‑250 (MQ=255)
                                      tAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGggttg                                                                                                                                                                                   >  1:247118/1‑251 (MQ=255)
                                                                    ggagtttAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAATTGAAACATCTGTATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCACTGCCGCTATCCTGGCGATGATCATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTGCGCTCGAAATCAACAAAAACATGCTGt                                                                                                                                                     <  1:293815/247‑1 (MQ=255)
                                                                                               aaaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAACAACATGCTGTTATGGATAAATGACGCGCTGATGGCgg                                                                                                                          >  1:49334/1‑251 (MQ=255)
                                                                                                       aaaaaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACTTCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGACCAATAGCGGCGCAACTAGTGGATGGTCTCAAGTCTTTCTGGAGCCGACGGTTCACCTCCGGGTTGGTTCACTCGAACTCAACAAACACATACTGTTATGGATAAATGACGCGCTGATGGCGGTAtttttt                                                                                                                  >  1:200302/1‑250 (MQ=255)
                                                                                                                                       aaaGTGAATCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACt                                                                                 <  2:66984/251‑1 (MQ=255)
                                                                                                                                             aaaCATCTGCATCTATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCACCCCGTGGATGGTATCACGACTTTCTGGAGACGCAGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGAGGTTATGGAAAAATGACGCCCTGATGGAGGAGTTTTTCCTGTTTGTAGGTGTGGAAGTTAACCGTTAAATGATGc                                                                            >  2:30785/1‑250 (MQ=255)
                                                                                                                                                           aTTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGCGGCCGGTTCAGCTCCGGGGTGGTTCACTCGAAATAAACAAAAAAATGCTGTTATGGATAAATGACGCGCTGATGGCGGTGTTTTTCCTGTTAGTCAGGTTGGAAGTTAAAAGCGAAACGATGCAAGGATGCCTAGccg                                                             >  2:365464/1‑250 (MQ=255)
                                                                                                                                                                        gATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGTTATCACGTCTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGTAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGc                                                <  1:314108/251‑1 (MQ=255)
                                                                                                                                                                                               tATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCCTTTCCAGTTATCGCCGCt                          >  1:2274/1‑250 (MQ=255)
                                                                                                                                                                                                         gCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGgat               <  1:285250/251‑1 (MQ=255)
                                                                                                                                                                                                                      tgatTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGcc  >  1:255254/1‑250 (MQ=255)
                                                                                                                                                                                                                      ggatTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCa  <  1:365464/250‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              
CATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCA  >  NC_000913/17351‑17815

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