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_000913	2,859,538	0	C	T	100.0%	145.3 / NA	45	L816L (CTG→TTG)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/24);  total (21/24)

ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859404‑2859670                                                                                                                                       |                                                                                                                                     aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                                  <  1:1011893/139‑1 (MQ=255) aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                                  <  1:1415279/139‑1 (MQ=255)         gCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCag                                                                                                                          <  1:1332328/139‑1 (MQ=255)                       ggCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGTTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCt                                                                                                            <  1:792375/139‑1 (MQ=255)                             cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCAAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTc                                                                                                      >  2:1342673/1‑139 (MQ=255)                                           ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaaaa                                                                                        >  2:757127/1‑139 (MQ=255)                                                       cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                                                            >  1:1469793/1‑139 (MQ=255)                                                       cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                                                            <  1:1394731/139‑1 (MQ=255)                                                        aCGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAt                                                                           >  1:943433/1‑139 (MQ=255)                                                            cAAAAGCTGCGTGAGCTAGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcac                                                                       >  2:601446/1‑139 (MQ=255)                                                                     cGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTAcc                                                                                                                            >  2:441505/1‑77 (MQ=255)                                                                     cGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTAcc                                                                                                                            <  1:441505/77‑1 (MQ=255)                                                                      gTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                                                             >  1:989107/1‑139 (MQ=255)                                                                        gAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                                                           >  2:748944/1‑139 (MQ=255)                                                                        gAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                                                           <  1:757127/139‑1 (MQ=255)                                                                             ggAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                      <  2:531632/139‑1 (MQ=255)                                                                             ggAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                      <  2:200922/139‑1 (MQ=255)                                                                               aaaGCATTTCGCCGAACGGCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGACGCt                                                    >  1:993326/1‑139 (MQ=255)                                                                                   cATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg                                                <  1:601446/139‑1 (MQ=255)                                                                                   cATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg                                                <  2:943433/139‑1 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                <  1:86299/54‑1 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                >  2:86299/1‑54 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                           >  2:388500/1‑139 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTCTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                           >  2:1066120/1‑139 (MQ=255)                                                                                                   gccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGa                                <  1:1066120/139‑1 (MQ=255)                                                                                                       cTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 >  2:718653/1‑86 (MQ=255)                                                                                                       cTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 <  1:718653/86‑1 (MQ=255)                                                                                                        tACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGc                           <  2:989107/139‑1 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc                          <  1:438/138‑1 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc                          >  2:438/1‑138 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                         <  2:158124/139‑1 (MQ=255)                                                                                                            cAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTgg                       <  1:1342673/139‑1 (MQ=255)                                                                                                               gTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                <  2:293410/79‑1 (MQ=255)                                                                                                               gTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                >  1:293410/1‑79 (MQ=255)                                                                                                                 ggATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGtaa                  <  1:748944/139‑1 (MQ=255)                                                                                                                  gATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGtaat                 >  2:22007/1‑139 (MQ=255)                                                                                                                            aaaTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTccc       >  2:548613/1‑139 (MQ=255)                                                                                                                              aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg                                                >  1:719841/1‑96 (MQ=255)                                                                                                                              aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAg                                                <  2:719841/96‑1 (MQ=255)                                                                                                                              aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     <  1:548613/139‑1 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         >  1:738963/1‑100 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         <  2:738963/100‑1 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     <  1:702552/136‑1 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     >  2:702552/1‑136 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg  >  2:1138598/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859404‑2859670

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