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%	164.3 / NA	43	T613P (ACC→CCC)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/21);  total (22/21)

ACCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATG  >  NC_000913/2858782‑2859076                                                                                                                                                    |                                                                                                                                                    aCCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATcccc                                                                                                                                                   <  2:130966/150‑1 (MQ=255) aCCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCg                                                                                                                                                  >  1:636183/1‑151 (MQ=255)      gAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCg                                                                                                                                             <  1:548285/151‑1 (MQ=255)        aCTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGaa                                                                                                                                           <  1:139290/151‑1 (MQ=255)        aCTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGaa                                                                                                                                           <  2:255961/151‑1 (MQ=255)            cacCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATg                                                                                                                                       <  2:172609/151‑1 (MQ=255)              ccTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATggg                                                                                                                                     >  1:343358/1‑151 (MQ=255)                tGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCg                                                                                                                                   >  1:437502/1‑151 (MQ=255)                    cGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTa                                                                                                                                >  1:43765/1‑150 (MQ=255)                     gACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTaa                                                                                                                               >  2:450774/1‑150 (MQ=255)                         ttCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGt                                                                                                                           >  1:118900/1‑150 (MQ=255)                          tCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGt                                                                                                                           <  1:215490/149‑1 (MQ=255)                            aTTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCt                                                                                                                       >  2:220442/1‑151 (MQ=255)                              tGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCtat                                                                                                                     >  2:88870/1‑151 (MQ=255)                                  aaaCCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATgc                                                                                                                 <  2:661086/151‑1 (MQ=255)                                     ccGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGcc                                                                                                              >  1:553327/1‑151 (MQ=255)                                      cGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCa                                                                                                             <  2:28928/151‑1 (MQ=255)                                          cATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGAcc                                                                                                         >  1:33323/1‑151 (MQ=255)                                               gCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACt                                                                                                    <  1:356939/151‑1 (MQ=255)                                                 aTTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTg                                                                                                   >  2:229732/1‑150 (MQ=255)                                                    aCCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTg                                                                                               >  1:282494/1‑151 (MQ=255)                                                        aaGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCt                                                                                           <  2:285424/151‑1 (MQ=255)                                                          ggTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGa                                                                                         <  2:658806/151‑1 (MQ=255)                                                            tCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATg                                                                                       >  2:250331/1‑151 (MQ=255)                                                             cGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATgg                                                                                      >  1:72603/1‑151 (MQ=255)                                                             cGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATgg                                                                                      >  2:112015/1‑151 (MQ=255)                                                                 aTCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCt                                                                                   >  2:251834/1‑150 (MQ=255)                                                                        aGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGc                                                                           >  2:333091/1‑151 (MQ=255)                                                                            gAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGc                                                                        <  2:644962/150‑1 (MQ=255)                                                                                    aCTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTAcc                                                                <  1:391340/150‑1 (MQ=255)                                                                                       gAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCgg                                                              <  2:39810/149‑1 (MQ=255)                                                                                            aGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaaa                                                       <  1:489202/151‑1 (MQ=255)                                                                                            aGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACaaaa                                                       <  1:361226/151‑1 (MQ=255)                                                                                              ccATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAg                                                     <  2:343369/151‑1 (MQ=255)                                                                                                 tttATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGCCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGACACAACAAgac                                                   >  1:57018/1‑148 (MQ=255)                                                                                                       gCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTg                                            <  2:540634/151‑1 (MQ=255)                                                                                                        ccAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTgg                                           >  1:321709/1‑151 (MQ=255)                                                                                                        ccAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTgg                                           <  2:67901/151‑1 (MQ=255)                                                                                                                  gAATCTGTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcga                                  <  2:218255/150‑1 (MQ=255)                                                                                                                  gAATCTGTCGCCGCAGAGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTatcg                                   >  1:384747/1‑149 (MQ=255)                                                                                                                        gTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGc                              <  1:546525/148‑1 (MQ=255)                                                                                                                        gTCGCCGCAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATc                           >  1:83627/1‑151 (MQ=255)                                                                                                                           gccgcAGCGCCGCATGTTGATCATCCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCttt                        <  1:346397/151‑1 (MQ=255)                                                                                                                                              atcatcCCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGc       >  1:623757/1‑149 (MQ=255)                                                                                                                                               tcatccCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGc       >  1:431463/1‑148 (MQ=255)                                                                                                                                                 atccCCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATg  >  1:546843/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACCCTGAACTACACCTGCCCGACCTTCATTGATAAACCGGGCATTCGCATTACCGAAGGTCGCCATCCGGTAGTTGAACAAGTACTGAATGAGCCATTTATCGCCAACCCGCTGAATCTGTCGCCGCAGCGCCGCATGTTGATCATCACCGGTCCGAACATGGGCGGTAAAAGTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATG  >  NC_000913/2858782‑2859076

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