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

AAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCT  >  NC_000913/2970217‑2970478                                                                                                                                       |                                                                                                                                aaaGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCAt                                                                                                                             >  2:541202/1‑139 (MQ=255)                     gcTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                                         <  1:541202/139‑1 (MQ=255)                         aCCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                     <  2:171021/139‑1 (MQ=255)                          ccGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGc                                                                                                    <  1:218769/139‑1 (MQ=255)                          ccGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGc                                                                                                    <  1:812414/139‑1 (MQ=255)                            gAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                                  <  2:1209402/139‑1 (MQ=255)                             aaGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCata                                                                                                 >  2:1316182/1‑139 (MQ=255)                                          ccGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTc                                                                                    <  1:1380894/139‑1 (MQ=255)                                             gggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGc                                                                                                                >  2:943990/1‑108 (MQ=255)                                             gggCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGc                                                                                                                <  1:943990/108‑1 (MQ=255)                                                cGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                                         >  2:51271/1‑112 (MQ=255)                                                cGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTc                                                                                                         <  1:51271/112‑1 (MQ=255)                                                      gATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                        <  2:200271/139‑1 (MQ=255)                                                                           cTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGAc                                                                                                                    >  1:1065196/1‑74 (MQ=255)                                                                           cTTTTATTTGAAGAAGAATTTCACCAGTGAACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGAc                                                                                                                    <  2:1065196/74‑1 (MQ=255)                                                                             tttATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCa                                                 <  1:1316182/139‑1 (MQ=255)                                                                                     aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                                     >  1:631992/1‑111 (MQ=255)                                                                                     aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAggg                                                                     <  2:631992/111‑1 (MQ=255)                                                                                     aagaagaaTTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAAc                                         >  2:1444872/1‑139 (MQ=255)                                                                                                     gTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                                  >  1:807897/1‑66 (MQ=255)                                                                                                     gTGCACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCat                                                                                                  <  2:807897/66‑1 (MQ=255)                                                                                                        cACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                        >  2:1460626/1‑89 (MQ=255)                                                                                                        cACTACTGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCa                                                                        <  1:1460626/89‑1 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                                   >  1:1005143/1‑57 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                                   <  2:384688/57‑1 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                                   >  1:384688/1‑57 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCa                                                                                                   <  2:1005143/57‑1 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGc                                                                         >  1:152433/1‑83 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGc                                                                         <  2:152433/83‑1 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAg                                                             <  2:749765/95‑1 (MQ=255)                                                                                                             ctGGCCCGTCATTTTCTGATCCAGTAAACATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAg                                                             >  1:749765/1‑95 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGcc                                                                                                               >  1:1068319/1‑42 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGcc                                                                                                               <  2:1068319/42‑1 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCttt                                                                                                          <  1:1087094/47‑1 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCttt                                                                                                          >  2:1087094/1‑47 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                     <  1:1129607/52‑1 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCggg                                                                                                     >  2:1129607/1‑52 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGc                        <  2:367223/129‑1 (MQ=255)                                                                                                                gCCCGTCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGc                        >  1:367223/1‑129 (MQ=255)                                                                                                                     tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCAttttt                                                                              <  1:928493/70‑1 (MQ=255)                                                                                                                     tCATTTTCTGATCCAGTAACCATCGCTTTGACCTGCCGCTTCCCGGGCATACAATTACCGCTTCAttttt                                                                              >  2:928493/1‑70 (MQ=255)                                                                                                                            ctaaTCCAGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCTCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacgct  >  2:873442/4‑139 (MQ=255)                                                                                                                                   aGTAACCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTg                                                   <  2:1026830/83‑1 (MQ=255)                                                                                                                                       |                                                                                                                                AAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCT  >  NC_000913/2970217‑2970478

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