Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 300,220 G→A 100% L372L (CTG→CTA mdlB → fused predicted multidrug transporter subunits and ATP binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE300,2200GA90.4% 158.2 / 6.6 52L372L (CTG→CTAmdlBfused predicted multidrug transporter subunits and ATP binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base A (3/44);  total (3/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGG  >  minE/300164‑300270
                                                        |                                                  
aaTCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTAGTCGGGCAt                                           >  1:2139452/1‑66 (MQ=255)
aaTCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTAGTCGGGCAt                                           >  1:434404/1‑66 (MQ=255)
aaTCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTAGTCGGGCATa                                          >  1:1308628/1‑67 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGTCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1684018/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:68637/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:2098875/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:2142873/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:277342/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:370576/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:397087/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:477579/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:526799/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:572088/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1042930/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:687395/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:702602/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:705182/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:751769/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:764243/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:788584/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:793246/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:833328/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:897446/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:992001/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1408668/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1001810/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1066149/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1163889/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1208476/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1244466/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1261087/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1262859/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1317373/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1350313/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1401389/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:2093737/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1469723/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1517227/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1524821/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1529337/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1585528/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1590128/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1849661/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1900365/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:1911468/66‑1 (MQ=255)
                                 ccTTCGCGCAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:2008824/66‑1 (MQ=255)
                                    tctcgcAATTTTGTGGCGCTAGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTa          <  1:343729/60‑1 (MQ=255)
                                        gcAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg  <  1:621648/67‑1 (MQ=255)
                                        gcAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg  <  1:699238/67‑1 (MQ=255)
                                        gcAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg  <  1:2116128/67‑1 (MQ=255)
                                        gcAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg  <  1:945356/67‑1 (MQ=255)
                                         cAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATggg  <  1:976551/66‑1 (MQ=255)
                                                        |                                                  
AATCTGGTGCTAAAGAACATTAATCTCTCTGTGCCTTCGCGCAATTTTGTGGCGCTGGTCGGGCATACCGGCAGTGGCAAAAGCACCCTCGCCAGTTTATTGATGGG  >  minE/300164‑300270

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