Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,763,6360AG21.9% 57.9 / 14.4 32noncoding (208/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (4/3);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  minE/1763571‑1763705
                                                                 |                                                                     
tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                                  >  1:1440639/1‑71 (MQ=255)
 ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc                                                                 >  1:1694180/1‑71 (MQ=255)
  gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGATGCCCGAGGGTccc                                                                 <  1:1380256/70‑1 (MQ=255)
     aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                             >  1:256722/1‑71 (MQ=255)
       ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                           >  1:227598/1‑71 (MQ=255)
        cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                          <  1:752168/71‑1 (MQ=39)
        cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt                                                           <  1:475515/70‑1 (MQ=40)
           tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                       >  1:1979190/1‑71 (MQ=39)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                      <  1:400846/70‑1 (MQ=255)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:2188607/71‑1 (MQ=39)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     >  1:1110209/1‑71 (MQ=39)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     >  1:1311223/1‑71 (MQ=39)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:1430830/71‑1 (MQ=39)
             ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                     <  1:1427045/71‑1 (MQ=39)
                caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                                  >  1:1822044/1‑71 (MQ=39)
                 accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGGCCCCCTCTTTGGTCTTGCg                                                 <  1:1754755/71‑1 (MQ=255)
                        aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                          >  1:680475/1‑71 (MQ=34)
                                     tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                    >  1:1912397/1‑64 (MQ=35)
                                     tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                    >  1:255925/1‑64 (MQ=35)
                                      gggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTAc                            >  1:2492912/1‑71 (MQ=34)
                                       ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCt                              <  1:2152553/68‑1 (MQ=35)
                                       ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                             <  1:690454/69‑1 (MQ=35)
                                            aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                      <  1:2367888/71‑1 (MQ=34)
                                                     cAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa             <  1:2680425/71‑1 (MQ=34)
                                                      aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCg                                       >  1:2854034/1‑44 (MQ=38)
                                                        gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                        <  1:1115698/41‑1 (MQ=38)
                                                           gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                             >  1:659568/1‑49 (MQ=38)
                                                           gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                             >  1:213967/1‑49 (MQ=38)
                                                           gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                   <  1:1614129/59‑1 (MQ=35)
                                                             ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCACTAGTTATCCCCCTc     <  1:2976372/71‑1 (MQ=255)
                                                              cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTcc    <  1:2542788/71‑1 (MQ=34)
                                                               gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGTGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCa   >  1:1942451/1‑71 (MQ=21)
                                                                aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt  >  1:576395/1‑71 (MQ=34)
                                                                 |                                                                     
TGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAT  >  minE/1763571‑1763705

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: