Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,885 T→C 11.6% intergenic (+359/+181) 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,8850TC11.6% 116.6 / 5.0 53intergenic (+359/+181)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/4);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGG  >  NC_000913/1286736‑1287031
                                                                                                                                                     |                                                                                                                                                  
gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGt                                                                                                                                                    <  1:1333923/150‑1 (MQ=255)
 cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAGCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGATTTCCGGGGCTTTTGAAAGTg                                                                                                                                                   <  1:1343367/150‑1 (MQ=255)
 cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTg                                                                                                                                                   <  2:1338493/150‑1 (MQ=255)
  ggcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGa                                                                                                                                                  <  2:1616503/149‑1 (MQ=255)
       gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtg                                                                                                                                             <  2:721988/150‑1 (MQ=255)
           cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtgg                                                                                                                                         <  2:862002/150‑1 (MQ=255)
           cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCAACTCGACCCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCCGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtgg                                                                                                                                         <  2:1213166/150‑1 (MQ=17)
           cGCCCTTCGGGTCGGTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTGCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGTTTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGAtggtggtgg                                                                                                                                         <  1:1915721/150‑1 (MQ=14)
            gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                                                                                                                        <  2:1062811/150‑1 (MQ=255)
             cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                                                                                       <  2:1218039/150‑1 (MQ=255)
             cccTCCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGCCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAATTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                                                                                       <  2:64798/150‑1 (MQ=255)
             cccATCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCATGGGCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTCATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                                                                                                                                       <  2:415960/150‑1 (MQ=14)
               cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGATTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGa                                                                                                                                     >  2:1081062/1‑150 (MQ=255)
                    ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                                                                <  1:965567/150‑1 (MQ=255)
                              tGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgcgc                                                                                                                      >  1:591668/1‑150 (MQ=255)
                                       gCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGGCTAGTGAAAGGGAAGGTCGTGGGGGAAGGCTAAATCCGCGCTgcgcggtg                                                                                                             >  1:604999/1‑147 (MQ=14)
                                                cTGGCGCTCAACCCGAACCTTGTGCGAATCTTCTCATCCTTCCCCGCTGTGGCAGATTCTTTGATTGCGGCTTCGTTTGAGAATTCCTGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg                                                                                                    <  2:1862315/150‑1 (MQ=14)
                                                                   ttGGGCTAAGATTCTCATCCTACCCCGCTTGGCCAGAAGCATAGATTGCCGAGTCGTATCAGAATTCCGTGTCTTATGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAAc                                                                                 <  1:858195/150‑1 (MQ=14)
                                                                     gcgcctaGCTTCTCCTCGTTACCCGCGGGGGCAGATGTTATCATTGGGTCTTCGTTTGAGTATTCCGGTTTTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                                               <  2:536382/144‑1 (MQ=11)
                                                                      gTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct                                                                              >  1:2005032/1‑150 (MQ=255)
                                                                        cGAACCTTCTCCTCCATCCCCGCGTGGGCAGACGCTTTCATTCGCGCCTCGTTTCAGTATTCCTGGGCTTTAGAAAGAGATGGTTGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCGCAGGCAACgaccc                                                                            <  1:1092022/150‑5 (MQ=9)
                                                                             cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGATTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc                                                                                                         <  2:1783192/116‑1 (MQ=255)
                                                                             cTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc                                                                                                         >  1:1783193/1‑116 (MQ=255)
                                                                              ttCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                      >  2:1172594/1‑150 (MQ=255)
                                                                                ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                                                                                                   >  2:72762/1‑103 (MQ=255)
                                                                                   aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                 >  1:1781025/1‑150 (MQ=255)
                                                                                       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGcac                                                                             >  2:823054/1‑132 (MQ=255)
                                                                                       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctc                                                                             <  1:823055/134‑1 (MQ=255)
                                                                                                tGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                        >  2:1188520/1‑130 (MQ=40)
                                                                                                      gAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTTACGCTCGAAGCGAACCTTAGTCg                                              >  2:654847/1‑150 (MQ=14)
                                                                                                          atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTggc                                                                   >  2:241245/1‑123 (MQ=39)
                                                                                                          atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTggc                                                                   <  1:241245/125‑3 (MQ=39)
                                                                                                          atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTATGTCGAAGc                                          >  2:233569/1‑150 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                                                                                                        <  2:343758/81‑1 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc                                                                                                        >  1:343758/1‑81 (MQ=255)
                                                                                                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAt                                   <  1:1172595/150‑1 (MQ=12)
                                                                                                                    ggctgcctttGTGTATTCCGGTGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                <  1:1591853/143‑1 (MQ=14)
                                                                                                                    ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCt                                >  1:850298/1‑150 (MQ=255)
                                                                                                                         cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                       <  2:1203147/106‑1 (MQ=255)
                                                                                                                         cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                       >  1:1203148/1‑106 (MQ=255)
                                                                                                                         cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtct                                      >  1:301278/1‑139 (MQ=12)
                                                                                                                          gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCATGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTccccc                          >  2:86481/1‑149 (MQ=255)
                                                                                                                            ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                        >  2:1057075/1‑150 (MQ=255)
                                                                                                                                gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg                                                                                            <  2:41038/78‑1 (MQ=255)
                                                                                                                                gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg                                                                                            >  1:41038/1‑78 (MQ=255)
                                                                                                                                       ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTggc                                                                   >  1:382648/1‑94 (MQ=14)
                                                                                                                                       ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTggc                                                                   <  2:382648/96‑3 (MQ=14)
                                                                                                                                           tttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg         <  2:723565/143‑1 (MQ=255)
                                                                                                                                             tttGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCGCGCATGGGCAGAATATTTggc       <  1:1550050/150‑3 (MQ=255)
                                                                                                                                             taaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                   >  1:383193/6‑122 (MQ=16)
                                                                                                                                             taaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                   <  2:383193/117‑1 (MQ=16)
                                                                                                                                                ggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                    <  1:1229126/116‑1 (MQ=255)
                                                                                                                                                ggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                    >  2:1229125/3‑118 (MQ=255)
                                                                                                                                                  aaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCgg  <  2:819136/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGG  >  NC_000913/1286736‑1287031

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