Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,155,197 G→A Q394Q (CAG→CAA)  mdtA → multidrug efflux pump membrane fusion protein MdtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,155,1970GA96.2% 72.4 / NA 26Q394Q (CAG→CAA) mdtAmultidrug efflux pump membrane fusion protein MdtA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/14);  minor base C (1/0);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGG  >  NC_000913/2154956‑2155428
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
tAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGAGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGGGGGGACAGACGGCATTGATCGCGTGACCGAAGGGGCGAAAGTGGACGGGGGGGAGGCCCACAGCGCCact                                                                                                                                                                                                                                 >  2:118288/1‑251 (MQ=255)
     gACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTcc                                                                                                                                                                                                                           >  1:122535/1‑251 (MQ=255)
         ccGAACAAACCGCCGTAGTGATCCCAACAGCCGACCTGCAAATGGGCAATGAAGGCCATGTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCTGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGTTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGaa                                                                                                                                                                                                                       <  2:37660/251‑1 (MQ=255)
         ccGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGGCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGAGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGCCGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCggc                                                                                                                                                                                                                        >  2:54931/1‑249 (MQ=255)
          cGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAag                                                                                                                                                                                                                      <  2:184299/251‑1 (MQ=255)
                 aaCGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGc                                                                                                                                                                                                               <  1:179651/251‑1 (MQ=255)
                                    cAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATAcgcg                                                                                                                                                                                            >  1:30681/1‑251 (MQ=255)
                                    cAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATAcgcg                                                                                                                                                                                            >  2:106657/1‑251 (MQ=255)
                                          cccTGCAACTGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGAGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCCTTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGAGACAGACGGCATTGATAGCCTGACCGAAGGGGCGACAGTGGAAGAGGTGGCAGCCCAAAGCGCCACT‑ACTACGGAAGAGACAGCCACCAGCCGCGAATCCGCGaaaac                                                                                                                                                                                       >  2:190285/1‑249 (MQ=255)
                                                               gCCATTTTGTCTGGGGGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAgg                                                                                                                                                                 <  2:123999/251‑1 (MQ=255)
                                                                           gggTGCTGAATAGCGAAAACAAGGTCAGCAAACAGCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTAccccc                                                                                                                                                       <  2:71383/249‑1 (MQ=255)
                                                                                      aGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggc                                                                                                                                          <  1:20493/251‑1 (MQ=255)
                                                                                           aaaCAAGGGAAGCAAACATCTGGTGACGACGGGAAGTCAGGACAGGCAGAAAGTGGTGATGAGTGCAGGTATTTCTGCGGGCGAGCGCGTGGTGACAGACGGCATTGAGCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACACCCGAGCAGCACAGGCGGccc                                                                                                                                     <  2:121981/251‑1 (MQ=255)
                                                                                               aaGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcg                                                                                                                                 >  1:84936/1‑251 (MQ=255)
                                                                                                aGGTCAGCAAACGTCTGGTGACGCAGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGGATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACA‑ACTCGGGAAGAGAAAGCCACCAGGCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTGGCCCCCGATGAGCACAGGCGGCCCGTcgc                                                                                                                                <  1:48889/251‑1 (MQ=255)
                                                                                                 ggTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCTAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcgcg                                                                                                                               <  2:3732/251‑1 (MQ=255)
                                                                                                  gTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCgcgc                                                                                                                              >  1:38507/1‑251 (MQ=255)
                                                                                                                                    ccAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGAGCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCGCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGTCCACCACgctgc                                                                                            <  2:144800/250‑1 (MQ=255)
                                                                                                                                                        cGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTc                                                                        <  2:116750/251‑1 (MQ=255)
                                                                                                                                                                                    tggtgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGc                                            >  1:221003/1‑251 (MQ=255)
                                                                                                                                                                                          cAGACGGCATTGATCGACTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCggtt                                      >  2:227301/1‑248 (MQ=255)
                                                                                                                                                                                                 cATTGAGCGCCTGACCGAAGTGGTGAAAGTGGTAGTGTTGGAAGCCCAACGTGCCCCT‑GCTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAATGGAGCTCGCACCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTAGTATGCGTCCTGTGGCCACCTCGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCgtctct                               <  2:184137/251‑6 (MQ=255)
                                                                                                                                                                                                     gatcgaatGACCAAAGGGGCGAACGTGCAAGTGGTGCAACCCCAAAGCGCCACT‑ACTCCGGACGAGAACGCCACCAGCCGCGACTACGCGAAAACAGGAGCACGCTCCTGATGCAGGTGCTACCCCCGCGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGAGTCCTGTGGCCACCACGCTGCTGCTGGTGGCGATCTTACTCGCCGGGATTAACGGTTATCTCGCCCCGCCCGGTTCGACGCTGCCg                           >  2:142311/8‑251 (MQ=255)
                                                                                                                                                                                                         gCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGa                         >  2:78911/1‑249 (MQ=255)
                                                                                                                                                                                                             gACCGAAGGGGCGAAAGGGGAAGGGGTGGAAGCCCAAAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTAGCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTgg                    <  2:204577/250‑1 (MQ=255)
                                                                                                                                                                                                                              gggaagtggTGGAAGCCCAAAGCGCCACTCAC‑CCGGAAGAAAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCTGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAATTGGACTATCCGACCATTCAgg  <  2:108925/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
TAGTCGACACCGAACAAAACGCCGTAGTGATCCCAACAGCCGCCCTGCAAATGGGCAATGAAGGCCATTTTGTCTGGGTGCTGAATAGCGAAAACAAGGTCAGCAAACATCTGGTGACGCCGGGCATTCAGGACAGTCAGAAAGTGGTGATCCGTGCAGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACT‑ACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGCCGGAAGTGGACTATCCGACCATTCAGG  >  NC_000913/2154956‑2155428

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: ‑