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,4770TC90.9% 48.0 / ‑3.1 22V330A (GTG→GCG) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (4/16);  total (4/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTT  >  NC_000913/18230‑18692
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                          
gcgaagcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggt                                                                                                                                                                                                                         <  1:232294/251‑1 (MQ=255)
 ggaagcgaCTGGTGCATTTGTTACCCCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCGCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGGCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTAtggcggtg                                                                                                                                                                                                                        <  2:129482/250‑1 (MQ=255)
    agcgaCTGGAGCATGTGTTACACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTgggg                                                                                                                                                                                                                     <  1:137060/251‑1 (MQ=255)
       gaCTTGAGCATGAGTTACACCCGTGGGTGGCGTATCTGATTTAGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGGGTCACGATGGCTGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGACTTGCTGCTTGGCAAACCGCTGGGGATTAGTCTGTTCTGATGGTTGGCACTGCGATTGAAACAGGCGCATGTGCCTGGGGGAACGACTTATCAGCATATTATGGCGGAggggttt                                                                                                                                                                                                                  >  2:285834/1‑248 (MQ=255)
            gagcttgtgttttgccagggggtggCTTATCTTGTTTTGCCGCTAGTTGCAATTGCTAATGCTTGCGTTTCACGTGAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGTCATTGTGGATCATCGCGGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTGTGTGCTGCGGGTTTGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTAGCAGCAAATTATGGCGGTGGGGATCCtgt                                                                                                                                                                                                              <  2:32725/231‑1 (MQ=255)
                       acACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCtg                                                                                                                                                                                                               <  1:379706/237‑1 (MQ=255)
                                                             cATTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGt                                                                                                                                                           <  1:285834/251‑1 (MQ=255)
                                                              aTTTGCTAATGCTGGCGTTTCACTGCAAGGTGTCACGGTGGATGGCTTGACCTCCATGCTGCCATTGAGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTa                                                                                                                                                          <  1:215110/251‑1 (MQ=255)
                                                                ttGCTAATGCTGGCGTTTCACGGCAAGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGTGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAg                                                                                                                                                         <  1:58575/250‑1 (MQ=255)
                                                                                         aaGGTGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTgg                                                                                                                               <  1:2817/251‑1 (MQ=255)
                                                                                                cACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGGGCCTGTGAGGGATCGGGTTTTTTATGTCTCTCTTTTTTTTCAGCCTGGGCTTTGTTGGCGTAGagaccggcatgattaactgtgggaca                                                                                                                          >  1:264370/1‑225 (MQ=255)
                                                                                                      ggATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTa                                                                                                                  <  1:337376/251‑1 (MQ=255)
                                                                                                               gACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCGGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTAGGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTc                                                                                                          <  1:218220/250‑1 (MQ=255)
                                                                                                                          ggccttgttgggtcatcGCGGTCTTGCGGATTGGCAAGCCTCTGGTGATTAGTCTGTTCTGCTGGTTGTCGCTGCGTTGGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAAGTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAGCTCGGTATCCTGGTCGGTTCTATctct                                                                                              <  2:86047/239‑1 (MQ=255)
                                                                                                                                                            cAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGGACCTGTGCGGGATCGGGTTTTATTTTTCTTTCTTTTTTTCCCGCCCGGCCTTTTGTTGCGTTAATCCAGAAATGATTAAATGTGGGGAAATCGGGATCCTGGTCGGTTTTTTCTCTTTGGCGGGTAATTggtaaaaatgggtatgcgttt                                                            >  1:349309/1‑232 (MQ=255)
                                                                                                                                                                ccGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGGGGGGGGCCTCTTCGGGTTCGGTTTTAATTTGTATTTCTTTATTGCCAGCCCGGCCTTTGGTTGCGTGGATCCAGAACCGCTTTAATGGGCGGAACCCGGGATCCTGGGCGGTTCTCTCTCTTATGCGGTGATTAGATACaacctgtttcacgttctctt                                                       >  1:137106/1‑232 (MQ=255)
                                                                                                                                                                     ggggatgtttctgtgtttctgGTTGTCGCTGCGTTGTAAATGTGCGCATCTGCCTGAGGGACCGTCTTGGCAGCAAATTAGGGCGGTGGTGATCCTGTGCGGTATCGGTGTTTCTATTGCTATCTTTATGGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGTGGTAATGGGATACAGCTGGTTACGCGTTCGTTGGCGTc                                                  <  2:351223/233‑1 (MQ=255)
                                                                                                                                                                                 gTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTg                                        <  1:163086/249‑1 (MQ=255)
                                                                                                                                                                                 gTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGAc                                      <  1:167775/251‑1 (MQ=255)
                                                                                                                                                                                      gctgGTTTGCGCGTCGTTTGTAACGGGCGCATCTGCCTTAGGGAACGACTTATCATCACATTATGGCGGTGGGGAGTCGGTGCGGTATCGGTTTTACTATGTCTATTTTTATTGCCGTCCTGGCCGTTGGTAGCGTAGATCCAGAACTGATGAACTGGGCGATATTCGGTATCCTGGTCGGTTCTATCTTTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCACCATTTTGACAgg                                   <  2:65068/249‑1 (MQ=255)
                                                                                                                                                                                                     gTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGc                  <  1:254433/251‑1 (MQ=255)
                                                                                                                                                                                                               ggCGCATCTGCCTGAGGGAACGTCTTATCAGCAAATTATGGCGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGc        <  1:131864/251‑1 (MQ=255)
                                                                                                                                                                                                                     atttGCCTGAGGGAACGACTTATCAGCAAAGTATGGTGGTGGGGATACTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACGCGGTATCCTGGTCGGTTCTATCTCTGCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCt   <  1:164385/248‑1 (MQ=255)
                                                                                                                                                                                                                         gccggAGGGAACGACTTATCAGCAAATTATGGCGGTGGGGATCCTGTTCTGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCtt  <  1:14414/244‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                          
GCGAAGCGACTGGAGCATGTGTTGCACCCGTGGGTGGCGTATCTGATTTTGCCGCTGTTTGCATTTGCTAATGCTGGCGTTTCACTGCAAGGCGTCACGCTGGATGGCTTGACCTCCATTCTGCCATTGGGGATCATCGCTGGCTTGCTGATTGGCAAACCGCTGGGGATTAGTCTGTTCTGCTGGTTGGCGCTGCGTTTGAAACTGGCGCATCTGCCTGAGGGAACGACTTATCAGCAAATTATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTT  >  NC_000913/18230‑18692

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: