Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,615,859 T→G T403P (ACC→CCC)  yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,615,8590TG96.0% 266.3 / 2.1 75T403P (ACC→CCC) yjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (38/34);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

GTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTT  >  NC_000913/4615770‑4615955
                                                                                         |                                                                                                
gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                        >  2:400650/1‑100 (MQ=255)
gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                        <  1:144329/100‑1 (MQ=255)
gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                        <  2:326606/100‑1 (MQ=255)
 ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                       >  1:658876/1‑100 (MQ=255)
   aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                     <  1:1745931/100‑1 (MQ=255)
   aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGAATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTgcc                                                                                     >  1:1246462/1‑98 (MQ=255)
    cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                    >  2:1242207/1‑100 (MQ=255)
     gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                   >  2:415244/1‑100 (MQ=255)
     gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                   >  2:1453866/1‑100 (MQ=255)
       ttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCAt                                                                                 >  1:1184291/1‑100 (MQ=255)
         gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                               >  2:823652/1‑100 (MQ=255)
         gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                               >  2:1354604/1‑100 (MQ=255)
         gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                               >  2:1714523/1‑100 (MQ=255)
           tGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGcc                                                                             <  2:1164763/100‑1 (MQ=255)
            ggTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCa                                                                            >  1:296958/1‑100 (MQ=255)
               ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATAc                                                                         <  1:1476693/100‑1 (MQ=255)
               ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATAc                                                                         <  2:1679111/100‑1 (MQ=255)
               ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATAc                                                                         <  2:827407/100‑1 (MQ=255)
                tcgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGTGGATCGGATCTGGCTCATCGCCATAc                                                                         <  1:1654386/99‑1 (MQ=255)
                     gAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaaa                                                                   <  2:1664010/100‑1 (MQ=255)
                     gAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaaa                                                                   <  2:1477585/100‑1 (MQ=255)
                         gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                               <  2:1218618/100‑1 (MQ=255)
                            aGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                             >  1:975010/1‑99 (MQ=255)
                            aGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                                            >  1:831924/1‑100 (MQ=255)
                             gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGcc                                                           >  2:315855/1‑100 (MQ=255)
                              aaTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGccc                                                          >  1:1911976/1‑100 (MQ=255)
                               aTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCg                                                         >  1:53128/1‑100 (MQ=255)
                                tGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCgg                                                        >  2:890022/1‑100 (MQ=255)
                                  tCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCg                                                      >  2:898013/1‑100 (MQ=255)
                                   cGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                     <  2:1498514/100‑1 (MQ=255)
                                    gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                    <  1:332883/100‑1 (MQ=255)
                                      tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggt                                                  >  1:1196875/1‑100 (MQ=255)
                                       gATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                                 >  2:161987/1‑100 (MQ=255)
                                           ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                             >  2:938238/1‑100 (MQ=255)
                                               aataataatAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATc                                         >  2:1887764/1‑100 (MQ=255)
                                                  aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAATGCGCGCCCGGCGTGGTGCCGATATCGGa                                      <  1:1125372/100‑1 (MQ=255)
                                                  aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCAGCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                       <  2:1369231/99‑1 (MQ=255)
                                                   ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                     >  1:809732/1‑100 (MQ=255)
                                                    taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    >  2:716850/1‑100 (MQ=255)
                                                    taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    <  1:415242/100‑1 (MQ=255)
                                                    taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    <  2:1571039/100‑1 (MQ=255)
                                                       taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                 <  2:853305/100‑1 (MQ=255)
                                                        aaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGAt                                <  2:447449/100‑1 (MQ=255)
                                                         aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATc                               >  2:1900520/1‑100 (MQ=255)
                                                            atgatgCGGGGCGACAGTTTGCAGATTGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccn                            <  2:330083/100‑2 (MQ=255)
                                                            atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                            >  1:1772229/1‑100 (MQ=255)
                                                             tccgatctggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                           <  2:786802/92‑1 (MQ=255)
                                                                  cGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgt                      >  1:1476441/1‑100 (MQ=255)
                                                                    gggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCCTACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGc                    <  1:1242203/100‑1 (MQ=255)
                                                                     ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGagatcgga                   >  1:786800/1‑93 (MQ=255)
                                                                     ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                   <  1:880403/100‑1 (MQ=255)
                                                                     ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                   <  2:831925/100‑1 (MQ=255)
                                                                        gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt                >  2:1877844/1‑100 (MQ=255)
                                                                        gACAGTTTGCAGATGGGGGATCGGATCTGGCTCAGCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt                <  1:1079093/100‑1 (MQ=255)
                                                                        gACAGTTTGCAGATGGGGGAGCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt                >  1:1062731/1‑100 (MQ=255)
                                                                         aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa               >  1:1483304/1‑100 (MQ=255)
                                                                         aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa               <  2:101252/100‑1 (MQ=255)
                                                                         aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACCGATCCCCGACTGTGCGTGTa               >  1:1393375/1‑100 (MQ=255)
                                                                          cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa              >  2:483442/1‑100 (MQ=255)
                                                                            gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa              <  2:1867902/98‑1 (MQ=255)
                                                                            gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa            >  2:1458351/1‑100 (MQ=255)
                                                                            gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa            <  1:1811359/100‑1 (MQ=255)
                                                                              ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATg          >  1:1630384/1‑100 (MQ=255)
                                                                               tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg         <  2:1468872/100‑1 (MQ=255)
                                                                               tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg         >  1:256564/1‑100 (MQ=255)
                                                                               tGCAGATGGAGGATAGGATCGGGCTCAACGCCATACGGGAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGGGCGTGTAACATgg         <  1:595823/100‑1 (MQ=255)
                                                                                gCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgc        >  2:103945/1‑100 (MQ=255)
                                                                                gCAGATGGGGGATCGGATCTGGCTCATCGCCATATGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgc        >  1:82554/1‑100 (MQ=255)
                                                                                 cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGGTATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg       <  2:1246466/100‑1 (MQ=255)
                                                                                 cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg       <  2:53128/100‑1 (MQ=255)
                                                                                 cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg       <  2:411185/100‑1 (MQ=255)
                                                                                 cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg       >  2:108949/1‑100 (MQ=255)
                                                                                 cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg       <  2:940988/100‑1 (MQ=255)
                                                                                   gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg     <  1:275986/100‑1 (MQ=255)
                                                                                   gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg     >  2:830994/1‑100 (MQ=255)
                                                                                      ggggGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGttt  >  1:405899/1‑100 (MQ=255)
                                                                                         |                                                                                                
GTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTT  >  NC_000913/4615770‑4615955

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