Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,481,608 A→G L190S (TTA→TCA)  emrY ← tripartite efflux pump membrane subunit EmrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,481,6080AG95.5% 65.9 / ‑5.3 22L190S (TTA→TCA) emrYtripartite efflux pump membrane subunit EmrY
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (10/11);  total (10/12)
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.33e-01

ACGGGACTTAAACAAACTGAGATCAAGAATCGGGTTCTCTGAGGTCGACTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTA  >  NC_000913/2481361‑2481841
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
aCGGGACTTAAACAAACTGAGATCAAGAATCGGGTTCTCTGAGGTCGACTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGc                                                                                                                                                                                                                                        >  1:409552/1‑251 (MQ=255)
  gggACTTAAACAAACTGAGATCAAGAATCGGGTTCTCTGAGGTCGACTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCAt                                                                                                                                                                                                                                      >  1:173840/1‑251 (MQ=255)
             aaacGGAGATCAAGAAGCGGGTTCTCTGAGGTCGACTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGTCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTGCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGAc                                                                                                                                                                                                                           <  2:89782/251‑1 (MQ=255)
                                               aCTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGACCTTTATCAAGCATAATTTGCAAGCCACCAACCCCGAGCCCTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATTTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATc                                                                                                                                                                                         <  1:147960/251‑1 (MQ=255)
                                                  cccAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGTACATTTATTAAAAATATcc                                                                                                                                                                                        <  2:369645/249‑1 (MQ=255)
                                                      aaTGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAAcccc                                                                                                                                                                                  >  2:476934/1‑251 (MQ=255)
                                                          aCTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCt                                                                                                                                                                              >  1:188693/1‑251 (MQ=255)
                                                              aaGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAg                                                                                                                                                                          <  1:448640/251‑1 (MQ=255)
                                                                     cAGAGAAATAACTGAGACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATca                                                                                                                                                                    <  2:365450/250‑1 (MQ=255)
                                                                      agagAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATcaca                                                                                                                                                                  <  1:476934/251‑1 (MQ=255)
                                                                          aaaTAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTGTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTAGCACAAata                                                                                                                                                              <  2:358476/251‑1 (MQ=255)
                                                                                                              tCGAGTTGAAACCATCAAGAACGCGACCTCTATCGAGCATAATTTGCAAGACAACAACCCCCAGCACTAACAGGGTCAGTCCTGGAAGATTCATTTTGACCGGTGAAGTATCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATACCCATCGGGACATTGATTAAAAATATCCAACCCCAGCTACAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATAGGAGCGAt                                                                                                                           >  2:527877/1‑250 (MQ=255)
                                                                                                                                    gcgcCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCcac                                                                                                                                        >  1:509394/1‑215 (MQ=255)
                                                                                                                                    gcgcCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCcac                                                                                                                                        <  2:509394/215‑1 (MQ=255)
                                                                                                                                                                  accaACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAatat                                                                                                                                                                                          >  1:186917/1‑135 (MQ=255)
                                                                                                                                                                  accaACACCGAGCACGAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAatat                                                                                                                                                                                          <  2:186917/135‑1 (MQ=255)
                                                                                                                                                                             gCACTAACAGGGTAAGTCCTGGGAGATTCATATTGACCGGTAAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATACAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGAACATAATGCCAGAGCAAATGTTCTTTTTTCTGGAGGATACTTCCTtaat                                                           >  2:311480/1‑251 (MQ=255)
                                                                                                                                                                                       ggTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAActct                                                 <  2:181400/251‑1 (MQ=255)
                                                                                                                                                                                        gTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCtg                                                <  2:132358/251‑1 (MQ=255)
                                                                                                                                                                                                               gACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAATGTTGAGCATAATGGCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGGGACAGTGGAATTAACGGCCCCg                         <  2:472769/251‑1 (MQ=255)
                                                                                                                                                                                                                        aGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAAcc                >  2:37800/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tCCTTTAAGTAAGGTTGAGCATAATGTCAGGACGATAATCCCCATACGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTa  >  1:88233/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
ACGGGACTTAAACAAACTGAGATCAAGAATCGGGTTCTCTGAGGTCGACTCCCAAATGACTAAAGAGATCAGAGAAATAACTGATACTACTGTTAATATTATTATTGTACTCGAGTTGAACCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTA  >  NC_000913/2481361‑2481841

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