Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349432,981,2530CT46.7% ‑6.0 / 18.4 15noncoding (1267/1550 nt)EO946_RS1543516S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base T (2/5);  total (2/13)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAG  >  NZ_CP034943/2981117‑2981355
                                                                                                                                        |                                                                                                      
tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCt                                                                                                     <  1:3044733/137‑1 (MQ=1)
tgaaCAAGGCCCGGGAACGTATTCAAAGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCt                                                                                                     >  2:3164302/4‑140 (MQ=0)
          cctggAACTTATTCACCGCGGCCTGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGtt                                                                                           <  1:2090034/137‑1 (MQ=0)
          cctggAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGtt                                                                                           <  2:945004/137‑1 (MQ=1)
                   taatCACCGCGGCATGCTGATCAGCTATTAAAATCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCAtt                                                                                  >  2:1506017/4‑140 (MQ=1)
                     gtccccGCGGCATGCAGATCCGCGAGTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGt                                                                                <  1:3246228/136‑1 (MQ=0)
                         acgcgtgaTGCTGCTCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGATTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCa                                                                            <  2:304078/133‑1 (MQ=0)
                               caggCTGATCCGCGATTACTCGCGATTCCAGCTTCCCGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACgtgtg                                                                      <  1:2265630/137‑1 (MQ=0)
                               caagCTGATCCGCGATTACGAGCGATTCCAGCTTCACGCACTCGACTTGCAGACTGCGCTCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACgtgtg                                                                      <  1:643587/137‑1 (MQ=1)
                                     gagccGCGATTACTCGCGAATCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGccc                                                                <  2:2189471/137‑1 (MQ=0)
                                                      gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACTTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                               <  1:1651926/137‑1 (MQ=1)
                                                      gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGAGTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                               <  2:2652518/137‑1 (MQ=1)
                                                                 caagCAGTCGAGTTGCAGAATGCTATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCACTTTGTTCTGTCGATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATTATTTGACg                                    <  1:2792465/137‑1 (MQ=1)
                                                                                                   cgaccAGACTTGCGGGATTGGCTAAACCTCGCGGTTTCGCGGCCCTTTGTTCTGTCCATTGGAGCACGTGTGTAGCCCAGGTCATACGGGGCATGATGATTTGACGTCATCCCCACCGTCCTCCGGTTTGTCACCGGCAg  <  2:2979382/136‑1 (MQ=1)
                                                                                                   agcaCAGATTTGTGGGAGTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg  <  1:1052160/137‑1 (MQ=0)
                                                                                                                                        |                                                                                                      
TGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAG  >  NZ_CP034943/2981117‑2981355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: