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

CTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAT  >  NC_000913/1287086‑1287343
                                                                                                                                  |                                                                                                                                
cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGtt                                                                                                                         >  2:832667/1‑139 (MQ=11)
           gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACtttt                                                                                                              >  2:427213/1‑139 (MQ=255)
             gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGa                                                                                                            <  1:427213/139‑1 (MQ=11)
               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGaaa                                                                                                          >  2:447272/1‑139 (MQ=11)
                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                      <  1:832667/139‑1 (MQ=255)
                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                      <  1:996973/139‑1 (MQ=255)
                          ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                               <  1:671698/139‑1 (MQ=255)
                                      cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                                                                >  1:104672/1‑126 (MQ=11)
                                      cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                                                                <  2:104672/126‑1 (MQ=255)
                                                             tCGAGCCGAAACGTGATCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAATATTTGATTGCGGATTCG‑CTAGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                                            <  1:576209/139‑1 (MQ=14)
                                                                      aCCTTAATCGAAGCTTCTCATCCTTCCCCGAATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                                                   <  2:901839/139‑1 (MQ=14)
                                                                                                   gCATGGGCAGTATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACTGTCTGCTCCCTTTGGCCGCTCGGGAAc                       <  1:695888/138‑1 (MQ=18)
                                                                                                   gCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcc                      <  2:1046083/139‑1 (MQ=21)
                                                                                                   gCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                       >  2:695888/1‑138 (MQ=18)
                                                                                                       gggCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                             >  1:203075/1‑112 (MQ=17)
                                                                                                       gggCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                             <  2:203075/112‑1 (MQ=17)
                                                                                                        ggCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          <  2:100285/66‑1 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          >  1:100285/1‑66 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          >  2:176679/1‑66 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          <  1:176679/66‑1 (MQ=255)
                                                                                                         gCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgtt                                                                                               >  2:1021993/1‑58 (MQ=255)
                                                                                                         gCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                               <  1:1021993/60‑1 (MQ=255)
                                                                                                                       tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTg                                                                                                       >  1:866391/1‑38 (MQ=255)
                                                                                                                       tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTg                                                                                                       <  2:866391/38‑1 (MQ=255)
                                                                                                                       tGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAt  >  2:979601/1‑139 (MQ=18)
                                                                                                                         cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg                                                         >  1:245389/1‑82 (MQ=14)
                                                                                                                         cGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg                                                         <  2:245389/82‑1 (MQ=17)
                                                                                                                         cGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                                                   >  2:845337/1‑88 (MQ=18)
                                                                                                                                  |                                                                                                                                
CTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAT  >  NC_000913/1287086‑1287343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑