Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,856,8780T.62.5% 25.3 / 51.0 32intergenic (‑177/+182)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base T (12/0);  new base . (0/20);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 4.43e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCC  >  minE/1856830‑1856919
                                                |                                         
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                   >  1:64861/1‑57 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                >  1:528179/1‑60 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:904499/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:610588/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:739746/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:78934/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:852934/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:173406/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:889250/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                              >  1:146612/1‑62 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                               >  1:560321/1‑61 (MQ=255)
gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                               >  1:1170720/1‑61 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:871258/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:883669/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:752794/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:728505/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:727056/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:71306/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:659533/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:580591/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:59354/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:1021218/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:505340/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:46503/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:452054/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:406386/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:388924/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:342770/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:227341/62‑1 (MQ=255)
                           cAGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:169259/62‑1 (MQ=255)
                            aGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:558650/61‑1 (MQ=255)
                            aGGTTGATAGAAATATCGCCGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc  <  1:105084/61‑1 (MQ=255)
                                                |                                         
GAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCC  >  minE/1856830‑1856919

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: