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%	138.6 / NA	44	L816L (CTG→TTG)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/19);  total (25/19)

ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGA  >  NC_000913/2859404‑2859667                                                                                                                                       |                                                                                                                                  aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTGCGCAAGTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                               <  1:641862/139‑1 (MQ=255)  cGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                >  1:1718454/1‑137 (MQ=255)  cGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                <  2:1718454/137‑1 (MQ=255)      cTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTAc                                                                                                                          <  1:7321/139‑1 (MQ=255)        ggCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCa                                                                                                                        >  1:771769/1‑139 (MQ=255)                                 ccAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAg                                                                                               <  2:771769/139‑1 (MQ=255)                                        agGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGa                                                                                        <  1:1548694/139‑1 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaaga                                                                                                                   >  1:1360421/1‑110 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaaga                                                                                                                   <  2:1360421/110‑1 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGc                                                                                                       <  2:827316/122‑1 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGc                                                                                                       >  1:827316/1‑122 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaaa                                                                                      >  2:863521/1‑139 (MQ=255)                                                   cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                             <  2:717988/139‑1 (MQ=255)                                                    gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATcc                                                                            >  2:883702/1‑139 (MQ=255)                                                         cGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                                                    <  1:815715/110‑1 (MQ=255)                                                         cGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGt                                                                                                    >  2:815715/1‑110 (MQ=255)                                                                  cTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCAccc                                                              >  1:1714333/1‑139 (MQ=255)                                                                    gCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcg                                                            >  2:1191989/1‑139 (MQ=255)                                                                    gCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcg                                                            >  2:393755/1‑139 (MQ=255)                                                                     cGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgc                                                           >  1:394969/1‑139 (MQ=255)                                                                             ggAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                   <  1:883702/139‑1 (MQ=255)                                                                                      ttCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                                                        >  1:1335089/1‑125 (MQ=255)                                                                                      ttCGCCGAACGCCGCCGCTACGCAAGTGGATGGGACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTc                                                        <  2:1335089/125‑1 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCt                                                                                                         <  2:339727/74‑1 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCt                                                                                                         >  1:339727/1‑74 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                        >  1:1211892/1‑139 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                        >  1:1291915/1‑139 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                        >  2:1290974/1‑139 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                        >  2:1337900/1‑139 (MQ=255)                                                                                         gCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAtt                                       >  2:1212995/1‑139 (MQ=255)                                                                                            gAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTAt                                    <  1:388125/139‑1 (MQ=255)                                                                                              aCGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCg                                  >  1:1497293/1‑139 (MQ=255)                                                                                                 ccgccgcTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCtt                               <  2:1211892/139‑1 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      <  1:1191989/139‑1 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      <  2:1714333/139‑1 (MQ=255)                                                                                                          cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      <  1:906470/139‑1 (MQ=255)                                                                                                                tGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                                                       <  1:419641/84‑1 (MQ=255)                                                                                                                tGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGAtt                                                                       >  2:419641/1‑84 (MQ=255)                                                                                                                       tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                                                          >  1:568924/1‑90 (MQ=255)                                                                                                                       tACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                                                          <  2:568924/90‑1 (MQ=255)                                                                                                                         cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAAtt       >  1:1317886/1‑139 (MQ=255)                                                                                                                         cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAAtt       >  1:62964/1‑139 (MQ=255)                                                                                                                            aaaTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTccc    >  1:1147400/1‑139 (MQ=255)                                                                                                                              aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa  >  2:340342/1‑139 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                   >  2:1072994/1‑87 (MQ=255)                                                                                                                                 tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGACCCGGATTCACTCACCCCGCGTCAGgcg                                                   <  1:1072994/87‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGA  >  NC_000913/2859404‑2859667

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: ‑