Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG53.1% ‑3.0 / 57.8 49intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (12/14);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGG  >  NC_000913/1286835‑1287077
                                                                                                                                  |                                                                                                                
gCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                          >  1:482695/1‑139 (MQ=255)
    aaTATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                      <  2:482695/139‑1 (MQ=255)
          tGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                              >  2:836217/1‑125 (MQ=255)
          tGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                              <  1:836217/125‑1 (MQ=255)
                cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                       <  1:358982/126‑1 (MQ=255)
                cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                       >  2:358982/1‑126 (MQ=255)
                                        cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        <  1:760092/101‑1 (MQ=35)
                                        cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        >  2:760092/1‑101 (MQ=35)
                                        cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                  <  2:1320176/139‑1 (MQ=37)
                                                   gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                  >  1:294009/1‑128 (MQ=16)
                                                   gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                  <  2:294009/128‑1 (MQ=16)
                                                     tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCa                                                                      >  1:1234446/1‑122 (MQ=255)
                                                     tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCa                                                                      <  2:1234446/122‑1 (MQ=255)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        >  1:440544/1‑82 (MQ=37)
                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        <  2:440544/82‑1 (MQ=38)
                                                            gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                  <  2:1014636/119‑1 (MQ=255)
                                                            gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATggg                                                                  >  1:1014636/1‑119 (MQ=255)
                                                               ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCg                                           >  2:748270/1‑139 (MQ=255)
                                                                 aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                         >  1:141051/1‑139 (MQ=255)
                                                                  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        <  1:1170651/75‑1 (MQ=39)
                                                                  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                        >  2:1170651/1‑75 (MQ=38)
                                                                  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                 >  1:1216445/1‑98 (MQ=255)
                                                                  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCAt                                                                                 <  2:1216445/98‑1 (MQ=255)
                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                              >  1:1312690/1‑115 (MQ=255)
                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                              <  2:1312690/115‑1 (MQ=255)
                                                                     aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                     <  2:243722/139‑1 (MQ=255)
                                                                      ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCgttgggaa                                    <  1:748270/139‑8 (MQ=11)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                               <  1:623489/120‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                               >  2:623489/1‑120 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                            <  1:425380/139‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            <  2:1227838/139‑1 (MQ=255)
                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            >  2:61624/1‑139 (MQ=11)
                                                                                 tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGaa                                                                                          <  2:909551/74‑1 (MQ=255)
                                                                                 tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGaa                                                                                          >  1:909551/1‑74 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                       >  1:51642/1‑139 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            >  2:775593/1‑134 (MQ=255)
                                                                                   cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            <  1:775593/134‑1 (MQ=255)
                                                                                        ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCa                                                                      <  1:575554/87‑1 (MQ=255)
                                                                                        ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCa                                                                      >  2:575554/1‑87 (MQ=255)
                                                                                             ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                >  2:59222/1‑88 (MQ=255)
                                                                                             ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                <  1:59222/88‑1 (MQ=255)
                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGG  <  2:141051‑M1/139‑1 (MQ=255)
                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGG  <  1:61624‑M1/139‑1 (MQ=255)
                                                                                                                  cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCg                                                                        >  2:168467/1‑59 (MQ=14)
                                                                                                                  cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCg                                                                        <  1:168467/59‑1 (MQ=14)
                                                                                                                  cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            >  2:701849/1‑103 (MQ=14)
                                                                                                                  cAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                            <  1:701849/103‑1 (MQ=11)
                                                                                                                       ctctctCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                <  1:454608/62‑1 (MQ=11)
                                                                                                                       ctctctCGCTGGCGCTCGAGCCGAACCTTAATCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                >  2:454608/1‑62 (MQ=14)
                                                                                                                                  |                                                                                                                
GCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGG  >  NC_000913/1286835‑1287077

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: ‑