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

ACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGT  >  NC_000913/1905697‑1905824                                                                 |                                                                aCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                              >  1:480606/1‑68 (MQ=255) aCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                              >  1:2537002/1‑68 (MQ=255) aCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGc                                                              >  1:8621/1‑68 (MQ=255)  cTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGcc                                                             <  1:418189/68‑1 (MQ=255)   tGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCa                                                            >  1:1201733/1‑68 (MQ=255)    gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                           <  1:2982788/68‑1 (MQ=255)    gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                           <  1:1267174/68‑1 (MQ=255)    gCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCAc                                                           <  1:3041970/68‑1 (MQ=255)     cTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATGAATCGGTAAAGATGCCAcc                                                          >  1:2221246/1‑68 (MQ=255)      tACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCCTTTGCTGCATCAATCGGTAAAGATGCCAccc                                                         >  1:1280859/1‑68 (MQ=255)       aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCt                                                        >  1:2449039/1‑68 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:1396660/68‑1 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:3336361/68‑1 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:2315271/68‑1 (MQ=255)         tGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                      <  1:889493/68‑1 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:1945899/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:1022004/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:685016/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:1619978/1‑68 (MQ=255)            tcttctTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATa                                                   >  1:1182041/1‑68 (MQ=255)               tcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                >  1:2692989/1‑68 (MQ=255)                    cGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGaaa                                           >  1:2651579/1‑68 (MQ=255)                           tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                    >  1:2184766/1‑68 (MQ=255)                           tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                    >  1:147893/1‑68 (MQ=255)                           tATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                    >  1:1436110/1‑68 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:2888655/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:1547015/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:3413647/68‑1 (MQ=255)                               tCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGc                                <  1:2535591/68‑1 (MQ=255)                                 gatggatgCATTTGCTGCATCAATCGGTAAAGAGGCCACCCTCCATAAACCGAAATTTTCTGAAGCAt                              <  1:3211250/68‑1 (MQ=255)                                    ggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGc                           >  1:800397/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATCAACCGAAATTTTCTGAAGCATTGCg                          >  1:1415150/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:197038/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:2621009/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCg                          >  1:4097/1‑68 (MQ=255)                                     gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAACTGTCTGAAGCATTGCg                          >  1:379331/1‑68 (MQ=255)                                       tgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGaa                        <  1:245062/68‑1 (MQ=255)                                           tttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg                    <  1:242951/68‑1 (MQ=255)                                           tttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCgg                    <  1:1935752/68‑1 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   >  1:1038350/1‑68 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   <  1:3322123/68‑1 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   <  1:1382819/68‑1 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   <  1:1657634/68‑1 (MQ=255)                                            ttGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                   >  1:1988691/1‑68 (MQ=255)                                             tgctgcATCAATCGGTAAAGATGCCCCCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGcc                  <  1:978433/68‑1 (MQ=255)                                              gctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCt                 >  1:487017/1‑68 (MQ=255)                                                   aTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttt            <  1:22356/68‑1 (MQ=255)                                                            tAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCg   <  1:2212863/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:1916324/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:764385/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:3351519/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:1006807/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:2666967/68‑1 (MQ=255)                                                             aaaGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGt  <  1:2783150/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 < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑