Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,6920AT21.4% 53.0 / 11.7 28intergenic (+353/‑163)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (10/12);  new base T (3/3);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTC  >  NC_000913/2653555‑2653835
                                                                                                                                         |                                                                                                                                               
gTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACAGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCgg                                                                                                                                      <  4:24916/149‑1 (MQ=255)
    taaaataaCCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGtttt                                                                                                                                  <  5:163894/149‑1 (MQ=255)
                gCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTgg                                                                                                                        <  7:123529/147‑1 (MQ=255)
                         tCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTcgcg                                                                                                             >  5:30891/1‑149 (MQ=255)
                                      aaaTATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttttt                                                                                                                                <  7:23001/117‑1 (MQ=38)
                                      aaaTATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttttt                                                                                                                                >  8:23001/1‑117 (MQ=38)
                                            aaataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttttt                                                                                                                                 >  6:5438/1‑110 (MQ=255)
                                            aaataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGttttt                                                                                                                                 <  5:5438/110‑1 (MQ=255)
                                            aaataaatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTa                                                                                          >  4:103166/1‑149 (MQ=255)
                                                 aatCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTAcc                                                                                     <  3:103166/149‑1 (MQ=255)
                                                  atCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTg                                                                                                                         <  1:168988/112‑1 (MQ=255)
                                                  atCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTg                                                                                                                         >  2:168988/1‑112 (MQ=255)
                                                         taaaaaTTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttttt                                                                                                                                <  3:66315/96‑1 (MQ=37)
                                                         tCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGtttttt                                                                                                                                >  4:66315/1‑98 (MQ=37)
                                                                               gTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCTGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGgaga                                                        >  7:72237/1‑148 (MQ=255)
                                                                                   aaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGttttt                                                                                                                                 >  7:71061/1‑71 (MQ=21)
                                                                                   aaGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCTGATTTTGGCGGttttt                                                                                                                                 <  8:71061/71‑1 (MQ=21)
                                                                                            cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcagc                                                                                                  <  5:11335/93‑1 (MQ=255)
                                                                                            cTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCcagc                                                                                                  >  6:11335/1‑93 (MQ=255)
                                                                                             ttAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGt                                          >  1:174658/1‑148 (MQ=255)
                                                                                                     gTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATAc                                                  <  7:102615/132‑1 (MQ=255)
                                                                                                     gTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATAc                                                  >  8:102615/1‑132 (MQ=255)
                                                                                                     gTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCtt                                 >  4:45716/1‑149 (MQ=255)
                                                                                                            cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTc                           <  6:30891/148‑1 (MQ=255)
                                                                                                            cACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTc                           <  8:72237/148‑1 (MQ=255)
                                                                                                                   cTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCaaaa                   <  3:45716/149‑1 (MQ=255)
                                                                                                                              aaaaaacCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGttatta        >  1:48591/1‑149 (MQ=255)
                                                                                                                                    ccGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTc  <  2:174658/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                               
GTAATAAAATAACCTGGCTTATTAGTCCGAATTCAGACAAATATAAATAAATCCTGCTCAAAATTAAAAATTCTAACCGGTAAAAGATATTACTTAAACATGTAAATTCACTTTCCTTTAAAAAACAAAAAACCGCCAAAATCAGGCGGTTTTTTGTTGCTGGTCCGGTTCGCGGCCTTTCCAGCAGGTTGTATTACCGTAGTAATGCAAGCGCGTCTCAGCGGAGACAATACTCGCCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTTATTATTGTTC  >  NC_000913/2653555‑2653835

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