Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5710TA58.3% 1.3 / 11.2 12D398E (GAT→GAA) icdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base A (3/4);  total (8/4)
Fisher's exact test for biased strand distribution p-value = 8.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAA  >  minE/775503‑775621
                                                                    |                                                  
aTTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGatggaagg                                                  >  1:93258/1‑71 (MQ=255)
             tGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAAc                                          >  1:244424/1‑66 (MQ=255)
               gAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                                    <  1:100605/70‑1 (MQ=255)
               gAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                                    <  1:114820/70‑1 (MQ=255)
               gAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                                    >  1:130610/1‑70 (MQ=255)
               gAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                                    <  1:15682/70‑1 (MQ=255)
               gAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                                    <  1:192978/70‑1 (MQ=255)
                    cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                              >  1:106174/1‑71 (MQ=255)
                    cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                              >  1:187712/1‑71 (MQ=255)
                         tCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCagag                         >  1:45239/1‑71 (MQ=255)
                             cgcgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                     >  1:218082/1‑71 (MQ=255)
                                                           tCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGaaaa  >  1:32063/1‑60 (MQ=255)
                                                           tCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGaaaa  >  1:4884/1‑60 (MQ=255)
                                                               atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgc              >  1:87643/1‑44 (MQ=255)
                                                                    |                                                  
ATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAA  >  minE/775503‑775621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑