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

GAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTG  >  NC_000913/2970216‑2970479                                                                                                                                        |                                                                                                                                 gAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCa                                                                                                                               >  2:344200/1‑139 (MQ=255)       agcGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCtt                                                                                                                         >  2:978975/1‑139 (MQ=255)         cGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTg                                                                                                                       >  2:927066/1‑139 (MQ=255)                        ttACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCg                                                                                                        <  1:639931/139‑1 (MQ=255)                         tACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCgg                                                                                                       >  2:752006/1‑139 (MQ=255)                         tACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCgg                                                                                                       <  2:14096/139‑1 (MQ=255)                          aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                      <  1:344200/139‑1 (MQ=255)                                     tGGTGCCCGGGGCGAACGGATTCTGACGCTGCAACGCGGATTTTATTTGAAGAAGAATTTCACCAGTGCACGACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATATTTAc                                                                                           <  2:915837/139‑1 (MQ=255)                                              gggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAtt                                                                                  >  2:1015599/1‑139 (MQ=255)                                               ggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAttt                                                                                 >  2:84228/1‑139 (MQ=255)                                                     cGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGTACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                                           <  1:522776/139‑1 (MQ=255)                                                     cGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                                           <  1:336298/139‑1 (MQ=255)                                                       gATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                         <  1:3474/139‑1 (MQ=255)                                                         ttCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGTCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAgg                                                                       <  1:1015599/139‑1 (MQ=255)                                                           cTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGc                                                                     <  1:84228/139‑1 (MQ=255)                                                                    ccACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGAt                                                            <  1:360957/139‑1 (MQ=255)                                                                    ccACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGAt                                                            >  2:596853/1‑139 (MQ=255)                                                                             ttttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGc                                                   <  1:835010/139‑1 (MQ=255)                                                                                tATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATg                                                <  1:713473/139‑1 (MQ=255)                                                                                                                 gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCgg               >  2:1325/1‑139 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAATAATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                                           >  2:350386/1‑74 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTgg                                                                           <  1:350386/74‑1 (MQ=255)                                                                                                                      tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgc          >  2:535563/1‑139 (MQ=255)                                                                                                                           ttCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacg     <  1:978975/139‑1 (MQ=255)                                                                                                                           ttCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacg     <  1:596853/139‑1 (MQ=255)                                                                                                                              tGATCCAGTAACCATCTCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacgctg  <  1:310641/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 GAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTG  >  NC_000913/2970216‑2970479

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