Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,511,468	C→T	*413* (TAG→TAA)	mntH ←	manganese/divalent cation transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,511,468	0	C	T	100.0%	135.4 / NA	43	*413* (TAG→TAA)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/24);  total (19/24)

TGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCAC  >  NC_000913/2511351‑2511592                                                                                                                      |                                                                                                                             tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCg                                                                                                              >  1:1007680/1‑134 (MQ=255) tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCg                                                                                                              <  2:1007680/134‑1 (MQ=255)      gATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                    <  2:662074/139‑1 (MQ=255)             cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAg                                                                                             >  2:332102/1‑139 (MQ=255)                   tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                           <  2:425075/103‑1 (MQ=255)                   tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                           >  1:425075/1‑103 (MQ=255)                              aGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                            <  1:337780/139‑1 (MQ=255)                               ggCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                           >  1:638275/1‑139 (MQ=255)                                gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCa                                                                          >  2:116146/1‑139 (MQ=255)                                 ccTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGGCAGATGCGACGCTTATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAg                                                                         <  2:748186/139‑1 (MQ=255)                                            gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                           <  2:360901/78‑1 (MQ=255)                                            gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                           >  1:360901/1‑78 (MQ=255)                                            gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                                               <  1:14977/106‑1 (MQ=255)                                            gCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                                               >  2:14977/1‑106 (MQ=255)                                               ctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAccc                                                           <  1:277501/139‑1 (MQ=255)                                                 ttAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                         <  1:116146/139‑1 (MQ=255)                                                         tAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                            <  1:20529/112‑1 (MQ=255)                                                         tAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                            >  2:20529/1‑112 (MQ=255)                                                          aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                                                >  2:202583/1‑139 (MQ=255)                                                            ggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGttt                                              >  2:927006/1‑139 (MQ=255)                                                                      gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                    >  1:702576/1‑139 (MQ=255)                                                                      gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                    <  2:608454/139‑1 (MQ=255)                                                                      gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                    <  2:638275/139‑1 (MQ=255)                                                                      gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                    <  2:702576/139‑1 (MQ=255)                                                                           tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGcc                                                                                 <  1:318874/89‑1 (MQ=255)                                                                           tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGcc                                                                                 >  2:318874/1‑89 (MQ=255)                                                                           tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                  >  1:415898/1‑136 (MQ=255)                                                                           tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                  <  2:415898/136‑1 (MQ=255)                                                                                  cATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTc                              <  1:330041/133‑1 (MQ=255)                                                                                  cATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTc                              >  2:330041/1‑133 (MQ=255)                                                                                         cAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCa                 <  1:166179/139‑1 (MQ=255)                                                                                                     aTGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                         <  2:14182/87‑1 (MQ=255)                                                                                                     aTGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                         >  1:14182/1‑87 (MQ=255)                                                                                                      tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTc    <  2:346168/139‑1 (MQ=255)                                                                                                        cGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCAc  <  1:202583/139‑1 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAg                                                                                             <  1:242547/47‑1 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAg                                                                                             >  2:242547/1‑47 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcac                                                                      <  1:986842/70‑1 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcac                                                                      >  2:986842/1‑70 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                  <  2:831811/106‑1 (MQ=255)                                                                                                         gACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                  >  1:831811/1‑106 (MQ=255)                                                                                                           cGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                  >  2:357539/1‑40 (MQ=255)                                                                                                           cGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                  <  1:357539/40‑1 (MQ=255)                                                                                                                      |                                                                                                                             TGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCAC  >  NC_000913/2511351‑2511592

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