Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,654,5071.T60.0% 56.8 / 43.4 30pseudogene (46/111 nt)ynfPECK1574:JW5257:b4534; hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/10);  new base T (18/0);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 2.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTG  >  W3110S.gb/1654461‑1654548
                                               |                                         
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTa                           >  1:1118367/1‑64 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTa                           >  1:1380398/1‑64 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTa                           >  1:2220109/1‑64 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTa                           >  1:1598275/1‑64 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:999835/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:983572/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:666121/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:528344/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:451438/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:327954/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:2211039/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:2079227/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:1151361/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:1757547/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:1213570/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:178083/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTTCACCAAATCCTTAg                          >  1:1263174/1‑65 (MQ=255)
aaaTCCATAAATTCAAGCGCAGTGCCCAGCCATCCAGATACTGCTGCTTTTCACCAAATCCTTa                           >  1:1874171/1‑64 (MQ=255)
         aaTTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTc                >  1:1954461/1‑65 (MQ=255)
         aaTTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTc                >  1:2265242/1‑65 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:1363723/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:635803/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:562389/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:1315664/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:253089/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:248535/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:1771441/65‑1 (MQ=255)
                       gCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCATTCGTGTTTTTCTATTg  <  1:1489582/65‑1 (MQ=255)
                        cccAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:1090513/64‑1 (MQ=255)
                                   cGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTg  <  1:357381/53‑1 (MQ=255)
                                               |                                         
AAATCCATAAATTCAAGCGCAGTGCCCAGCCATCCCGATACTGCTGCTTTCACCAAATCCTTAGTGCTTCTTTCGTGTTTTTCTATTG  >  W3110S.gb/1654461‑1654548

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: