Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,902,0910CA74.2% 40.8 / 23.7 31E652* (GAG→TAG) nrdDanaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base A (0/23);  total (0/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAG  >  minE/2902031‑2902161
                                                            |                                                                      
gggAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTaaa                                                               <  1:1184777/70‑1 (MQ=255)
gggAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTaaa                                                               <  1:1040359/70‑1 (MQ=255)
gggAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTaaa                                                               <  1:866258/70‑1 (MQ=255)
       gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACt                                                       <  1:923598/71‑1 (MQ=255)
       gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACt                                                       <  1:2403705/71‑1 (MQ=255)
       gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACt                                                       <  1:2178666/71‑1 (MQ=255)
       gTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACt                                                       <  1:2889021/71‑1 (MQ=255)
                          gggCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTaaa                                                               <  1:2697960/44‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCTGTGTATTGGTGCCGTAATACGGTa             <  1:1493867/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:330572/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:3016793/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:2990933/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:571469/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:2607427/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:2554600/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:696493/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:906986/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:2310164/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:245192/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1976972/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1827716/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1670196/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1589645/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1489949/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1422419/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1385964/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1287235/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1265914/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1195679/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1155546/70‑1 (MQ=255)
                                                  cACTCGAACTAACCGGTAAAGCCACACTCGTAGCACTCATAAATCGGTGTATTGGTGCCGTAATACGGTa             <  1:1505292/70‑1 (MQ=255)
                                                            cACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAg  >  1:2521465/1‑71 (MQ=255)
                                                            |                                                                      
GGGAGGCGTCATGGTTACCACATTTCGGGCAAGTGAAGCCTTTGCTGGTGCACTCGAACTCACCGGTAAAGCCACACTCGTAGCACTCATCAATCGGTGTATTGGTGCCGTAATACGGTACATGCTGATAG  >  minE/2902031‑2902161

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: