Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,674,851 T→C 26.5% H304R (CAT→CGT)  yphC ← putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,674,8510TC26.5% 117.2 / 48.6 79H304R (CAT→CGT) yphCputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base C (4/17);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 5.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CATAGGCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCAATGTA  >  NC_000913/2674753‑2674948
                                                                                                  |                                                                                                 
cATAGGCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATg                                                                                                  >  1:367945/1‑100 (MQ=255)
 aTAGGCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATgg                                                                                                 <  1:223049/100‑1 (MQ=255)
 aTAGGCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATgg                                                                                                 >  2:142246/1‑100 (MQ=255)
    ggCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAAc                                                                                              >  1:500495/1‑100 (MQ=255)
    ggCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAAc                                                                                              >  1:1371908/1‑100 (MQ=255)
      cATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAg                                                                                            >  2:1548405/1‑100 (MQ=255)
       aTCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGa                                                                                           >  2:557050/1‑100 (MQ=255)
        tCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGAc                                                                                          >  2:879322/1‑100 (MQ=255)
         cACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACt                                                                                         >  2:180192/1‑100 (MQ=255)
           cctgcctgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTgg                                                                                       <  1:89899/100‑1 (MQ=255)
           cctgcctgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTgg                                                                                       >  2:1587176/1‑100 (MQ=255)
             tgcctgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTc                                                                                     <  2:2059382/100‑1 (MQ=255)
              gcctgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCa                                                                                    >  2:1206221/1‑100 (MQ=255)
              gcctgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCa                                                                                    >  1:2135905/1‑100 (MQ=255)
                ctgTTCCAGCGAGAAGCGATGGGTAATGGTGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCAcc                                                                                  <  2:1890422/100‑1 (MQ=255)
                 tgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCAccc                                                                                 <  1:25597/100‑1 (MQ=255)
                 tgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCACATTTTTCCATACGGAACAGACTGGTCAccc                                                                                 <  2:912403/100‑1 (MQ=255)
                 tgTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTCCCAGTCCGTCAGACCATGGGCGCATTTTTCCATATGGAACAGACTGGacacgc                                                                                 >  1:977701/1‑98 (MQ=255)
                  gTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCa                                                                                >  2:395855/1‑100 (MQ=255)
                      cAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCTGAACATGGGCGCATTTTTCGTTATGGAACACACTGGTCACGCAgg                                                                              >  2:1881975/1‑98 (MQ=255)
                        gCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGcc                                                                          <  2:615515/100‑1 (MQ=255)
                             aaGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATaa                                                                     <  2:500496/100‑1 (MQ=255)
                               gCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGACAATc                                                                   >  1:1893772/1‑100 (MQ=255)
                                cGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATcc                                                                  >  1:138572/1‑100 (MQ=255)
                                    gggTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATccgccg                                                              <  1:1266179/100‑1 (MQ=255)
                                      gTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGtt                                                            >  1:1112875/1‑100 (MQ=255)
                                          tGGCGTTACTCGGCCAATGCTTCCCGTCCGTCAGATCATGGACGCTTTTTTCCATACGGAACAGACTGGTCACCCAGTAGCCGATAATCCGCCGTTGa                                                          <  2:1437654/98‑1 (MQ=255)
                                           ggCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCTTATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGACGTTGATgg                                                       >  2:1068836/1‑100 (MQ=255)
                                           ggCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATgg                                                       <  2:465473/100‑1 (MQ=255)
                                           ggCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATgg                                                       <  2:2237638/100‑1 (MQ=255)
                                            gCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGt                                                      >  2:1275571/1‑100 (MQ=255)
                                             cGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTg                                                     <  2:1800860/100‑1 (MQ=255)
                                             cGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTg                                                     <  2:1163844/100‑1 (MQ=255)
                                              gTTACGGGGCGACAGCTTCCAGTCCGTCAGATCGGGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGc                                                    <  1:2124818/100‑1 (MQ=255)
                                              gTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGc                                                    <  1:1941412/100‑1 (MQ=255)
                                               ttACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCa                                                   <  1:1420292/100‑1 (MQ=255)
                                               ttACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACACAGGAGCCGATAATCCGCCGTTGATGGTGCa                                                   <  1:1749784/100‑1 (MQ=255)
                                               ttACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCa                                                   <  1:1600322/100‑1 (MQ=255)
                                                tACGCGGCCACATCTTCCAGTTCGTCATATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATTCGCCGTTGATGGGGCAt                                                  <  2:1211013/100‑1 (MQ=255)
                                                tACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCAt                                                  >  2:1938892/1‑100 (MQ=255)
                                                 aCGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATc                                                 >  1:1996662/1‑100 (MQ=255)
                                                 aCGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATc                                                 <  1:1548416/100‑1 (MQ=255)
                                                   gcgGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCa                                                <  2:270734/99‑1 (MQ=255)
                                                    cgGCCACAGCTCCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGa                                              >  2:919842/1‑100 (MQ=255)
                                                     ggCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGa                                              >  1:896599/1‑99 (MQ=255)
                                                      gttacaGCTTCCAGTCCGTCAGATCATGTGCGTATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATc                                            <  2:2085455/97‑1 (MQ=255)
                                                      gCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATc                                            >  1:102608/1‑100 (MQ=255)
                                                           aGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGgcgc                                       <  1:483926/100‑1 (MQ=255)
                                                           aGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGgcgc                                       <  1:531357/100‑1 (MQ=255)
                                                             cTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTg                                     <  1:1592476/100‑1 (MQ=255)
                                                               tCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGAc                                   <  1:557049/100‑1 (MQ=255)
                                                                ccAGTCCGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGAcc                                  >  1:397230/1‑100 (MQ=255)
                                                                 cAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCt                                 >  2:445189/1‑100 (MQ=255)
                                                                   gTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCg                               <  2:1722928/100‑1 (MQ=255)
                                                                      cGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAAt                            <  1:1844694/100‑1 (MQ=255)
                                                                      cGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAAt                            <  1:138388/100‑1 (MQ=255)
                                                                      cGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAAt                            <  2:1797512/100‑1 (MQ=255)
                                                                      cGTCAGATCATGGGCGCATTTTTCCATACGGAACAGACGGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAAt                            <  1:227583/100‑1 (MQ=255)
                                                                      cGTCAGATCATGGACGCATTTTTTCATATGGTACAGACGGGTCACCCAGGAGCCGATAATACGCCGTTGATGGTGCATCAGATCGGCGCTGCCCTCGAAt                            >  2:1383228/1‑100 (MQ=255)
                                                                       gTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAAtt                           <  2:1510245/100‑1 (MQ=255)
                                                                        tCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTc                          <  2:1531324/100‑1 (MQ=255)
                                                                         cAGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTcc                         <  1:1225377/100‑1 (MQ=255)
                                                                          aGATCATGGTCTTTTTTTTCCATACGGAACAGACTGGTCACCCAGGAGCAGATAAGCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCa                        <  1:715696/100‑1 (MQ=255)
                                                                          aGATCATGGGCGCATTTTTCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCa                        >  2:1071041/1‑100 (MQ=255)
                                                                           ggtCATGGGCGCATTTTTCCATATGGATCAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCa                        <  2:1191648/97‑1 (MQ=255)
                                                                           gATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCAc                       >  1:660914/1‑100 (MQ=255)
                                                                           gATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCAc                       >  1:1531061/1‑100 (MQ=255)
                                                                           gATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCAc                       >  1:735308/1‑100 (MQ=255)
                                                                             tCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACtt                     >  1:1033923/1‑100 (MQ=255)
                                                                                  ggCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCg                >  2:1554611/1‑100 (MQ=255)
                                                                                   gcgcATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTcc                 <  1:2101401/98‑1 (MQ=255)
                                                                                   gcgcATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCgg               >  1:2097133/1‑100 (MQ=255)
                                                                                     gcATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGtt             <  2:971527/100‑1 (MQ=255)
                                                                                         ttttCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCAc         <  2:1889989/100‑1 (MQ=255)
                                                                                         ttttCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACCTTTCCGGTTTCAc         <  1:1524809/100‑1 (MQ=255)
                                                                                           ttCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCa       <  1:921323/100‑1 (MQ=255)
                                                                                           ttCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCa       <  1:391964/100‑1 (MQ=255)
                                                                                           ttCCATACGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCa       >  1:1849701/1‑100 (MQ=255)
                                                                                              cATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCAATg    >  1:1330450/1‑100 (MQ=255)
                                                                                                tatGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCAATGTa  <  1:234972/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
CATAGGCATCACCTGCCTGTTCCAGCGAGAAGCGATGGGTAATGGCGTTACGCGGCCACAGCTTCCAGTCCGTCAGATCATGGGCGCATTTTTCCATATGGAACAGACTGGTCACCCAGGAGCCGATAATCCGCCGTTGATGGTGCATCAGATCGGCGCTGACCTCGAATTCCACTTTTCCGGTTTCACCAATGTA  >  NC_000913/2674753‑2674948

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