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%	131.9 / NA	41	*413* (TAG→TAA)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/22);  total (19/22)

TCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAAT  >  NC_000913/2511333‑2511602                                                                                                                                        |                                                                                                                                       tCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCACCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                                     >  2:737281/1‑139 (MQ=255)                  aTGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                    <  1:109948/139‑1 (MQ=255)                   tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTGATTCAATTACAATCCCAGCGCCGTcccc                                                                                                                   >  1:665814/1‑139 (MQ=255)                     ccTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccac                                                                                                                 <  2:585234/139‑1 (MQ=255)                      cTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                <  2:691174/139‑1 (MQ=255)                          tGCGACGCTAGTGGGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                            <  1:509798/139‑1 (MQ=255)                          tGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                            <  1:686006/139‑1 (MQ=255)                                                aGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGAGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                                      <  2:665814/139‑1 (MQ=255)                                                 ggCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                            >  1:17657/1‑116 (MQ=255)                                                 ggCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                            <  2:17657/116‑1 (MQ=255)                                                      aCAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                     >  1:358850/1‑134 (MQ=255)                                                      aCAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                     <  2:358850/134‑1 (MQ=255)                                                            cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                          >  1:691658/1‑139 (MQ=255)                                                            cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                          >  2:307026/1‑139 (MQ=255)                                                                       cgcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctg                                                               <  2:691658/139‑1 (MQ=255)                                                                            aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                                                          >  2:226560/1‑139 (MQ=255)                                                                              ggCGGATAAGGTGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                                                         <  2:890803/90‑1 (MQ=255)                                                                              ggCGGATAAGGTGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                                                         >  1:890803/1‑90 (MQ=255)                                                                                      aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGata                                                                                                    >  1:3695/1‑87 (MQ=255)                                                                                      aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGata                                                                                                    <  2:3695/87‑1 (MQ=255)                                                                                      aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                                                <  1:154893/139‑1 (MQ=255)                                                                                        gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              <  1:307026/139‑1 (MQ=255)                                                                                        gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              <  1:446675/139‑1 (MQ=255)                                                                                        gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              >  2:239946/1‑139 (MQ=255)                                                                                        gCGTTTACGCCGCATCCACCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              <  1:737281/139‑1 (MQ=255)                                                                                           tttACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCcca                                                                                                                  <  2:492148/68‑1 (MQ=255)                                                                                           tttACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCcca                                                                                                                  >  1:492148/1‑68 (MQ=255)                                                                                                   gcATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                    >  1:1081768/1‑58 (MQ=255)                                                                                                   gcATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                    <  2:1081768/58‑1 (MQ=255)                                                                                                     aTCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGAt                                 <  2:153553/139‑1 (MQ=255)                                                                                                                 tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCa                                                                                                <  2:334851/64‑1 (MQ=255)                                                                                                                 tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCa                                                                                                >  1:334851/1‑64 (MQ=255)                                                                                                                     aGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCt                 <  2:883168/139‑1 (MQ=255)                                                                                                                     aGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCt                 >  1:883168/1‑139 (MQ=255)                                                                                                                        tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              <  2:159737/107‑1 (MQ=255)                                                                                                                        tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                              >  1:159737/1‑107 (MQ=255)                                                                                                                          cGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCAc            >  1:588217/1‑139 (MQ=255)                                                                                                                          cGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCAc            >  2:114286/1‑139 (MQ=255)                                                                                                                             cGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                            >  2:1041983/1‑104 (MQ=255)                                                                                                                             cGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                            <  1:1041983/104‑1 (MQ=255)                                                                                                                                    tCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAAt  >  1:588891/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       TCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAAT  >  NC_000913/2511333‑2511602

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