Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG23.1% 56.9 / 12.4 45intergenic (+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 (16/14);  new base G (5/4);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGC  >  NC_000913/1286819‑1287093
                                                                                                                                                  |                                                                                                                                
aTCCTTCACCGTTGGGGCAGAATTTTTAATTTCGGAGTAGTTAGAGTTTTGGTGTTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGc                                                                                                                               <  2:710667/150‑1 (MQ=11)
  ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                             <  1:660916/150‑1 (MQ=255)
       cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                        <  1:1512605/150‑1 (MQ=255)
            ttGGGCAGAATATTTGATTGCGGCGTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                   <  1:1816743/150‑1 (MQ=255)
             tGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGACGCTAGAGACGAAccg                                                                                                                  >  1:603651/1‑149 (MQ=255)
                           gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGACGAACCTTGGTCGAAGCTtca                                                                                                    >  1:49693/1‑149 (MQ=255)
                                       gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgtc                                                                                                                        >  1:1183888/1‑116 (MQ=25)
                                       gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgtc                                                                                                                        <  2:1183888/118‑3 (MQ=25)
                                        tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGACGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGc                                                                                       >  1:92935/1‑150 (MQ=12)
                                           gagaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                    >  2:495864/1‑150 (MQ=255)
                                             gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                              >  2:1613061/1‑106 (MQ=25)
                                             gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                              <  1:1613061/106‑1 (MQ=25)
                                             gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                  >  1:1866136/1‑150 (MQ=255)
                                                  ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCg                                                                                                                                <  2:1565700/99‑1 (MQ=17)
                                                  ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCg                                                                                                                                >  1:1565700/1‑99 (MQ=17)
                                                        tttttGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                       <  1:1510018/149‑1 (MQ=38)
                                                            aggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTg                                                                   <  2:1305001/147‑1 (MQ=255)
                                                              ttaGTGATGGTGGTGGGGTAATGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCT‑‑TCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                               <  2:1215952/148‑1 (MQ=14)
                                                                    aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                                                                                  <  2:1057888/95‑1 (MQ=37)
                                                                    aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                                                                                  >  1:1057888/1‑95 (MQ=37)
                                                                    aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                           <  1:495864/150‑1 (MQ=39)
                                                                           tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCCAGACGAACCTTAGTCGAATCTTCTCATCCTTCACCGCATGGGCAAAATATATGATTGCGGATTAGCTTgaga                                                    >  2:403859/1‑150 (MQ=11)
                                                                              gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                 >  2:1903188/1‑150 (MQ=255)
                                                                                         cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCACACTCGCTGACGCGCGAGACGAACCTTAGTCGAAGCTTCTAATACATCCCCCCATTGGCAGAATATTTTATTGAGGATTAGGTTGTGAGTTAAGGGAAtttt                                      >  2:1962817/1‑150 (MQ=14)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                     <  1:403859/150‑1 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                     <  2:1542602/150‑1 (MQ=255)
                                                                                           cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAAAGTTCAGGGACTTTTGa                                    >  2:1038830/1‑150 (MQ=255)
                                                                                           cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGGTTGATAGTTCAGGGACTTTAGa                                    >  2:249923/1‑150 (MQ=14)
                                                                                              cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                 <  1:975684/150‑1 (MQ=255)
                                                                                                       cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGGTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtgg                        >  2:753307/1‑150 (MQ=255)
                                                                                                       cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtgg                        <  2:1866136/150‑1 (MQ=255)
                                                                                                           gCCCTTCGTGTCGTTTTCTTCGGCAACGCTCTCTCGCTTACGCTCGAGCCGTACCTTAGTCGAAGCTTCACATCCTTCCCCGCATGTGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  2:1559100/150‑1 (MQ=14)
                                                                                                           gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  2:1898945/150‑1 (MQ=255)
                                                                                                                 ttgcgatttGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg              <  1:873429/143‑1 (MQ=11)
                                                                                                                 cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCAAAATACCCAGCAAGGGCAGTATATTATATTGCGGATTAGCTTGAGATTTCAGGGGCGTTGGAAAGTTATATTATTGGGGGAAgg              >  2:1739437/1‑150 (MQ=9)
                                                                                                                    tgtgTTGCCGGCGGCAACGCTCTCTCGCTGTCGCTCTAGGCGAACCTTTTTCGAAGATTCTCATGCTTAGCCGCTTTGGCAGAATATTTGATTGCGGATGCGCTTGAGAGTTCAGGGAATTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           <  2:848803/147‑1 (MQ=11)
                                                                                                                      cGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCCAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc         <  1:1049775/150‑1 (MQ=255)
                                                                                                                       gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGATGGATTACt        >  2:1558020/1‑150 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGCCGCTCTAGCCGAACCTTATGCGAAGCTTCTCATCCTTCCCCACATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc       <  1:493248/150‑1 (MQ=14)
                                                                                                                         tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGGGTGGAGGGAGGTGAGATCTGGAGGGGGGTGGGGGAAGGATTTCTCa      >  1:1889068/1‑150 (MQ=14)
                                                                                                                           ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGATTGATAGTTCAGAGACGATTCAAAGTGATGTTGGTCAGGGAATAATTACTCata    >  2:622846/1‑148 (MQ=14)
                                                                                                                             tGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAgagc  >  1:658912/1‑150 (MQ=255)
                                                                                                                               cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                              <  1:1084083/88‑1 (MQ=255)
                                                                                                                               cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                              >  2:1084083/1‑88 (MQ=255)
                                                                                                                                ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   >  2:577130/1‑82 (MQ=255)
                                                                                                                                                  |                                                                                                                                
ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGC  >  NC_000913/1286819‑1287093

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