Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,970,351	G→A	*230* (TAG→TAA)	mutH →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,970,351	0	G	A	100.0%	70.2 / NA	23	*230* (TAG→TAA)	mutH	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/12);  total (11/12)

ACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAG  >  NC_000913/2970241‑2970485                                                                                                               |                                                                                                                                       aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                            <  2:122172/139‑1 (MQ=255)                              gATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                               <  1:638347/139‑1 (MQ=255)                                                     tttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCa                                                        <  1:222096/139‑1 (MQ=255)                                                             aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                                           <  2:58378/96‑1 (MQ=255)                                                             aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                                           >  1:58378/1‑96 (MQ=255)                                                                          ccAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGcc                                   <  1:701640/139‑1 (MQ=255)                                                                                cACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGc                                                                           >  2:433694/1‑93 (MQ=255)                                                                                cACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGc                                                                           <  1:433694/93‑1 (MQ=255)                                                                                   tactGGCGCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTgcg                          <  2:20251/139‑1 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                            <  2:774659/52‑1 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                            >  1:774659/1‑52 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTa                                                                     <  2:234342/91‑1 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTa                                                                     >  1:234342/1‑91 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                         <  1:647785/103‑1 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                         >  2:647785/1‑103 (MQ=255)                                                                                        gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCgg                     >  1:174284/1‑139 (MQ=255)                                                                                            gTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTcc                                                                                                               <  1:783805/45‑1 (MQ=255)                                                                                            gTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTcc                                                                                                               >  2:783805/1‑45 (MQ=255)                                                                                             tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                            >  2:900386/1‑47 (MQ=255)                                                                                             tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                            <  1:900386/47‑1 (MQ=255)                                                                                             tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgc                >  1:305829/1‑139 (MQ=255)                                                                                                       aTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCACATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACgctgct      >  2:846909/1‑139 (MQ=255)                                                                                                           aGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGgag  >  2:666637/1‑139 (MQ=255)                                                                                                               |                                                                                                                                       ACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAG  >  NC_000913/2970241‑2970485

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