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

ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859404‑2859670                                                                                                                                         |                                                                                                                                     aCGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTc                                                                                                                                  <  1:18371/139‑1 (MQ=255)     ccTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                <  2:152392/139‑1 (MQ=255)       tGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTAcc                                                                                                                            <  2:1606933/139‑1 (MQ=255)           ggTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaa                                                                                                                        <  1:943276/139‑1 (MQ=255)           ggTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCagaa                                                                                                                        <  2:602493/139‑1 (MQ=255)                        gCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTATGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTg                                                                                                           <  2:815546/139‑1 (MQ=255)                              gTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                     <  1:1280072/139‑1 (MQ=255)                                  cAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGc                                                                                                 <  2:32814/139‑1 (MQ=255)                                      agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTg                                                                                             >  2:928786/1‑139 (MQ=255)                                        agGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGa                                                                                           <  2:29264/139‑1 (MQ=255)                                          gTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGaaa                                                                                         <  2:1139183/139‑1 (MQ=255)                                                  gcgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 >  2:569275/1‑139 (MQ=255)                                                  gcgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGAt                                                                                 <  1:815796/139‑1 (MQ=255)                                                   cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                <  1:1537364/123‑1 (MQ=255)                                                   cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGct                                                                                                >  2:1537364/1‑123 (MQ=255)                                                   cgcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATc                                                                                >  2:1027089/1‑139 (MQ=255)                                                    gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                   <  2:911464/87‑1 (MQ=255)                                                    gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTtgttgt                                                                                                                                   >  1:911464/1‑87 (MQ=255)                                                    gcgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATcc                                                                               >  1:1531219/1‑139 (MQ=255)                                                     cgcACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                                              >  1:111012/1‑139 (MQ=255)                                                                          gCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAg                                                         >  2:961905/1‑139 (MQ=255)                                                                           cTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                                                             <  2:238763/134‑1 (MQ=255)                                                                           cTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCcgcg                                                             >  1:238763/1‑134 (MQ=255)                                                                             ggAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGgcg                                                      <  1:1525534/139‑1 (MQ=255)                                                                                aaGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTg                                                   <  2:1531219/139‑1 (MQ=255)                                                                                  gCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGa                                                 <  1:1531632/139‑1 (MQ=255)                                                                                    aTTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTAc                                                                                                                             <  1:876209/61‑1 (MQ=255)                                                                                    aTTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTATTTGTTGTCAGTAc                                                                                                                             >  2:876209/1‑61 (MQ=255)                                                                                      ttCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTgg                                             >  1:182753/1‑139 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                <  1:805798/54‑1 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAg                                                                                                                                >  2:805798/1‑54 (MQ=255)                                                                                        cGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAt                                           >  1:1351621/1‑139 (MQ=255)                                                                                          ccGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         <  1:1027089/139‑1 (MQ=255)                                                                                                    ccgcTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa                               <  2:111012/139‑1 (MQ=255)                                                                                                    ccgcTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGaa                               <  1:928786/139‑1 (MQ=255)                                                                                                         aCGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATTGCTTGAAGAGcc                          <  1:1240794/139‑1 (MQ=255)                                                                                                         aCGCAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGcc                          <  1:911337/139‑1 (MQ=255)                                                                                                            cAA‑‑GTGGATGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTgg                       <  2:1351621/139‑1 (MQ=255)                                                                                                                      tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                         <  1:243737/129‑1 (MQ=255)                                                                                                                      tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCt                         >  2:243737/1‑129 (MQ=255)                                                                                                                      tGGTACGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGtaataa               >  1:1320243/1‑139 (MQ=255)                                                                                                                           cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATTTTGATTCGGATTCACTCCACCCGCGTCCGGCGGTGGAGGGGATTTATCGCTTGAAGAGCCTGGTGGAATAACaaat          >  1:195460/1‑137 (MQ=255)                                                                                                                           cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGTTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAAtt          >  1:1045277/1‑139 (MQ=255)                                                                                                                           cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGtaataa               >  1:678584/1‑134 (MQ=255)                                                                                                                           cGCAAATGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGtaataa               <  2:678584/134‑1 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                     <  2:891438/43‑1 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCg                                                                                                     >  1:891438/1‑43 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTACTAACAATTCCCGa     <  2:517055/139‑1 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     <  1:316097/139‑1 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     >  2:1169602/1‑139 (MQ=255)                                                                                                                                aTGTCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     <  2:195460/139‑1 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                                              <  1:1454837/63‑1 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCg                                                                              >  2:1454837/1‑63 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         >  2:835541/1‑100 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGAttt                                         <  1:835541/100‑1 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     <  2:311424/136‑1 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa     >  1:311424/1‑136 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg  >  1:185695/1‑139 (MQ=255)                                                                                                                                   tCTTTGTTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg  >  2:1320116/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                     ACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAA‑‑GTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAG  >  NC_000913/2859404‑2859670

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