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

TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGG  >  NC_000913/1905699‑1905818                                                               |                                                          tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATg                                                         >  1:1006327/1‑65 (MQ=255)  gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                        <  1:4494390/65‑1 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                     <  1:823551/65‑1 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                     >  1:1195976/1‑65 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                     <  1:4387423/65‑1 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                     >  1:2259614/1‑65 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTAGGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                     <  1:3859589/65‑1 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCCTCAATCGGTAAAGATGCCAccc                                                   >  1:2967521/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                   >  1:4247607/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                   >  1:3792884/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                   >  1:2254971/1‑65 (MQ=255)         ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCACTCGGTCAAGATGACACCCTc                                                 >  1:4383859/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                >  1:1034852/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                >  1:2929129/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                >  1:2862790/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                <  1:2438203/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                >  1:1858530/1‑65 (MQ=255)             tcttGCGTTTGGTATGTCGATGGATGCAATTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                             >  1:781223/1‑65 (MQ=255)                         tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAtttt                                 >  1:3368740/1‑65 (MQ=255)                           tGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCt                               >  1:129254/1‑65 (MQ=255)                              cGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGaa                            >  1:2627858/1‑65 (MQ=255)                              cGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGaa                            >  1:4519513/1‑65 (MQ=255)                               gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAg                           >  1:1169754/1‑65 (MQ=255)                               gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAg                           >  1:2033813/1‑65 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                          <  1:1712260/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                          <  1:4347728/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                          <  1:3983895/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                          <  1:1517993/65‑1 (MQ=255)                                 tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCa                         >  1:1412868/1‑65 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                       >  1:2038734/1‑65 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                       <  1:1111436/65‑1 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                       >  1:994258/1‑65 (MQ=255)                                     tgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGc                     >  1:830732/1‑65 (MQ=255)                                       cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGa                   >  1:686899/1‑65 (MQ=255)                                         tttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                 >  1:1440280/1‑65 (MQ=255)                                          ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                >  1:2910767/1‑65 (MQ=255)                                           tgctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCg               >  1:2050516/1‑65 (MQ=255)                                            gctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg              <  1:2546075/65‑1 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc            >  1:4258741/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc            >  1:1510970/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc            >  1:3827853/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc            >  1:3572484/1‑65 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc            >  1:2297240/1‑65 (MQ=255)                                                  tCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt        <  1:3226465/65‑1 (MQ=255)                                                  tCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt        >  1:2270603/1‑65 (MQ=255)                                                        ggTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTgg  >  1:2081109/1‑65 (MQ=255)                                                               |                                                          TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGG  >  NC_000913/1905699‑1905818

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