Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,157,782 T→C 73.7% L128S (TTG→TCG)  mdtB → multidrug efflux system, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,157,7820TC73.7% 65.2 / 23.8 38L128S (TTG→TCG) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base C (18/10);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 2.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAT  >  W3110S.gb/2157718‑2157845
                                                                |                                                               
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGa                                                            >  1:138795/1‑70 (MQ=255)
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAg                                                           >  1:189239/1‑71 (MQ=255)
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAg                                                           >  1:964183/1‑71 (MQ=255)
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAg                                                           >  1:267850/1‑71 (MQ=255)
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAg                                                           >  1:1366915/1‑71 (MQ=255)
aCATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAg                                                           >  1:1360317/1‑71 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTTCGACCAACTTGTCGCCGAGCg                                                         <  1:510697/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCg                                                         <  1:1493663/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCg                                                         <  1:1691909/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:87498/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:85002/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:585818/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:449927/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:255777/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:1800399/71‑1 (MQ=255)
  aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:1905762/71‑1 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:1340801/71‑1 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:1127813/71‑1 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        >  1:40415/1‑71 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:289827/71‑1 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:1360748/71‑1 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        >  1:1337460/1‑71 (MQ=255)
   ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:1111091/71‑1 (MQ=255)
                   ccGAGCATGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCg                                                         <  1:34815/54‑1 (MQ=255)
                            aaGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAg                                >  1:877308/1‑70 (MQ=255)
                            aaGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAg                                >  1:385033/1‑70 (MQ=255)
                            aaGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAg                                >  1:1665029/1‑70 (MQ=255)
                            aaGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAg                                >  1:579997/1‑70 (MQ=255)
                                 cAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGt                          >  1:5152/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:452105/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:1352487/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:1306180/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:150046/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:446254/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:113084/1‑71 (MQ=255)
                                                gCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATcc           >  1:1765156/1‑71 (MQ=255)
                                                      aCCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGAt     <  1:1392980/71‑1 (MQ=255)
                                                         aaCTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAt  <  1:1759713/71‑1 (MQ=255)
                                                                |                                                               
ACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAT  >  W3110S.gb/2157718‑2157845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: