Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,859,176	T→G	I695S (ATC→AGC)	mutS →	DNA mismatch repair protein MutS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,859,176	0	T	G	100.0%	129.2 / NA	37	I695S (ATC→AGC)	mutS	DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/20);  total (17/20)

CCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAG  >  NC_000913/2859098‑2859252                                                                               |                                                                             ccTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCg                                                                            <  1:515295/81‑1 (MQ=255)  cTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCgg                                                                           <  1:515296/81‑1 (MQ=255)  cTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCgg                                                                           <  2:515297/81‑1 (MQ=255)     tATGGTGGAGATGACTGAAACCGCCAATATTTCACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCg                                                                        >  2:515366/1‑81 (MQ=255)      aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCg                                                                        <  2:515298/80‑1 (MQ=255)        ggtggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTg                                                                      >  2:515339/1‑80 (MQ=255)             agaTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGt                                                                <  2:515299/81‑1 (MQ=255)               aTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTcc                                                              <  2:515300/81‑1 (MQ=255)               aTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTc                                                               >  1:515351/1‑80 (MQ=255)                tGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCa                                                             <  2:515301/81‑1 (MQ=255)                              aaTATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTg                                               >  2:515356/1‑81 (MQ=255)                               atatTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGt                                              >  2:515357/1‑81 (MQ=255)                               atatTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGt                                              >  1:515350/1‑81 (MQ=255)                                          aaCGCCACCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgg                                   >  2:515352/1‑81 (MQ=255)                                               cACCAAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcgtg                              <  1:515302/81‑1 (MQ=255)                                                aCCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc                             <  2:515303/81‑1 (MQ=255)                                                aCCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc                             <  1:515304/81‑1 (MQ=255)                                                aCCGAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc                             <  1:515305/81‑1 (MQ=255)                                                   gAATACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGcgcg                          >  2:515364/1‑81 (MQ=255)                                                    aaTACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGcgcg                          >  2:515370/1‑80 (MQ=255)                                                    aaTACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCgg                         <  2:515306/81‑1 (MQ=255)                                                      tACAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGa                        >  2:515347/1‑80 (MQ=255)                                                       aCAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGCACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaa                       >  2:515371/1‑80 (MQ=255)                                                        cAGTCTGGTGTTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaaa                      <  2:515307/80‑1 (MQ=255)                                                               gtgtTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGc               <  1:515308/80‑1 (MQ=255)                                                               gtgtTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTCGc               <  2:515309/80‑1 (MQ=255)                                                                 gtTAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa            >  1:515362/1‑81 (MQ=255)                                                                  ttAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa            >  1:515361/1‑80 (MQ=255)                                                                   tAATGGATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAAt           >  2:515374/1‑80 (MQ=255)                                                                        gATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt      >  2:515359/1‑80 (MQ=255)                                                                        gATGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt      >  1:515369/1‑80 (MQ=255)                                                                         aTGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTa    <  1:515312/81‑1 (MQ=255)                                                                         aTGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTa    <  2:515311/81‑1 (MQ=255)                                                                         aTGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTa    <  2:515310/81‑1 (MQ=255)                                                                          tGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTaa   <  1:515313/81‑1 (MQ=255)                                                                          tGAGAGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTaa   <  2:515314/81‑1 (MQ=255)                                                                           gagaGCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAg  >  2:515377/1‑81 (MQ=255)                                                                               |                                                                             CCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAG  >  NC_000913/2859098‑2859252

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