Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 18,477 T→C V330A (GTG→GCG)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,4770TC96.0% 73.1 / ‑6.5 25V330A (GTG→GCG) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (3/21);  total (3/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGA  >  NC_000913/18241‑18718
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
ggAGCCTGTGTTAGACCCGTGGGTGGGGTATCTGATTTTGCCGCTTGTTT‑CATTTGCTAATGCTGGCGTTTCACGGCAAGGTGTCCCGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCtgt                                                                                                                                                                                                                                     <  2:237485/251‑1 (MQ=255)
ggAGCATGTGTTACACCCGTGGGAGGCGTATCTGATTTTGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTTGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTgggggatgtgt                                                                                                                                                                                                                                     >  1:286349/1‑244 (MQ=255)
            acACCCGTGGGTGGCGTATCTGATTTTGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGtt                                                                                                                                                                                                                         <  1:449121/250‑1 (MQ=255)
                         gCGTATCTGATTTTGCCGCT‑GTGTGCATTTGCTAATGCTGGCGTTTCACGGCAAGGTGTCACGCTGGATGGGTGGACCTCCATGCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCGGTTCTGCGGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGGGGGGATCCTGTGCGGTATCGGGTTTACTATGTCTAt                                                                                                                                                                                                            <  2:282436/251‑1 (MQ=255)
                              tCTGATTTTGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCCGTGCGGGATCGGTTTTTTTATGTCtttttttt                                                                                                                                                                                                       >  1:64736/1‑244 (MQ=255)
                                     ttGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGTCTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCa                                                                                                                                                                                                 <  1:91719/250‑1 (MQ=255)
                                      tGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAg                                                                                                                                                                                                <  1:353803/250‑1 (MQ=255)
                                                      ttGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAg                                                                                                                                                                                 <  2:121523/250‑1 (MQ=255)
                                                            aTGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGAt                                                                                                                                                                          <  1:264479/251‑1 (MQ=255)
                                                                                               ggCTTGCCCCCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                                                                                                                                        <  2:37044/250‑1 (MQ=255)
                                                                                                                    ttGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGGGCCTGTGCGGGATCGGTTTTTTTATTTCTATCTTTTTTTCCAGCCTGGCCTTTTGTTGCGTAGATACAGCACTTATTAACTGGGCGAAACTAGGGATTCTGGTCgggtttttttcttcggc                                                                                                                  >  1:13129/1‑236 (MQ=255)
                                                                                                                                                 cAAACCGCTGGGGATTAGGCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGcgttc                                                                                     <  2:280983/251‑1 (MQ=255)
                                                                                                                                                  aaaCCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGTTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGcgttcg                                                                                    <  2:322053/251‑1 (MQ=255)
                                                                                                                                                   aaCCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGcgttcgt                                                                                   <  1:5397/251‑1 (MQ=255)
                                                                                                                                                            ggATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTTGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCa                                                                          <  1:473159/251‑1 (MQ=255)
                                                                                                                                                                           gctgGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAgg                                                             <  1:34204/249‑1 (MQ=255)
                                                                                                                                                                                  ggcgtaggggtggaaacagggtcttCTGCCTGTGGGAACGTGTTAGCAGCAAATTATTGCGTTGGGGATCCGGTGCGGTATCGGTTTTACTATGTCTTGCTTTAGTGCCAGCCAGGCCTTTTGTAGCGTAGTTCCAGAACTGATTAACTGGGCGAAACTCGGTATCTAGGTTGGTGTTATCTCTTCGGCGGTAATTTGATACAGCTGGTGACGTGTTCGTTTGCGGCCATCAGTTTGTCAGGACGTTTTa                                                     <  2:487514/226‑1 (MQ=255)
                                                                                                                                                                                   gctcTGCGTTTGTAACTGGCGCATCTGCCTGAGTGTACGCCTTATCAGCAAATTATGGCGGTGGGGTTCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCATCCTGGCCTTTTGTAGCGTAGATCCAGAACTGACTAACTTGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACATCTGGTTACGCGTTCGTTTGCGACCATCAGTTTGACATGACTGTTTAcc                                                   <  1:109267/248‑1 (MQ=255)
                                                                                                                                                                                                aacgggCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATTGCGGTGGGGATACTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATaaa                                      <  1:18921/247‑1 (MQ=255)
                                                                                                                                                                                                             gCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTc                         <  1:12370/251‑1 (MQ=255)
                                                                                                                                                                                                                 gAGGGAACGACTTATCAGCAAATTATGGCGGTGTGGATCCTGTGCGTTATCGGGTTTACTATGGCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTTGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCAtt                      <  1:343622/250‑1 (MQ=255)
                                                                                                                                                                                                                    ggAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTa                  <  1:423143/251‑1 (MQ=255)
                                                                                                                                                                                                                     gAACGTCTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTAt                 <  2:49480/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gaaatttATGGTGGTGGGGATCCTGTGCGGTTTCTGTTTTACTATGCCTACCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGACCTGATTAACGGTGCGAAACTCGGGATCCTGGTCGGTTTTATCTCTTCGGCGGTAATTGGATGCTTCTGGTTACGCTTTCGGTTGCGTCCATCAGTTTTACAGGACGGGTTACAGTGGAGCCATAAACGGCTCCATTTTCATTGTTATCatgttgagaaat    <  2:49336/246‑12 (MQ=255)
                                                                                                                                                                                                                                    aaaTTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTTCCAGCCTTGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGa  <  1:434512/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
GGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCT‑GTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGA  >  NC_000913/18241‑18718

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