Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,859,538 C→T L816L (CTG→TTG)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,859,5380CT100.0% 190.8 / NA 60L816L (CTG→TTG) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/29);  total (31/29)

GGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859406‑2859670
                                                                                                                                    |                                                                                                                                    
ggCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                >  1:333983/1‑139 (MQ=255)
      gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCag                                                                                                                          >  2:1401889/1‑139 (MQ=255)
      gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCag                                                                                                                          >  2:1593274/1‑139 (MQ=255)
                          cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAgaagaa                                                                                                                     <  2:1661121/124‑1 (MQ=255)
                          cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAgaagaa                                                                                                                     >  1:1661121/1‑124 (MQ=255)
                            tGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGa                                                                                                    >  1:1349512/1‑139 (MQ=255)
                             gCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                                  >  2:822446/1‑108 (MQ=255)
                             gCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                                  <  1:822446/108‑1 (MQ=255)
                              ccAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                                                  >  1:1095794/1‑139 (MQ=255)
                               cAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGc                                                                                                 <  1:508134/139‑1 (MQ=255)
                                       gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGATAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaag                                                                                                                       >  2:1766635/1‑109 (MQ=255)
                                       gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGATAGCATTTCGCCGAACGCCGCCGCTACGAAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaag                                                                                                                       <  1:1766635/109‑1 (MQ=255)
                                       gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGcc                                                                                                             >  2:607803/1‑119 (MQ=255)
                                       gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGcc                                                                                                             <  1:607803/119‑1 (MQ=255)
                                         tattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAAt                                                                                       >  1:1164570/1‑139 (MQ=255)
                                              tgcgcgcACGGCAAAAGCTGCGTTAGCTGGAAAGCATTTCGCCGATCGCCGCCGCTACGCACGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAACCTTCGCCTCCGGTCGAAGCTCTGGAAAATCTTGa                                                                                  >  1:622496/2‑139 (MQ=255)
                                                cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                     <  2:992007/118‑1 (MQ=255)
                                                cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                     >  1:992007/1‑118 (MQ=255)
                                                cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                <  1:1668862/139‑1 (MQ=255)
                                                  cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                >  1:384621/1‑89 (MQ=255)
                                                  cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                <  2:384621/89‑1 (MQ=255)
                                                  cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                                              >  1:1492033/1‑139 (MQ=255)
                                                  cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                                              >  2:1788217/1‑139 (MQ=255)
                                                    cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTg                                                                                                           >  1:127344/1‑108 (MQ=255)
                                                    cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTg                                                                                                           <  2:127344/108‑1 (MQ=255)
                                                    cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                                                            <  2:1620291/139‑1 (MQ=255)
                                                    cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                                                            <  1:949191/139‑1 (MQ=255)
                                                     aCGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAt                                                                           >  2:1344102/1‑139 (MQ=255)
                                                          aaaaGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcact                                                                      <  2:333983/139‑1 (MQ=255)
                                                                 gCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaa                                                                                          >  1:556387/1‑112 (MQ=255)
                                                                 gCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaa                                                                                          <  2:556387/112‑1 (MQ=255)
                                                                                  tttCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTg                                              >  2:1061584/1‑139 (MQ=255)
                                                                                       ccGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         <  2:1349512/139‑1 (MQ=255)
                                                                                          aaCGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                <  1:1445061/49‑1 (MQ=255)
                                                                                          aaCGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                >  2:1445061/1‑49 (MQ=255)
                                                                                               cgccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTg                                 >  2:1217433/1‑139 (MQ=255)
                                                                                                gccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGa                                <  2:1492033/139‑1 (MQ=255)
                                                                                                 ccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa                               <  2:583717/139‑1 (MQ=255)
                                                                                                 ccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTCATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa                               <  1:1061584/139‑1 (MQ=255)
                                                                                                       cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                         <  2:1095794/139‑1 (MQ=255)
                                                                                                       cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                         <  1:1401889/139‑1 (MQ=255)
                                                                                                         cAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTgg                       <  2:1836274/139‑1 (MQ=255)
                                                                                                            gTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAgg                                                        <  1:850667/103‑1 (MQ=255)
                                                                                                            gTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAgg                                                        >  2:850667/1‑103 (MQ=255)
                                                                                                             tGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGCCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGta                   <  1:1217433/139‑1 (MQ=255)
                                                                                                                 tGGTACGCAAATGTCTTTGTTGTCAGTAGCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 <  1:1297206/73‑1 (MQ=255)
                                                                                                                 tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 >  2:1297206/1‑73 (MQ=255)
                                                                                                                    tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                >  2:1576490/1‑55 (MQ=255)
                                                                                                                    tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                <  1:1576490/55‑1 (MQ=255)
                                                                                                                    tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                <  2:452617/71‑1 (MQ=255)
                                                                                                                    tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                >  1:452617/1‑71 (MQ=255)
                                                                                                                    tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACaa            >  1:652107/1‑139 (MQ=255)
                                                                                                                          aaTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                      >  2:1470549/1‑91 (MQ=255)
                                                                                                                          aaTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                      <  1:1470549/91‑1 (MQ=255)
                                                                                                                           aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     >  1:1753783/1‑139 (MQ=255)
                                                                                                                              tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                >  1:1102711/1‑45 (MQ=255)
                                                                                                                              tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                <  2:1102711/45‑1 (MQ=255)
                                                                                                                              tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg  >  1:870860/1‑139 (MQ=255)
                                                                                                                                tttGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGa                                <  2:346483/107‑1 (MQ=255)
                                                                                                                                tttGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGa                                >  1:346483/1‑107 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859406‑2859670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑