Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0440GA41.4% 23.4 / 27.0 29intergenic (+518/+22)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base A (5/7);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAA  >  NC_000913/1286913‑1287156
                                                                                                                                   |                                                                                                                
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                                                           <  1:425380/139‑1 (MQ=255)
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           >  2:61624/1‑139 (MQ=11)
cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           <  2:1227838/139‑1 (MQ=255)
     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                                                                      >  1:51642/1‑139 (MQ=255)
     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           >  2:775593/1‑134 (MQ=255)
     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           <  1:775593/134‑1 (MQ=255)
                          ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGG                                                                                 <  2:141051‑M1/139‑1 (MQ=255)
                          ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGG                                                                                 <  1:61624‑M1/139‑1 (MQ=255)
                                    cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           >  2:701849/1‑103 (MQ=14)
                                    cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                           <  1:701849/103‑1 (MQ=11)
                                                 ctggcgctcGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                                                            >  1:1041976‑M1/5‑121 (MQ=255)
                                                 ctggcgctcGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                                                            <  2:1041976‑M1/117‑1 (MQ=255)
                                                          gAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTg                                                                                                    <  1:125533/88‑1 (MQ=14)
                                                          gAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTg                                                                                                    >  2:125533/1‑88 (MQ=11)
                                                                  ccTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttaca                                         <  2:53534‑M1/139‑34 (MQ=255)
                                                                          gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctc                                 >  1:1135185‑M1/1‑98 (MQ=255)
                                                                            aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccc                               >  2:521319‑M1/1‑96 (MQ=255)
                                                                                          ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                                                                                <  2:470848‑M1/76‑1 (MQ=255)
                                                                                          ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                                                                                >  1:470848‑M1/1‑76 (MQ=255)
                                                                                              gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                                                  >  1:796207/1‑54 (MQ=255)
                                                                                              gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                                                  <  2:796207/54‑1 (MQ=255)
                                                                                                   ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccacc        <  1:521319‑M1/139‑67 (MQ=255)
                                                                                                    gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcg                                                                >  2:394997‑M1/1‑72 (MQ=255)
                                                                                                    gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcg                                                                <  1:394997‑M1/82‑11 (MQ=255)
                                                                                                         atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccgctcgggaaccccaccaggggt  <  2:1135185‑M1/139‑73 (MQ=255)
                                                                                                                    cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccg                       <  2:16114‑M1/107‑52 (MQ=255)
                                                                                                                    cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctccctttggccg                       >  1:16114‑M1/1‑56 (MQ=255)
                                                                                                                    cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGG                                                                                          >  2:475450‑M1/1‑40 (MQ=255)
                                                                                                                    cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGG                                                                                          <  1:475450‑M1/40‑1 (MQ=255)
                                                                                                                                   |                                                                                                                
CGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAA  >  NC_000913/1286913‑1287156

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