Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,858,929	A→C	T613P (ACC→CCC)	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,858,929	0	A	C	100.0%	184.8 / NA	50	T613P (ACC→CCC)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/28);  total (22/28)

CTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGC  >  NC_000913/2858785‑2859071                                                                                                                                                 |                                                                                                                                               cTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTc                                                                                                                                          >  1:201028/1‑151 (MQ=255)   gAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTcc                                                                                                                                         <  1:413828/150‑1 (MQ=255)   gAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTc                                                                                                                                          <  1:19644/149‑1 (MQ=255)    aaCTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGa                                                                                                                                       <  1:624247/151‑1 (MQ=255)     aCTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGa                                                                                                                                       >  2:369202/1‑150 (MQ=255)        acacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAAc                                                                                                                                     <  1:29859/149‑1 (MQ=255)                ccGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTa                                                                                                                           <  1:663779/151‑1 (MQ=255)                            gATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCtat                                                                                                                <  1:609706/150‑1 (MQ=255)                                aaCCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATgcg                                                                                                           <  2:194834/151‑1 (MQ=255)                                        aTTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCg                                                                                                   <  1:616435/151‑1 (MQ=255)                                             cATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACt                                                                                               >  1:27478/1‑150 (MQ=255)                                               ttACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGAt                                                                                            >  2:100633/1‑151 (MQ=255)                                                tACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGAtt                                                                                           >  2:563857/1‑151 (MQ=255)                                                     aaGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCt                                                                                      <  1:554351/151‑1 (MQ=255)                                                                   gTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAACAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCg                                                                        <  2:185906/151‑1 (MQ=255)                                                                   gTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCg                                                                        <  1:492042/151‑1 (MQ=255)                                                                     aGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCgg                                                                       <  2:294975/150‑1 (MQ=255)                                                                     aGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGc                                                                      <  1:332147/151‑1 (MQ=255)                                                                       ttGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGc                                                                      >  1:176670/1‑149 (MQ=255)                                                                        tGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAg                                                                    <  2:458515/150‑1 (MQ=255)                                                                          aaCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCt                                                                  >  2:373904/1‑150 (MQ=255)                                                                          aaCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCt                                                                  >  2:245660/1‑150 (MQ=255)                                                                          aaCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTa                                                                 <  1:249070/151‑1 (MQ=255)                                                                           aCAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTAt                                                                >  1:395453/1‑151 (MQ=255)                                                                            cAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATg                                                               <  2:630741/151‑1 (MQ=255)                                                                                  cTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCgg                                                         >  2:621463/1‑151 (MQ=255)                                                                                       tGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGcaca                                                     >  2:344838/1‑150 (MQ=255)                                                                                        gAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaa                                                   <  1:676022/151‑1 (MQ=255)                                                                                         aGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaa                                                   >  2:321778/1‑150 (MQ=255)                                                                                           ccATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAg                                                <  2:263606/151‑1 (MQ=255)                                                                                            cATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCTTCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGt                                               <  2:211705/151‑1 (MQ=255)                                                                                               ttATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCg                                             <  2:474276/150‑1 (MQ=255)                                                                                                tATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCg                                             >  2:70200/1‑149 (MQ=255)                                                                                                   cGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGAt                                         <  2:277573/150‑1 (MQ=255)                                                                                                     ccAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTgg                                      <  1:435053/151‑1 (MQ=255)                                                                                                       aaCCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGAc                                    >  1:213281/1‑151 (MQ=255)                                                                                                         cccGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCt                                  >  2:412163/1‑151 (MQ=255)                                                                                                               gAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcga                             <  2:129070/150‑1 (MQ=255)                                                                                                               gAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcg                              <  1:560266/149‑1 (MQ=255)                                                                                                                  tCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGc                         >  1:372061/1‑151 (MQ=255)                                                                                                                   cTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGc                         >  2:457442/1‑150 (MQ=255)                                                                                                                     gTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCa                        <  2:201023/149‑1 (MQ=255)                                                                                                                      tCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCt                     >  2:281662/1‑151 (MQ=255)                                                                                                                      tCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCt                     >  1:140579/1‑151 (MQ=255)                                                                                                                       cgccgcAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCt                     <  1:136008/150‑1 (MQ=255)                                                                                                                               cgccgcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCcgcg            <  1:47651/151‑1 (MQ=255)                                                                                                                                gccgcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGt           >  1:380417/1‑151 (MQ=255)                                                                                                                                 ccgcATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGt           >  1:549400/1‑150 (MQ=255)                                                                                                                                     aTGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGgcg      <  1:155197/151‑1 (MQ=255)                                                                                                                                      tGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGgcgc     >  1:216340/1‑151 (MQ=255)                                                                                                                                             catccCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGc  >  1:46555/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                               CTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGC  >  NC_000913/2858785‑2859071

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