Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,905,761	G→A	G25D (GGT→GAT)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,905,761	0	G	A	100.0%	150.6 / NA	43	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/18);  total (25/18)

TACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAT  >  NC_000913/1905702‑1905812                                                            |                                                    tACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCa                                                <  1:3620567/65‑1 (MQ=255)  aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGccct                                               >  1:3783707/1‑63 (MQ=255)  aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                               >  1:2831633/1‑65 (MQ=255)  aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                               <  1:1986148/65‑1 (MQ=255)  aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                               >  1:1823882/1‑65 (MQ=255)   cTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAcc                                              <  1:2390989/65‑1 (MQ=255)    tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                             >  1:3663322/1‑65 (MQ=255)    tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                             >  1:979727/1‑65 (MQ=255)      ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTc                                           <  1:1004719/65‑1 (MQ=255)      ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTc                                           <  1:465427/65‑1 (MQ=255)      ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTc                                           <  1:888062/65‑1 (MQ=255)      ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTc                                           <  1:808129/65‑1 (MQ=255)       tcttctTGCGTTTGGTCTGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          <  1:88177/65‑1 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          <  1:82896/65‑1 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          <  1:3679021/65‑1 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          >  1:3844468/1‑65 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          >  1:3907191/1‑65 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          >  1:1344120/1‑65 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          >  1:4037678/1‑65 (MQ=255)       tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                          >  1:1893726/1‑65 (MQ=255)                      tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAtttt                           >  1:1800751/1‑65 (MQ=255)                         gTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                        >  1:2266430/1‑65 (MQ=255)                           cGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGaa                      >  1:2118246/1‑65 (MQ=255)                            gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAg                     >  1:3374654/1‑65 (MQ=255)                            gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAg                     >  1:3708064/1‑65 (MQ=255)                             atggatgCGTTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                    <  1:4270702/65‑1 (MQ=255)                             atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                    <  1:1338003/65‑1 (MQ=255)                             agggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                    <  1:2733065/63‑1 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 <  1:688272/65‑1 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:769762/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:52611/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:1753518/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:1219518/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:4135141/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:4077773/1‑65 (MQ=255)                                gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                 >  1:1358770/1‑65 (MQ=255)                                       ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc          >  1:4256053/1‑65 (MQ=255)                                          ctgcATCATTCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc       <  1:610957/65‑1 (MQ=255)                                          ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc       <  1:3634614/65‑1 (MQ=255)                                          ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc       <  1:1450111/65‑1 (MQ=255)                                           tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACTGGcc      >  1:2540605/1‑65 (MQ=255)                                           tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc      >  1:1230851/1‑65 (MQ=255)                                           tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc      >  1:3680851/1‑65 (MQ=255)                                           tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc      >  1:2516679/1‑65 (MQ=255)                                              aTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTa   <  1:4050962/65‑1 (MQ=255)                                               tCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt  <  1:3212339/65‑1 (MQ=255)                                                            |                                                    TACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAT  >  NC_000913/1905702‑1905812

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