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%	123.7 / NA	39	*230* (TAG→TAA)	mutH	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/22);  total (17/22)

GAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGG  >  NC_000913/2970216‑2970466                                                                                                                                        |                                                                                                                    gAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCa                                                                                                                  <  1:390961/139‑1 (MQ=255)         cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGTCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                          >  1:1555304/1‑139 (MQ=255)         cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                          >  1:857649/1‑139 (MQ=255)         cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                          >  1:29161/1‑139 (MQ=255)             tGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCt                                                                                                      >  1:327333/1‑139 (MQ=255)                        ttACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                           <  2:1141386/139‑1 (MQ=255)                             gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                      <  2:1470728/139‑1 (MQ=255)                             gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                      <  2:1299785/139‑1 (MQ=255)                             gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                      <  1:1006235/139‑1 (MQ=255)                                           ccGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                        >  1:1661088/1‑139 (MQ=255)                                               ggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAttt                                                                    >  1:961046/1‑139 (MQ=255)                                                  gAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAtttttt                                                                 >  1:893647/1‑139 (MQ=255)                                                  gAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAtttttt                                                                 >  1:1173884/1‑139 (MQ=255)                                                       gATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                            <  2:1661088/139‑1 (MQ=255)                                                                 cTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTa                                                  >  1:9865/1‑139 (MQ=255)                                                                  tGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAg                                                 <  2:29161/139‑1 (MQ=255)                                                                  tGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAg                                                 <  2:893647/139‑1 (MQ=255)                                                                           gCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                              >  1:420875/1‑117 (MQ=255)                                                                           gCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                              <  2:420875/117‑1 (MQ=255)                                                                             ttttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                      <  2:1555304/139‑1 (MQ=255)                                                                             ttttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                      <  2:857649/139‑1 (MQ=255)                                                                             ttttATTTGAAGAAGAATTCCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                      <  2:265760/139‑1 (MQ=255)                                                                                tATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATg                                   <  2:620151/139‑1 (MQ=255)                                                                                            aaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                       >  1:1638109/1‑75 (MQ=255)                                                                                            aaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                       <  2:1638109/75‑1 (MQ=255)                                                                                                       tGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTa                                                                               >  2:1367828/1‑72 (MQ=255)                                                                                                       tGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTa                                                                               <  1:1367828/72‑1 (MQ=255)                                                                                                            tactGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTgcg       <  2:961046/139‑1 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                           >  1:1465137/1‑50 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                           <  2:1465137/50‑1 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                      >  2:1177387/1‑103 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                      <  1:1177387/103‑1 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATc                        >  1:428950/1‑117 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATc                        <  2:428950/117‑1 (MQ=255)                                                                                                                 gCCCGTCATCTTCTGATCCAGTAACCATCGCTTTTACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCgg  <  2:1173884/139‑1 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCttttt                                                   >  2:293033/1‑85 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCttttt                                                   <  1:293033/85‑1 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGAt                                >  2:1030697/1‑104 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGAt                                <  1:1030697/104‑1 (MQ=255)                                                                                                                                        |                                                                                                                    GAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGG  >  NC_000913/2970216‑2970466

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