Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,2190GA35.7% 27.7 / 17.6 28noncoding (18/171 nt)rttRrtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base G (9/9);  new base A (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.

CTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAT  >  NC_000913/1287086‑1287343
                                                                                                                                      |                                                                                                                            
cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtt                                                                                                                         >  2:832667/1‑139 (MQ=11)
           gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                                                                                              >  2:427213/1‑139 (MQ=255)
             gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                                                                                            <  1:427213/139‑1 (MQ=11)
               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                                                                                          >  2:447272/1‑139 (MQ=11)
                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                      <  1:832667/139‑1 (MQ=255)
                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                      <  1:996973/139‑1 (MQ=255)
                          ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                               <  1:671698/139‑1 (MQ=255)
                                      cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                                                                >  1:104672/1‑126 (MQ=11)
                                      cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg                                                                                                <  2:104672/126‑1 (MQ=255)
                                                             tCGAGCCGAAACGTGATCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAATATTTGATTGCGGATTCGCTAGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                                            <  1:576209/139‑1 (MQ=14)
                                                                      aCCTTAATCGAAGCTTCTCATCCTTCCCCGAATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                                                   <  2:901839/139‑1 (MQ=14)
                                                                                                   gCATGGGCAGTATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACTGTCTGCTCCCTTTGGCCGCTCGGGAAc                       <  1:695888/138‑1 (MQ=18)
                                                                                                   gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcc                      <  2:1046083/139‑1 (MQ=21)
                                                                                                   gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                       >  2:695888/1‑138 (MQ=18)
                                                                                                       gggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                             >  1:203075/1‑112 (MQ=17)
                                                                                                       gggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                             <  2:203075/112‑1 (MQ=17)
                                                                                                        ggCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          <  2:100285/66‑1 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          >  1:100285/1‑66 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          >  2:176679/1‑66 (MQ=255)
                                                                                                        ggCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                          <  1:176679/66‑1 (MQ=255)
                                                                                                         gCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgtt                                                                                               >  2:1021993/1‑58 (MQ=255)
                                                                                                         gCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                               <  1:1021993/60‑1 (MQ=255)
                                                                                                                       tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                                                                                                       >  1:866391/1‑38 (MQ=255)
                                                                                                                       tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                                                                                                       <  2:866391/38‑1 (MQ=255)
                                                                                                                       tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAt  >  2:979601/1‑139 (MQ=18)
                                                                                                                         cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg                                                         >  1:245389/1‑82 (MQ=14)
                                                                                                                         cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg                                                         <  2:245389/82‑1 (MQ=17)
                                                                                                                         cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                                                   >  2:845337/1‑88 (MQ=18)
                                                                                                                                      |                                                                                                                            
CTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAT  >  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: