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

GGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGA  >  NC_000913/2859406‑2859667                                                                                                                                     |                                                                                                                                  ggCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                             <  1:1194581/139‑1 (MQ=255)                      gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTg                                                                                                        >  2:1313253/1‑139 (MQ=255)                                        gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaaa                                                                                      >  1:885912/1‑139 (MQ=255)                                         ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                <  2:136577/96‑1 (MQ=255)                                         ttattaAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                >  1:136577/1‑96 (MQ=255)                                                gcgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                              >  1:369489/1‑139 (MQ=255)                                                 cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                  >  2:8068/1‑118 (MQ=255)                                                 cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                  <  1:8068/118‑1 (MQ=255)                                                     cACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGa                                                                         <  1:125040/139‑1 (MQ=255)                                                          cAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCCTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATtcac                                                                    >  1:484933/1‑139 (MQ=255)                                                                ctgtgTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGa                                                                                        >  1:303711/5‑113 (MQ=255)                                                                ctgtgTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGa                                                                                        <  2:303711/109‑1 (MQ=255)                                                                    gTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTa                                                                                                                           >  2:120972/1‑74 (MQ=255)                                                                    gTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTa                                                                                                                           <  1:120972/74‑1 (MQ=255)                                                                        gCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTAcc                                                                                                                         <  2:558894/72‑1 (MQ=255)                                                                        gCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTAcc                                                                                                                         >  1:558894/1‑72 (MQ=255)                                                                           ggAAAGCATTTCGCCGAACGCCGCCGCGACGCAAGTGGATGGTACTCAAAGGTGTTTGTTGTCACTACCAGAAGAAACGTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCCCCCCGCGTCAGgcg                                                   <  1:1097940/139‑1 (MQ=255)                                                                             aaaGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGGAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCt                                                 >  2:390005/1‑139 (MQ=255)                                                                                 cATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaa                                                                                       <  2:12439/97‑1 (MQ=255)                                                                                 cATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaa                                                                                       >  1:12439/1‑97 (MQ=255)                                                                                                     cTACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAagag                         <  2:1284358/139‑1 (MQ=255)                                                                                                      tACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                  >  1:1032922/1‑65 (MQ=255)                                                                                                      tACGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                  <  2:1032922/65‑1 (MQ=255)                                                                                                       aCGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc                       >  1:1165702/1‑139 (MQ=255)                                                                                                       aCGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc                       <  2:369489/139‑1 (MQ=255)                                                                                                        cGCAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGGAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      <  1:390005/139‑1 (MQ=255)                                                                                                          cAAGTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTgg                    <  1:2741/139‑1 (MQ=255)                                                                                                                  tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                      <  1:148568/113‑1 (MQ=255)                                                                                                                  tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                      >  2:148568/1‑113 (MQ=255)                                                                                                                          aaaTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTccc    >  1:1159541/1‑139 (MQ=255)                                                                                                                            aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa  <  1:1121737/139‑1 (MQ=255)                                                                                                                               tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                                                  <  1:340640/56‑1 (MQ=255)                                                                                                                               tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCt                                                                                  >  2:340640/1‑56 (MQ=255)                                                                                                                               tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      <  2:931201/116‑1 (MQ=255)                                                                                                                               tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCACGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                      >  1:931201/1‑116 (MQ=255)                                                                                                                                     |                                                                                                                                  GGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGA  >  NC_000913/2859406‑2859667

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