Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,157,782 T→C 14.8% 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,7820TC14.8% 62.9 / 4.8 27L128S (TTG→TCG) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base C (2/2);  total (13/14)
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.50e-01

TTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAT  >  W3110S.gb/2157721‑2157845
                                                             |                                                               
ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:3244689/71‑1 (MQ=255)
ttACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGa                                                        <  1:2336729/71‑1 (MQ=255)
             tCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAAccc                                            <  1:1310226/70‑1 (MQ=255)
             tCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAAccc                                            <  1:546452/70‑1 (MQ=255)
             tCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAAccc                                            <  1:544209/70‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:902213/63‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:1025896/63‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:573632/63‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:2764572/63‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:2675085/63‑1 (MQ=255)
                           gTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGtt                                     <  1:1618973/63‑1 (MQ=255)
                             gCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAg                           >  1:2003414/1‑71 (MQ=255)
                             gCAGGCCGCGATTAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAg                           >  1:463751/1‑71 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:2218677/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:624931/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:1614971/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:477613/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:163075/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:2834388/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:2800156/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:203289/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                               >  1:2222595/1‑40 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAg                           <  1:1242450/60‑1 (MQ=255)
                                        ttAACGCTGCGACCAACTTGTCGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAg                           <  1:2200189/60‑1 (MQ=255)
                                                    ccAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGa              >  1:1241494/1‑61 (MQ=255)
                                                    ccAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGAt     >  1:1929306/1‑70 (MQ=255)
                                                      aaCTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAt  <  1:3245414/71‑1 (MQ=255)
                                                             |                                                               
TTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCAT  >  W3110S.gb/2157721‑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: