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

ACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGT  >  NC_000913/1905697‑1905824                                                                 |                                                                ccTGCTACTGTTCTTCTTGCGCTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                              >  1:79889/2‑68 (MQ=255)   tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCa                                                            >  1:2335514/1‑68 (MQ=255)    gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                           <  1:1496028/68‑1 (MQ=255)      tACTGTACTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAccc                                                         >  1:3591116/1‑68 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      >  1:3878801/1‑68 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      >  1:3607371/1‑68 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:100391/68‑1 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      >  1:3208851/1‑68 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:1983780/68‑1 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:2209335/68‑1 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGTTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:3198475/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:1569418/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:741219/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:2327371/1‑68 (MQ=255)               tcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                <  1:684736/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:216185/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                >  1:644324/1‑68 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:3568577/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:44086/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:148817/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                >  1:3876684/1‑68 (MQ=255)                                 gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAt                              <  1:1671972/68‑1 (MQ=255)                                   tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTg                            >  1:480251/1‑68 (MQ=255)                                   tggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTg                            >  1:1654929/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:606337/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:3821191/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:36507/1‑68 (MQ=255)                                      atgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGa                         <  1:358761/68‑1 (MQ=255)                                       tgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGaa                        <  1:532834/68‑1 (MQ=255)                                        gCATTTGCTGCATCAATCGGTGAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAc                       >  1:3333712/1‑68 (MQ=255)                                         cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                      >  1:3914286/1‑68 (MQ=255)                                           tttGCGGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg                    <  1:1709309/68‑1 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   >  1:791278/1‑68 (MQ=255)                                                tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTa               >  1:458210/1‑68 (MQ=255)                                                tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCGTa               >  1:3611627/1‑68 (MQ=255)                                                tgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCCTa               >  1:3610534/1‑68 (MQ=255)                                                          ggTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGc     >  1:528562/1‑68 (MQ=255)                                                            tAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCg   <  1:1603669/68‑1 (MQ=255)                                                            tAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCg   <  1:1413983/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:3738879/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:1437349/68‑1 (MQ=255)                                                                 |                                                                ACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGT  >  NC_000913/1905697‑1905824

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: ‑