Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,415,887 (G)6→7 coding (1613/1626 nt) aidB → DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,415,8811.G97.9% 167.5 / ‑2.5 48coding (1607/1626 nt)aidBDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (24/23);  total (25/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 = 9.15e-01

CATCAGCTATTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGT  >  NC_000913/4415733‑4416007
                                                                                                                                                     |                                                                                                                              
catcaGCTATTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACgg                                                                                                                              >  2:504332/1‑151 (MQ=255)
 atcaGCTATTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACggg                                                                                                                             <  2:75400/151‑1 (MQ=255)
       tATTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACggggggg                                                                                                                         >  1:49621/1‑148 (MQ=255)
        aTTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCTAGAATGATTTATTGCTGCGGGCGACggggggga                                                                                                                        <  2:293205/149‑3 (MQ=255)
          tCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAgtgt                                                                                                                     >  2:401677/1‑151 (MQ=255)
           ccTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAgtgtg                                                                                                                    >  2:309836/1‑151 (MQ=255)
                      gCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCg                                                                                                          <  1:63296/150‑1 (MQ=255)
                         gTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGtata                                                                                                      <  2:308160/151‑1 (MQ=255)
                           gcgcAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACg                                                                                                    <  1:232633/151‑1 (MQ=255)
                            cgcAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGa                                                                                                   <  2:220261/151‑1 (MQ=255)
                              cAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACt                                                                                                 >  2:274227/1‑151 (MQ=255)
                               aaaTGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACt                                                                                                 <  1:107018/150‑1 (MQ=255)
                                 ctGTTGAAATATGCTTCTCCGCCAATGTCGCAGGCGTGGTGTCAGGTGATGTTAGATACGAGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTg                                                                                                >  2:244870/2‑149 (MQ=255)
                                   gTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGc                                                                                            <  1:312717/151‑1 (MQ=255)
                                   gTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCg                                                                                           >  1:226882/1‑151 (MQ=255)
                                     tGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCg                                                                                           <  1:135949/150‑1 (MQ=255)
                                       aaaTATGCTTCTGCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACg                                                                                        >  2:536673/1‑151 (MQ=255)
                                              cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGttt                                                                                 >  2:5240/1‑151 (MQ=255)
                                                tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGttt                                                                                 <  2:24879/149‑1 (MQ=255)
                                                tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCg                                                                               >  2:95469/1‑151 (MQ=255)
                                                  tcCGCCAATGCAGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGACGGCGGACGGTTGTCAGAGCAGATCAAGAATGATTTATTGCTGAGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTAAACGACTGATGCGTCGCTGGTTTCGa                                                                              >  2:447612/1‑150 (MQ=255)
                                                     gCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTa                                                                           >  1:431740/1‑150 (MQ=255)
                                                       cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACt                                                                        >  1:275616/1‑151 (MQ=255)
                                                         aTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTaa                                                                      <  2:586020/151‑1 (MQ=255)
                                                           ggCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTaaa                                                                     <  2:279491/150‑1 (MQ=255)
                                                           ggCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGATGGTTTCGATTAACTAAAt                                                                    <  1:17863/151‑1 (MQ=255)
                                                               cAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGaa                                                                 >  2:147745/1‑150 (MQ=255)
                                                                           cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                     <  1:583063/150‑1 (MQ=255)
                                                                              gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                  >  1:223928/1‑150 (MQ=255)
                                                                              gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTa                                                 <  2:355832/151‑1 (MQ=255)
                                                                                gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAgg                                               <  1:100485/151‑1 (MQ=255)
                                                                                gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAg                                                <  2:642217/150‑1 (MQ=255)
                                                                                gATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAg                                                <  1:164900/150‑1 (MQ=255)
                                                                                        aTACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACa                                        >  2:465976/1‑150 (MQ=255)
                                                                                         tACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACa                                        >  2:125597/1‑149 (MQ=255)
                                                                                          aCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                      <  2:309785/150‑1 (MQ=255)
                                                                                                   ggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcgcg                             <  1:546310/150‑1 (MQ=255)
                                                                                                    gcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgcc                           >  1:502592/1‑151 (MQ=255)
                                                                                                        aCGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATGGAAGCCCGGACAAGGGGCCCGCGCGGCa                        >  2:309843/1‑150 (MQ=255)
                                                                                                               tCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAt                >  2:4331/1‑151 (MQ=255)
                                                                                                                 agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg              >  2:263795/1‑151 (MQ=255)
                                                                                                                 agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg              <  2:269014/151‑1 (MQ=255)
                                                                                                                 agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg              >  2:513216/1‑151 (MQ=255)
                                                                                                                 agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg              <  2:354559/151‑1 (MQ=255)
                                                                                                                  gagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGt             >  2:596417/1‑151 (MQ=255)
                                                                                                                    gCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTATGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGtt            <  1:144849/150‑1 (MQ=255)
                                                                                                                      aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt         >  1:395918/1‑151 (MQ=255)
                                                                                                                      aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAt         <  1:48945/151‑1 (MQ=255)
                                                                                                                        aTCCAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAttt       <  1:85530/151‑1 (MQ=255)
                                                                                                                           cAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCt    >  1:324379/1‑151 (MQ=255)
                                                                                                                           cAGAATGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCt    >  2:645292/1‑151 (MQ=255)
                                                                                                                                tGATTTATTGCTGCGGGCGACGGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGt  >  2:666356/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                              
CATCAGCTATTCCTGCTGGGCTGCGGTGCGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGT  >  NC_000913/4415733‑4416007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: