Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,530 T→C 8.1% intergenic (‑166/‑38) moeA ← / → iaaA molybdopterin molybdotransferase/isoaspartyl dipeptidase proenzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913866,5300TC 8.1% 158.4 / 3.2 60intergenic (‑166/‑38)moeA/iaaAmolybdopterin molybdotransferase/isoaspartyl dipeptidase proenzyme
Reads supporting (aligned to +/- strand):  ref base T (37/17);  new base C (3/2);  total (40/20)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGAGGCGTT  >  NC_000913/866399‑866662
                                                                                                                                   |                                                                                                                                    
cTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAACATTCTGATATTTATAGCAAAAGTgg                                                                                                                           >  2:442564/1‑143 (MQ=255)
 tattatGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAACGTGGc                                                                                                                          >  1:809964/1‑143 (MQ=255)
 tattatGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGc                                                                                                                          >  2:795206/1‑143 (MQ=255)
  attatGTGAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCg                                                                                                                         >  2:381909/1‑143 (MQ=255)
    tatGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTTCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGaa                                                                                                                       <  2:910819/143‑1 (MQ=255)
     atGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGaaa                                                                                                                      >  2:570140/1‑142 (MQ=255)
        tCAGAAAAACGCCACAGACTGTATGCCACCTCGCGCGTAGCGCTGGATCCTGCCTTTACATGCCACACCCATCTCTCTATCTTCAACAATTGAATGAGGACTTCATAAAAATTCTGACATCTACAGCAAAAGTAGCGCaccca                                                                                                                   >  2:122359/1‑141 (MQ=255)
          aGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAAATCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAacaacca                                                                                                                 >  2:701785/1‑142 (MQ=255)
             aaaaCGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTa                                                                                                              <  1:980424/143‑1 (MQ=255)
                 cGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAACATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATgg                                                                                                          >  1:817641/1‑143 (MQ=255)
                   ccACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAACACTTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCCAACCACCCTTAATGGAc                                                                                                        >  2:416108/1‑143 (MQ=255)
                   ccACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTACAGCAAAAGTGGCGAACCACCCTTAATGGAc                                                                                                        >  1:998751/1‑143 (MQ=255)
                   ccACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTACTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTACAGCAAAAGTGGCGAACCACCCTTAATGGAc                                                                                                        >  1:201299/1‑143 (MQ=255)
                           tGTATGCCACCTCGGGCGTAGCGCTGGGTACTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAGTTGAATGAGTAATTCATAAAAATTCTGATATTTATCGCAAGAGTGGGGAACCAGCCTTCATGGACGAATAcga                                                                                                >  1:738973/1‑141 (MQ=255)
                           tGTATGCCACCTCGGGCGTAGCGCTGGATCCTCCCTTTACGTGTCATATGCATCTTCCTATATTCAAAAATTGAAGAATTAGTTGATAATAATTCTGATATTGATAGCAAGAGTGGTGAAAGACCCTTAAAGGACGAATCCTa                                                                                                >  2:179095/1‑143 (MQ=255)
                             tATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATCTCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAGTAATGCTGATGTTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATg                                                                                              >  1:340157/1‑143 (MQ=255)
                             tATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATg                                                                                              <  2:652636/143‑1 (MQ=255)
                                 ccACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTACTGGACGAGTACTATGGGCa                                                                                          >  1:70470/1‑143 (MQ=255)
                                            gTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCAtt                                                                               <  1:578186/143‑1 (MQ=255)
                                            gTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCAtt                                                                               <  1:802688/143‑1 (MQ=255)
                                             tAGCGCTGGGTCCTGCATTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTg                                                                              <  2:951603/143‑1 (MQ=255)
                                                     ggTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAAGTCATAAAAATTCTTATCTTTATAGCAAAAGGGGCGAACCACCCATAATGGACGAATACTAAGAGCAAAGCAGTCATTGCAATTCaa                                                                      >  2:918639/1‑142 (MQ=255)
                                                      gTCCTGTCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCAtg                                                                     <  2:200461/143‑1 (MQ=255)
                                                       tCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCAtgg                                                                    >  2:61696/1‑143 (MQ=255)
                                                         cTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGACCCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCAtggtg                                                                  >  1:558756/1‑143 (MQ=255)
                                                                   aTGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGc                                                        >  2:844083/1‑143 (MQ=255)
                                                                    tGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCa                                                       >  2:728353/1‑143 (MQ=255)
                                                                    tGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCa                                                       >  1:940212/1‑143 (MQ=255)
                                                                     ggccaatccatccGTCTATATTCAAAAATTGAATGAGTAATTCACAAAAATGCTGATATCTAAAGCAAATTTGGCGAACCACCCTTAATGGACGAAGACTAGGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCaa                                                      <  2:790352/138‑1 (MQ=255)
                                                                         tatCCATCTTTGTATATGCTAAAATTGAATGAGTAATTCATAAAAATTCTGAGATTTATAGCAAAAGTGGCGAACCACCCTTAATTGACGAAGACGATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAg                                                  <  2:434701/143‑1 (MQ=255)
                                                                          atCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAAGTAGc                                                 >  2:560344/1‑143 (MQ=255)
                                                                           tCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCGGGTGCAATTAGcc                                                >  1:912209/1‑143 (MQ=255)
                                                                              atgtgtCTATCTGCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCcgcg                                             <  2:284997/138‑1 (MQ=255)
                                                                              aTCTTTCTATATTCAAAAATTTAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGGGGCGAACCAACATTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCcgcg                                             <  2:820220/143‑1 (MQ=255)
                                                                                    cTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTCATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCAGAGCGCata                                       >  1:405431/1‑141 (MQ=255)
                                                                                    cTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGa                                       >  1:405468/1‑143 (MQ=255)
                                                                                     tataTTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGAt                                      <  1:103067/143‑1 (MQ=255)
                                                                                     tataTTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGAt                                      >  2:70962/1‑143 (MQ=255)
                                                                                     tataTTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGAt                                      >  2:189051/1‑143 (MQ=255)
                                                                                       tatGCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGa                                    <  1:249216/143‑1 (MQ=255)
                                                                                          tCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTc                                 >  1:771010/1‑143 (MQ=255)
                                                                                          tCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAAATCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTc                                 >  2:312753/1‑143 (MQ=255)
                                                                                            aaaaaTTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTg                               >  1:997678/1‑143 (MQ=255)
                                                                                            aaaaaTTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTg                               >  1:661449/1‑143 (MQ=255)
                                                                                             aaaaTTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGc                              >  2:605985/1‑143 (MQ=255)
                                                                                                 ttGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACa                          >  1:608241/1‑143 (MQ=255)
                                                                                                    aaTGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGTAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGa                       <  1:381909/143‑1 (MQ=255)
                                                                                                      tGAGTAATTCATAAAAATTCCGATATTTATAGCCAAAGTGGCGAACCACCCTTAATGGACGAATCCTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAAt                     <  1:563866/143‑1 (MQ=255)
                                                                                                         gTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCATATGAGTCTGAAACAGGAAtgaa                  >  1:399118/1‑140 (MQ=255)
                                                                                                           aaTTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAAGCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAAGGAGCCGCGGGCAGATGAGTCTGCAACAGGAAGTACGc                >  2:1030572/1‑143 (MQ=255)
                                                                                                           aaTTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATCGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGCTGCAATTAGCCGCTCGCAGATGAGGCTGGAACAGGACTTACGc                >  2:1030577/1‑143 (MQ=255)
                                                                                                             ttCATAAAAATTCTGATATTTACAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTa              <  2:201299/143‑1 (MQ=255)
                                                                                                             ttCATAAAAATTCTGATATTTACAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTa              <  2:998751/143‑1 (MQ=255)
                                                                                                                aTAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAAGTAGCCGCGCGGAGATGAGTCTGCAACAGGAGTTACGCTAGAt           >  1:464978/1‑143 (MQ=255)
                                                                                                                 tAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACAta          >  1:272190/1‑142 (MQ=255)
                                                                                                                 tAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATc          <  1:578792/143‑1 (MQ=255)
                                                                                                                   aaaaTTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTAGGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGa        >  2:1003889/1‑143 (MQ=255)
                                                                                                                    aaaTTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATAGTGGCGCAGGTGCAGTTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGAg       >  2:968853/1‑143 (MQ=255)
                                                                                                                       ttCTGATGTTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGAGGCg    <  1:616223/143‑1 (MQ=255)
                                                                                                                         cTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGAGGCGtt  <  2:809964/143‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAAAGTGGCGAACCACCCTTAATGGACGAATACTATGGGCAAAGCAGTCATTGCAATTCATGGTGGCGCAGGTGCAATTAGCCGCGCGCAGATGAGTCTGCAACAGGAATTACGCTACATCGAGGCGTT  >  NC_000913/866399‑866662

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