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

TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGC  >  NC_000913/1905699‑1905821                                                               |                                                             tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATg                                                            <  1:1273579/65‑1 (MQ=255) tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATg                                                            >  1:795325/1‑65 (MQ=255) tGCTACTGTTCTTCTTGCGTTTGGTATGTCAATGGATGCATTTGCTGCATCAATCGGTAAAGATg                                                            >  1:3511794/1‑65 (MQ=255)  gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                           >  1:2132456/1‑65 (MQ=255)  gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                           >  1:1902676/1‑65 (MQ=255)   cTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCGTTTGCTGCATCAATCGGTAAAGATGcc                                                          <  1:901182/65‑1 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGGTGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                        >  1:2730035/1‑65 (MQ=255)     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGGATCAATCGGTAAAGATGCCAc                                                        >  1:3844845/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                      <  1:441512/65‑1 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                      >  1:1284339/1‑65 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                      <  1:3214300/65‑1 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                      <  1:2954288/65‑1 (MQ=255)       tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                      >  1:2668413/1‑65 (MQ=255)         ttcttcTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACtctc                                                    >  1:1754167/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCCCCCTcc                                                   <  1:2464136/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACTCTcc                                                   >  1:1111975/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                   >  1:3647658/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                   >  1:2957083/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                   <  1:978116/65‑1 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                   >  1:1550804/1‑65 (MQ=255)          tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                   <  1:1673066/65‑1 (MQ=255)            ttcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCAt                                                 >  1:953472/1‑65 (MQ=255)            ttcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCAt                                                 >  1:3648218/1‑65 (MQ=255)            ttcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCAt                                                 >  1:175601/1‑65 (MQ=255)             tcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                >  1:1029906/1‑65 (MQ=255)               ttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATaaa                                              >  1:506736/1‑65 (MQ=255)               ttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATaaa                                              >  1:1687738/1‑65 (MQ=255)                 gCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAcc                                            <  1:3864458/65‑1 (MQ=255)                        gTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAttt                                     >  1:3343266/1‑65 (MQ=255)                         tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAtttt                                    >  1:557495/1‑65 (MQ=255)                           tGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCt                                  <  1:2539176/65‑1 (MQ=255)                             tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGa                                >  1:2462916/1‑65 (MQ=255)                             tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGa                                <  1:1312877/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                             <  1:2472053/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                             <  1:361760/65‑1 (MQ=255)                                atggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                             <  1:3887027/65‑1 (MQ=255)                                 tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCa                            >  1:2588487/1‑65 (MQ=255)                                  ggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAt                           <  1:3226035/65‑1 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                          <  1:283032/65‑1 (MQ=255)                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACACTCCATAAACCGAAATTTTCTGAAGCAtt                          >  1:3279633/1‑65 (MQ=255)                                   gatgCATTGGCTGCATCAATGGGTAAAGATGCCTCCCTCCATAAACCGAAAGTTTCTGAAGCAtt                          <  1:2186096/65‑1 (MQ=255)                                      gCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                       >  1:16652/1‑65 (MQ=255)                                         tttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                    >  1:3826459/1‑65 (MQ=255)                                         tttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                    >  1:596406/1‑65 (MQ=255)                                          ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                   >  1:26625/1‑65 (MQ=255)                                            gctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg                 >  1:499983/1‑65 (MQ=255)                                            gctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg                 >  1:1648199/1‑65 (MQ=255)                                            gctgcATCAATCGGTAAAGATGCCACCCTCACTAAACCGAACTTTTCTGAAGCATTGCGAACCgg                 >  1:48886/1‑65 (MQ=255)                                             ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                <  1:1575598/65‑1 (MQ=255)                                             ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                <  1:200534/65‑1 (MQ=255)                                             ctacATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                <  1:3110789/62‑1 (MQ=255)                                             agacATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                <  1:1628263/62‑1 (MQ=255)                                              tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc               >  1:2646955/1‑65 (MQ=255)                                               gcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCt              <  1:252564/65‑1 (MQ=255)                                                 aTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTa            <  1:2408620/65‑1 (MQ=255)                                                      tCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttttt       >  1:499541/1‑65 (MQ=255)                                                      tCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttttt       >  1:1056140/1‑65 (MQ=255)                                                           aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGc  >  1:686918/1‑65 (MQ=255)                                                               |                                                             TGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGC  >  NC_000913/1905699‑1905821

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