Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,1430GA40.0% 0.1 / 21.2 25noncoding (94/171 nt)
P12S (CCA→TCA) 
rttR
tpr
rtT sRNA, processed from tyrT transcript
protamine‑like protein
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base A (3/7);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA  >  NC_000913/1287009‑1287260
                                                                                                                                      |                                                                                                                     
aTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                   <  2:1400624/139‑1 (MQ=14)
   ggCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACTGCTCCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCgc                                                                                                                >  1:942684/1‑138 (MQ=11)
                  tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCttaatcg                                                                                                 >  2:1725386/1‑134 (MQ=11)
                  tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCGAGCCGAACCttaatcg                                                                                                 >  1:895599/1‑134 (MQ=11)
                    cGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACTCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGaa                                                                                               >  2:806401/1‑139 (MQ=11)
                    cGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCt                                                                                                       <  1:1083847/131‑1 (MQ=11)
                    cGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGATCTCTCGCTGGCGCTCGAGCCGAACCt                                                                                                       >  2:1083847/1‑131 (MQ=11)
                     ggattcgttgggaagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAg                                                                                              <  1:237183/126‑1 (MQ=11)
                                       aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                            <  2:577358/139‑1 (MQ=2)
                                       aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGc                                                                            <  1:806401/139‑1 (MQ=2)
                                                      cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgcc                                                                                                             <  2:1006628/90‑3 (MQ=17)
                                                      cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgcc                                                                                                             >  1:1006628/2‑89 (MQ=17)
                                                                     aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgc                                                                                                              >  1:577045/1‑74 (MQ=17)
                                                                     aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgc                                                                                                              <  2:577045/75‑2 (MQ=21)
                                                                         aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTATTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgggt                                          <  1:1411973/139‑2 (MQ=255)
                                                                             cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                               >  2:1330820/1‑66 (MQ=35)
                                                                             cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                               <  1:1330820/66‑1 (MQ=35)
                                                                                                   ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtgg                <  1:1492201/139‑1 (MQ=255)
                                                                                                      gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg             <  2:894941/139‑1 (MQ=11)
                                                                                                                 tGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa  <  2:1193350/139‑1 (MQ=255)
                                                                                                                  gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg       <  1:396378/133‑1 (MQ=255)
                                                                                                                  gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg       >  2:396378/1‑133 (MQ=11)
                                                                                                                         cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                 <  1:179654/100‑1 (MQ=11)
                                                                                                                         cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                 >  2:179654/1‑100 (MQ=11)
                                                                                                                               ctcGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                   <  1:755147/92‑1 (MQ=11)
                                                                                                                               ctcGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                   >  2:755147/1‑92 (MQ=11)
                                                                                                                                      |                                                                                                                     
ATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA  >  NC_000913/1287009‑1287260

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