Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 883,396 C→T intergenic (‑8/‑277) ybjG ← / → mdfA undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913883,3960CT94.4% 49.5 / NA 18intergenic (‑8/‑277)ybjG/mdfAundecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/10);  minor base G (0/1);  total (7/11)
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.93e-01

GTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTT  >  NC_000913/883171‑883634
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
gTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTGTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCGGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTGTCCAGCAGAACTGCCGGACGCAAATGGATTAAAAGGGAGcc                                                                                                                                                                                                                       <  2:113791/251‑1 (MQ=255)
  ccAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCa                                                                                                                                                                                                                      >  2:147126/1‑250 (MQ=255)
            aaCAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCa                                                                                                                                                                                                             >  1:79374/1‑249 (MQ=255)
                                                            tgtgCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCACCGGCACCACGGTAACCAAATCTTTAGCAATAAAAATCGCCACCGAGATCCTCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCtt                                                                                                                                                           >  1:29848/1‑251 (MQ=255)
                                                                  gTAAGCCCCCACAACCAAAGTACCACGGGCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTc                                                                                                                                                      <  1:213417/250‑1 (MQ=255)
                                                                   tAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTccc                                                                                                                                                    >  2:54234/1‑251 (MQ=255)
                                                                                   aaGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGttt                                                                                                                                     >  2:202141/1‑250 (MQ=255)
                                                                                           cGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgggt                                                                                                                            <  2:185775/251‑1 (MQ=255)
                                                                                                 gCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGt                                                                                                                      <  1:13868/251‑1 (MQ=255)
                                                                                                      ggCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTGCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTCCGGTGGGTGATTGTAACgg                                                                                                                 <  2:150027/251‑1 (MQ=255)
                                                                                                             gggTAGTCAAATCTTTAGCAATAAAAAACGCCAACGAGATCAGCCACGGATCCGAGTCTGGCGTCGTGTTAATAAGAGAGAATAGAGAGAGAATCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAATTCGTTTAGAATTATGTGCAGGGTTAACTGGGTAAAGCGGTATCTGCTTATTTCCCTCAATCGGTCTGTTTACGGTGGGTGTTTGTAACGGGCATatg                                                                                                          <  2:234609/250‑3 (MQ=255)
                                                                                                             cGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAgg                                                                                                          <  2:2842/251‑1 (MQ=255)
                                                                                                                             aGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGGTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCGTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTaa                                                                                           <  2:191498/250‑1 (MQ=255)
                                                                                                                                                                   aGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGc                                                    <  1:53072/251‑1 (MQ=255)
                                                                                                                                                                        ggCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATACCTTCCTGACGCAAAGTGATTAAAAGGGACCCAATACAGGCACGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAACCTTAACGAAAGCGTAGCTATACTCGTAATAATGTAAGACTGTGCTTAAc                                               >  1:134412/1‑251 (MQ=255)
                                                                                                                                                                                                  agagagaTTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAAtttgaat                     >  1:189532/1‑246 (MQ=255)
                                                                                                                                                                                                        aTTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAAtttgaatctctct               <  2:189532/251‑12 (MQ=255)
                                                                                                                                                                                                                     aGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCtttt  <  1:109021/251‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
GTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTT  >  NC_000913/883171‑883634

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