Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,020,867 A→G 100% S159P (TCC→CCC)  sulA ← SOS cell division inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,020,8670AG76.5% 29.0 / 5.6 17S159P (TCC→CCC) sulASOS cell division inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base G (8/5);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAG  >  W3110S.gb/1020798‑1020923
                                                                     |                                                        
aCAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGaa                                                         <  1:1903148/71‑1 (MQ=255)
aCAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGaa                                                         <  1:2282256/71‑1 (MQ=255)
  aaaaaaaGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAg                                                       <  1:3041643/71‑1 (MQ=255)
      aaaGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGt                                                   <  1:2551606/71‑1 (MQ=255)
             cAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTg                                             <  1:2240448/70‑1 (MQ=255)
                          aaaTTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAgg                                   >  1:2231479/1‑67 (MQ=255)
                               tACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCtt                           <  1:891058/70‑1 (MQ=255)
                                        aTACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCa                 <  1:2393545/71‑1 (MQ=255)
                                          aCAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCat                >  1:3070548/1‑70 (MQ=255)
                                           cAAATTAGAGTGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGa                                  >  1:2080206/1‑51 (MQ=255)
                                                     tGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata            <  1:145375/63‑1 (MQ=255)
                                                     tGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata            <  1:1933156/63‑1 (MQ=255)
                                                     tGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCataata            <  1:3626713/63‑1 (MQ=255)
                                                     tGAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCa     >  1:213480/1‑70 (MQ=255)
                                                      gAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGa                     >  1:908715/1‑53 (MQ=255)
                                                      gAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGa                     >  1:1101261/1‑53 (MQ=255)
                                                      gAATTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGa                     >  1:3679987/1‑53 (MQ=255)
                                                       aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGagag                                    <  1:1391928/37‑1 (MQ=255)
                                                       aaTTTTTAGCCCGGGAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAg  >  1:1902048/1‑71 (MQ=255)
                                                                     |                                                        
ACAAAAAAAGTTCCAGGATTAATCCTAAATTTACTTAATGATACAAATTAGAGTGAATTTTTAGCCCGGAAAGTTGTCTCGTGGCGTGAGAGGATGCGCTTACCGGACGCATAATAAACCCCATAG  >  W3110S.gb/1020798‑1020923

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: