Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,132,6111.G33.3% 17.9 / 37.0 24coding (337/1890 nt)deaDATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base . (8/8);  new base G (0/8);  total (8/16)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCATTGCTTCAGCAACC  >  minE/2132561‑2132681
                                                   |                                                                        
aaCGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAAcc                                                     >  1:3467221/1‑70 (MQ=255)
aaCGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACca                                                    >  1:3218671/1‑71 (MQ=255)
aaCGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACca                                                    >  1:1773043/1‑71 (MQ=255)
aaCGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACca                                                    >  1:428098/1‑71 (MQ=255)
aaCGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACca                                                    >  1:2524998/1‑71 (MQ=255)
aaCGAACTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACca                                                    >  1:3728556/1‑71 (MQ=255)
                    gCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATgtgt                                >  1:2812753/1‑71 (MQ=255)
                    gCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCCCCGTACAGAGCAACCACATTTACGCCGCGCATgtgt                                >  1:1463690/1‑71 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:1003428/67‑1 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:3443032/67‑1 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:3189821/67‑1 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:2418061/67‑1 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:62038/67‑1 (MQ=255)
                       gggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:1654163/67‑1 (MQ=255)
                        ggCGCGTAATTGCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:1490038/66‑1 (MQ=255)
                         tcgcgTAATTTCACGTCATAACGCTGGCCGCCGCCGTACAGAGCAACCACATTTACGCCGCGCAt                                    <  1:4342503/64‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:634550/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:635546/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:2551933/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:2451516/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:1702978/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:1491678/71‑1 (MQ=255)
                                     aCGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:1268574/71‑1 (MQ=255)
                                      cGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCAt               <  1:4078052/70‑1 (MQ=255)
                                                  g‑‑‑gCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCATTGCTTCAGCAAcc  >  1:3498279/1‑71 (MQ=255)
                                                   |                                                                        
AACGATCTGCGGCCCCTGACGCAGGGCGCGTAATTGCACGTCATAACGCTG‑‑‑GCCGCCGTACAGAGCAACCACATTTACGCCGCGCATGTGTTTAGAGAAATCCGTCATTGCTTCAGCAACC  >  minE/2132561‑2132681

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: