Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,886,1630AC54.2% 15.7 / 24.3 24intergenic (‑288/+71)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base A (7/4);  new base C (13/0);  total (20/4)
Fisher's exact test for biased strand distribution p-value = 3.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGT  >  W3110S.gb/2886106‑2886223
                                                         |                                                            
gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAgg                                                            <  1:4538605/60‑1 (MQ=255)
gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAgg                                                            <  1:2523744/60‑1 (MQ=255)
    ccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAAtgct                                                      >  1:98435/1‑62 (MQ=255)
       aCGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATgctgct                                                   >  1:172757/1‑62 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAAtgctgat                                                   >  1:1797400/2‑59 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:3156590/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:4590883/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:4149187/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:1062876/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:2195486/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:2109773/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:2048109/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:1924695/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:1648787/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:1133035/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCTGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:1103721/2‑61 (MQ=255)
        tgCCATTCTATACCTGTTAATTCTCCCAGCGACGTTCTCGTTTAAATGCCGGAATgctgct                                                   >  1:3966204/2‑61 (MQ=255)
                   aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa                                       >  1:3485424/1‑62 (MQ=255)
                   aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa                                       >  1:4529951/1‑62 (MQ=255)
                   aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaaa                                       >  1:725551/1‑62 (MQ=255)
                                  ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaa                                        >  1:3419283/1‑46 (MQ=255)
                                  ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATaaa                                        >  1:1549579/1‑46 (MQ=255)
                                   cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCAtt                         <  1:907924/60‑1 (MQ=255)
                                                        cAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGt  <  1:78443/62‑1 (MQ=255)
                                                         |                                                            
GCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGT  >  W3110S.gb/2886106‑2886223

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: