Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,2150TC39.3% 24.0 / 25.3 28noncoding (22/171 nt)rttRrtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (5/12);  new base C (7/4);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT  >  NC_000913/1287091‑1287345
                                                                                                                             |                                                                                                                                  
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                         >  2:714683/1‑136 (MQ=255)
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                         <  1:714683/136‑1 (MQ=255)
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                                                                      <  2:3959/139‑1 (MQ=255)
                                                            gccGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGa                                                          >  2:709237/3‑139 (MQ=11)
                                                                aaccttaatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                                                                                             >  1:273367/9‑84 (MQ=255)
                                                                aaCCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                                                                                             <  2:273367/84‑1 (MQ=255)
                                                                      aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGAACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCt                                                >  1:479332/3‑139 (MQ=14)
                                                                       atCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                               <  2:479332/138‑1 (MQ=14)
                                                                          gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcg                                            <  2:101688‑M2/139‑63 (MQ=255)
                                                                                      ccTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctct                                <  2:394493‑M2/139‑75 (MQ=255)
                                                                                         tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCgg                             >  1:29800/1‑139 (MQ=18)
                                                                                          ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCggg                            >  2:120446/1‑139 (MQ=18)
                                                                                          ccccGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgct                            <  1:100578‑M2/139‑79 (MQ=255)
                                                                                          ccccGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgct                            <  1:639048‑M2/139‑79 (MQ=255)
                                                                                             cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAGCCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                         <  1:120446/139‑1 (MQ=18)
                                                                                             cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                         <  2:566558/139‑1 (MQ=21)
                                                                                                   ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGAc                                                                >  2:237293/1‑94 (MQ=11)
                                                                                                   ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGAc                                                                <  1:237293/94‑1 (MQ=11)
                                                                                                    gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCccacca                  >  2:586083/1‑139 (MQ=18)
                                                                                                    gCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgag                  <  1:391249‑M2/139‑89 (MQ=255)
                                                                                                             ttgaTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAAtt         >  2:461199/4‑139 (MQ=17)
                                                                                                             ttgaTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAAtt         >  2:697032/4‑139 (MQ=18)
                                                                                                                  tGCGGATTCGTTGGGAAGTTTAGGGACTTTTGAAAGTgatggtggtgggggaaggattcctcagttctgcgcgcttcgcccttcgggtagttacctgcgacaacgctctctcgctggcgctcgagccgaaccttaatcg    <  2:568327‑M2/139‑103 (MQ=255)
                                                                                                                  tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcg    <  1:676214‑M2/139‑103 (MQ=255)
                                                                                                                    cGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcgaa  >  2:391249‑M2/1‑35 (MQ=255)
                                                                                                                    cGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcgttgcctgcggcaacgctctctcgctggcgctcgagccgaaccttaatcgaa  <  1:620776‑M2/139‑105 (MQ=255)
                                                                                                                         tCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcg                                                     <  1:684154‑M2/83‑54 (MQ=255)
                                                                                                                         tCGTTGGGAAGTTCAGGGACTTTTGAAAGTgatggtggtgggggaaggattactcagcgctgcgcgcttcgcccttcgggtcg                                                     >  2:684154‑M2/1‑30 (MQ=255)
                                                                                                                             |                                                                                                                                  
CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT  >  NC_000913/1287091‑1287345

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