Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,557 G→T 100% G49V (GGT→GTT)  mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,5570GT98.3% 256.7 / ‑3.4 82G49V (GGT→GTT) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (42/38);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397461‑4397652
                                                                                                 |                                                                                               
gtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgg                                                                                              <  1:45622/100‑1 (MQ=255)
  gaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTgggg                                                                                            <  2:1694020/100‑1 (MQ=255)
    ggtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCg                                                                                          >  2:2079864/1‑100 (MQ=255)
       ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                       >  1:936780/1‑100 (MQ=255)
       ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGaaa                                                                                       >  2:643742/1‑100 (MQ=255)
         tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACt                                                                                     >  1:1053739/1‑100 (MQ=255)
         tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACt                                                                                     <  2:393635/100‑1 (MQ=255)
              cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATcc                                                                                <  2:611123/100‑1 (MQ=255)
               gACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCAACCCGTATCCATATTAATATCGAAGGCGTTGGGGCGAAACTTATACg                                                                               >  2:1515027/1‑100 (MQ=255)
                aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGc                                                                              >  2:1951496/1‑100 (MQ=255)
                aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGc                                                                              <  1:2002627/100‑1 (MQ=255)
                 ccTGNGTCGGTAGTCAAAGAACCAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCg                                                                               <  1:1839125/98‑1 (MQ=255)
                     cGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCg                                                                         >  2:2115346/1‑100 (MQ=255)
                       tCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTg                                                                       <  2:1205441/100‑1 (MQ=255)
                        cGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGa                                                                      >  1:904805/1‑100 (MQ=255)
                          gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     >  1:1352596/1‑99 (MQ=255)
                          gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATa                                                                    >  2:2022613/1‑100 (MQ=255)
                           tAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCATGTGCGACTCGTATCGATATTGATATCGAACGTGTTGGGGCGAAACTTATCCGCATTCGTGAt                                                                     <  2:770150/98‑1 (MQ=255)
                            aGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAAc                                                                  >  2:1176811/1‑100 (MQ=255)
                             gTCAAATAACTAGTTGAAAACAGCGTTGGTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 <  2:1946506/100‑1 (MQ=255)
                             gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 <  1:264387/100‑1 (MQ=255)
                             gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:1786581/1‑100 (MQ=255)
                             gTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACg                                                                 >  1:1396531/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1971290/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1500756/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:1393870/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  1:1067380/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:814117/1‑100 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               <  2:952445/100‑1 (MQ=255)
                               cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGc                                                               >  2:983645/1‑100 (MQ=255)
                                  aGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGc                                                            >  2:211060/1‑100 (MQ=255)
                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTTGTGATAACGGCTGCGGTa                                                        >  2:676470/1‑100 (MQ=255)
                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                        >  1:1047757/1‑100 (MQ=255)
                                      cTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTa                                                        >  2:223338/1‑100 (MQ=255)
                                        aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                                       >  1:366618/1‑99 (MQ=255)
                                        aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATc                                                      <  1:2072437/100‑1 (MQ=255)
                                         gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  2:1115380/1‑100 (MQ=255)
                                         gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTGGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     >  1:1050719/1‑100 (MQ=255)
                                          tGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaa                                                    <  1:12518/100‑1 (MQ=255)
                                           gtaaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   <  2:1356145/98‑1 (MQ=255)
                                           ggAAAACAGCCTCCATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGGGAAACTTAACCGCATTCGTGATAACGGCTGCGGTATCaaa                                                   >  2:1047679/1‑100 (MQ=255)
                                            gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCa                                                     <  2:988545/97‑1 (MQ=255)
                                             aaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                                 <  1:62924/100‑1 (MQ=255)
                                              aaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  2:1363596/100‑1 (MQ=255)
                                               aaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAg                                               <  1:584235/100‑1 (MQ=255)
                                                aCAGCCTCTATGCAGTTGCGAAGCGTATCGATATTGATATCGAACGCGTTGGGGAGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATAAAAAAAg                                               <  2:989181/99‑1 (MQ=255)
                                                aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGTGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaaa                                                <  1:983647/98‑1 (MQ=255)
                                                aCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                              >  1:1386313/1‑100 (MQ=255)
                                                  acccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            >  1:1521892/3‑100 (MQ=255)
                                                  aGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:1137000/100‑1 (MQ=255)
                                                   gCCTCGTTGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATg                                            <  1:1646276/99‑1 (MQ=255)
                                                   gCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGa                                           <  1:650877/100‑1 (MQ=255)
                                                          tGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcg                                    >  2:1633531/1‑100 (MQ=255)
                                                           gCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                   <  1:211060/100‑1 (MQ=255)
                                                               gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:496723/100‑1 (MQ=255)
                                                               gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:366619/100‑1 (MQ=255)
                                                               gTGCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggc                               <  2:236074/100‑1 (MQ=255)
                                                                cccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGGTAACGGCTGCGGTTTCAAAAAAGATCAgctcgagctgg                                >  1:154891/5‑98 (MQ=16)
                                                                 gCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             >  2:1102147/1‑100 (MQ=255)
                                                                 gCGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgc                             <  1:1259844/100‑1 (MQ=255)
                                                                 ccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                              >  1:1814954/4‑99 (MQ=37)
                                                                 ccaaCGCGTTCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGCCTGCGGTATCAAAAAAGATGCgctggcgctggcg                              >  2:1783408/4‑99 (MQ=25)
                                                                  cGACGCGTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgct                            <  1:588891/100‑1 (MQ=255)
                                                                     cgcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  1:1029933/100‑1 (MQ=255)
                                                                     cgcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  2:24759/100‑1 (MQ=255)
                                                                      gcgTATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         <  2:1324565/99‑1 (MQ=255)
                                                                      cgtgttCGATATCAATATCGATACGAGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATAAgctggcgctggcgctg                           >  2:1265577/6‑98 (MQ=255)
                                                                      cgcgttCGATATCAATATCGATACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggc                         >  1:1048975/6‑100 (MQ=37)
                                                                         tATCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt                     <  2:1603632/100‑1 (MQ=255)
                                                                          aTCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                    <  2:464299/100‑1 (MQ=255)
                                                                          aTCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTc                    <  1:1926240/100‑1 (MQ=255)
                                                                           tCGATATTGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                   >  1:1644257/1‑100 (MQ=255)
                                                                                  tGATATCGCACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:1499663/100‑1 (MQ=255)
                                                                                  tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  2:429800/100‑1 (MQ=255)
                                                                                  tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            <  1:1804140/100‑1 (MQ=255)
                                                                                  tGATATCGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacc            >  2:1259188/1‑100 (MQ=255)
                                                                                        cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  2:2149064/1‑100 (MQ=255)
                                                                                        cGAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaa      >  2:2121478/1‑100 (MQ=255)
                                                                                         gAACGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa     >  2:714547/1‑100 (MQ=255)
                                                                                          aaCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAAt    >  1:132122/1‑100 (MQ=255)
                                                                                          aaCGCGTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAAt    >  1:2079441/1‑100 (MQ=255)
                                                                                             cgcgTTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCg  <  2:120173/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCG  >  NC_000913/4397461‑4397652

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