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%	123.4 / NA	36	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/12);  total (24/12)

TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTG  >  NC_000913/1905699‑1905817                                                               |                                                         tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATg                                                        >  1:1665148/1‑65 (MQ=255)  gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                       >  1:1135186/1‑65 (MQ=255)     gcTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                    <  1:3794973/64‑1 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                    >  1:699397/1‑65 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                    >  1:1542718/1‑65 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATAGGTAGAGATGCAAc                                                    >  1:1978881/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                  >  1:358265/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                  >  1:359100/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                  >  1:2684960/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               >  1:1698422/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               <  1:1988000/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               <  1:639922/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               <  1:2651453/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               <  1:3614340/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               >  1:2821498/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                               >  1:1125530/1‑65 (MQ=255)                        gTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAttt                                 >  1:2059716/1‑65 (MQ=255)                           tGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCt                              >  1:120492/1‑65 (MQ=255)                              cGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGaa                           >  1:2863301/1‑65 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                         <  1:856888/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                         <  1:2271383/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                         <  1:1004750/65‑1 (MQ=255)                                 tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCa                        >  1:3162720/1‑65 (MQ=255)                                 tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCa                        >  1:3321217/1‑65 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                      >  1:422795/1‑65 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                      >  1:3550563/1‑65 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                      <  1:2953604/65‑1 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGGAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                      >  1:2676672/1‑65 (MQ=255)                                       cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGa                  <  1:2350228/65‑1 (MQ=255)                                       cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGa                  >  1:3724923/1‑65 (MQ=255)                                       cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGa                  >  1:751615/1‑65 (MQ=255)                                             ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc            <  1:1064195/65‑1 (MQ=255)                                             ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc            <  1:367079/65‑1 (MQ=255)                                             ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc            <  1:1459294/65‑1 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc           >  1:552195/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc           >  1:566131/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc           >  1:2114868/1‑65 (MQ=255)                                                 aTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTa        <  1:3023587/65‑1 (MQ=255)                                                      tCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttttt   >  1:622679/1‑65 (MQ=255)                                                       cGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTg  >  1:935831/1‑65 (MQ=255)                                                               |                                                         TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTG  >  NC_000913/1905699‑1905817

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