Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,2190GA46.2% 3.8 / 27.2 28noncoding (18/171 nt)rttRrtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base G (6/8);  new base A (8/4);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT  >  NC_000913/1287101‑1287345
                                                                                                                       |                                                                                                                              
ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGaa                                                                                                             >  2:573217/1‑138 (MQ=11)
ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGaa                                                                                                             <  1:573217/138‑1 (MQ=11)
                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                            <  1:117264/139‑1 (MQ=255)
                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                            <  2:903/139‑1 (MQ=255)
                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                            <  2:153044/139‑1 (MQ=255)
                               ctctctCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                        >  2:364414/1‑112 (MQ=11)
                               ctctctCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                                        <  1:364414/112‑1 (MQ=11)
                                              tCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAggg                                                                                                                      >  2:786397/1‑83 (MQ=14)
                                              tCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAggg                                                                                                                      <  1:786397/83‑1 (MQ=14)
                                               cGAGACGAACCGTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATATGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGc                                                             >  1:466155/1‑139 (MQ=11)
                                                       aCCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc                                                     <  1:734265/139‑1 (MQ=14)
                                                            aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                                                                              <  2:459382/91‑1 (MQ=255)
                                                            aatCGAAGCTTCTCATCATTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACGTTTGAAAGTGATGGTGGTggg                                                                                              >  1:459382/3‑93 (MQ=255)
                                                             atCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                               <  2:772571/138‑1 (MQ=14)
                                                                  aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTccc                                          >  1:42973/1‑139 (MQ=17)
                                                                       ctcATCCTTCCCCGCATTGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGATCCCTTTgg                                     >  1:399581/1‑139 (MQ=17)
                                                                               tCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTG‑‑AG‑‑‑A‑GGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGc                                                             >  2:10182/1‑101 (MQ=14)
                                                                               tCCCCGCATGGGCAGAATATTTGATTGCGGATTCG‑CTTG‑‑AG‑‑‑A‑GGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGc                                                             <  1:10182/101‑1 (MQ=14)
                                                                                      aTGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGc                                    <  1:484993/125‑1 (MQ=18)
                                                                                      aTGGGCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATGCGAACCTTCGAAGTCGATGACGGCAGACTTACAGTCTGCTCCCTTTGGc                                    >  2:484993/1‑125 (MQ=17)
                                                                                         ggCAGAATATTTGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCccacc                   >  2:264288/1‑139 (MQ=18)
                                                                                                  tttgaTTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                         <  1:231752/120‑1 (MQ=17)
                                                                                                  tttGATTGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAc                         >  2:231752/1‑124 (MQ=17)
                                                                                                        tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                                                                                                      <  2:468586/41‑1 (MQ=255)
                                                                                                        tGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                                                                                                      >  1:468586/1‑41 (MQ=255)
                                                                                                        tGCGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAt    >  1:427418/1‑139 (MQ=255)
                                                                                                          cGGATTCGCTTGAG‑AGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAttt  >  2:455358/1‑139 (MQ=18)
                                                                                                          cGGATTCG‑CTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAAttt  >  1:737585/1‑139 (MQ=255)
                                                                                                                       |                                                                                                                              
TTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCG‑TTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTT  >  NC_000913/1287101‑1287345

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