Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→G 14.4% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG14.4% 92.4 / 4.8 50intergenic (+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/26);  new base G (5/3);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCC  >  NC_000913/1286819‑1287106
                                                                                                                                                  |                                                                                                                                             
aTCCTTCCCCGCTTGGGCAGAATATGTGATTGCGGATTCGTGTGATAATTCCGGGGCTTTTGCACGTGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGAcc                                                                                                                                             >  1:767705/1‑148 (MQ=14)
  ccTTCCACGTGTGGGCATCATGTTTGATTGGGGCGGCGTTTGAGAATTTCGTGGTTTTTGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                           <  2:1339650/149‑1 (MQ=14)
    ttaccagCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGGCACGCGCTCTGGCTGGCGCTCg                                                                                                                                         >  2:1260363/7‑149 (MQ=14)
     tCCCCGCTTGGGCAGAATACTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                        <  1:371864/149‑1 (MQ=255)
         cGCTTGGGCAGAATATTTGATTGCGGATCCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCCCGCCGACGCTCGAGCCg                                                                                                                                    >  1:803758/1‑149 (MQ=255)
                gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAACTTa                                                                                                                             >  1:367445/1‑149 (MQ=255)
                 cAGAATATCGGATTGCGGATTCGTTTGAGATGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                                                                                            <  1:1272812/149‑1 (MQ=255)
                        tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGc                                                                                                                     >  1:1022952/1‑149 (MQ=255)
                           gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGTCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                  <  1:895860/149‑1 (MQ=255)
                            aTTGCGGCTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtct                                                                                                                 <  1:426753/149‑1 (MQ=255)
                                 ggATGCGTTTGAGAATTCCGGGACTTATGAAAGCGATGGTGGTGGGGGAAGGCCTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                            <  1:1425686/149‑1 (MQ=14)
                                  gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCACGCGGCAACGCTCTCTCGCTGACGCCCGAGCCCAACCCTCGTCGAAGCTTCTCATccc                                                                                                           >  2:1319259/1‑148 (MQ=14)
                                        tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGGTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGGGGGAACGCGGTGTCGCTGGCGCCCGAGCCGAACCTAGGTCGAAGCTTCTCATCCTTCCCCg                                                                                                     >  2:427382/1‑149 (MQ=11)
                                             ggcttCCGGGGGTTTTGCAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTCCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                <  2:106340/146‑1 (MQ=14)
                                             gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                >  1:635118/1‑149 (MQ=255)
                                                ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGcc                                                                                             >  1:1047485/1‑148 (MQ=255)
                                                   cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGctcca                                                                                          >  2:1100873/1‑145 (MQ=255)
                                                     ggggTTGTGGAAGGGATGGTGGTGGGGGAAGGATTACTCTGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                        <  2:749716/149‑1 (MQ=11)
                                                        cTTTCGAAACCGGTGGTGGCGGGGGAACCCTTCCCCACCGCCGCCCGCTCCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                     <  2:1184482/149‑1 (MQ=11)
                                                        cTTATGGAAGCGGTGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                     <  2:238001/149‑1 (MQ=11)
                                                               aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                              <  1:807400/149‑1 (MQ=255)
                                                                 gTGATGGTGGTGGGGCAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGACGACGCTCGCGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                            <  1:707096/149‑1 (MQ=255)
                                                                        tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                     <  1:842469/149‑1 (MQ=255)
                                                                          gtggGGGAAGGATTACTCAGCGCTGCGCGATTCGTGATATGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTCCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtggg                                                                   >  1:1440198/1‑146 (MQ=14)
                                                                             ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGggag                                                                >  1:181548/1‑149 (MQ=255)
                                                                              gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCCCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgcgact                                                               >  2:687255/1‑147 (MQ=255)
                                                                              gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTGGCCTGCGGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGt                                                               <  1:582396/149‑1 (MQ=255)
                                                                               ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTTTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                              >  2:670851/1‑149 (MQ=255)
                                                                                  aGGATTCGTGAGCGCTGCGCGTTTCGCCGTTCGGGTCGTTGCCGGGGGCTACGCTCTCTCGTTCACGCTTGAGCCGAACCTCAGTCGTAGGTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                           <  2:767705/149‑1 (MQ=17)
                                                                                  aGGATTACTCAGCGCGGCGCGCTCCGCCCCTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                           <  2:638758/149‑1 (MQ=255)
                                                                                        cccgaccgccgcGCGCTCCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                     <  2:1418401/143‑1 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACtttt                                                   >  1:1182782/1‑149 (MQ=255)
                                                                                               gcTGCGCGCTTCCCCCTCCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                              <  1:1286086/149‑1 (MQ=255)
                                                                                                      gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCGGGGACTTTTGAAAGTGAtggt                                       >  2:784142/1‑149 (MQ=255)
                                                                                                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                   <  2:1047485/149‑1 (MQ=255)
                                                                                                            cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                     >  1:1349238/1‑97 (MQ=255)
                                                                                                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                             <  1:769352/149‑1 (MQ=255)
                                                                                                                  gggTCGTTGCCTGCGGCGACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                           <  1:1394932/149‑1 (MQ=255)
                                                                                                                  gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                           <  1:178087/149‑1 (MQ=255)
                                                                                                                   ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                          <  1:1285618/149‑1 (MQ=255)
                                                                                                                   cgccgctGCCTGCGGCTACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGa                          <  2:1010934/143‑1 (MQ=14)
                                                                                                                      cggtGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                       <  1:1127100/146‑1 (MQ=255)
                                                                                                                      cGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACTGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                       <  1:1127222/149‑1 (MQ=255)
                                                                                                                      cGTTACCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                       <  1:670851/149‑1 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                               >  1:978796/1‑123 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                               <  2:978796/123‑1 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGTAGGAtggct                     >  1:908606/1‑145 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTAca                     >  2:391537/1‑148 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCCAACCTTAGTCGAAGCTTCTCATCCTGCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                     >  1:945304/1‑149 (MQ=255)
                                                                                                                                           cgcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGTGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGTAAGGATTACTCAGCGCTGCGCGCTTCGcc  <  1:947519/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCC  >  NC_000913/1286819‑1287106

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