Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 366,370 Δ1 bp coding (365/1083 nt) lacI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb366,3700A.92.0% 86.2 / 5.1 25coding (365/1083 nt)lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (0/23);  total (2/23)
Fisher's exact test for biased strand distribution p-value = 3.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGAT  >  W3110S.gb/366330‑366397
                                        |                           
aGCAATGGCATCCTGGTCATCCATCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1195666/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:797847/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1154175/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:786859/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:725407/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:541957/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:412569/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:348769/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:307339/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:266229/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:189170/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1637967/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1632483/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:156247/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1549077/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1524971/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1467668/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1350877/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1332521/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1320921/62‑1 (MQ=255)
aGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1285605/62‑1 (MQ=255)
      ggCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGAt  >  1:1428276/1‑62 (MQ=255)
      ggCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGAt  >  1:601110/1‑62 (MQ=255)
           ccTGGTCATCCAGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:289009/51‑1 (MQ=255)
                      aGCGGATAGTTAATGATCGCCCACTGACGCGTTGCGCgag       <  1:1270620/40‑1 (MQ=39)
                                        |                           
AGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGAT  >  W3110S.gb/366330‑366397

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