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

GTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATC  >  NC_000913/2511332‑2511603                                                                                                                                         |                                                                                                                                        gTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTa                                                                                                                                       >  2:805242/1‑139 (MQ=255)                   aTGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                     <  2:296577/139‑1 (MQ=255)                    tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTcccc                                                                                                                    >  1:1274234/1‑139 (MQ=255)                             cGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCaacaac                                                                                                           <  2:390525/139‑1 (MQ=255)                                      tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            >  1:246937/1‑113 (MQ=255)                                      tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            <  2:246937/113‑1 (MQ=255)                                          tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                               >  2:495237/1‑122 (MQ=255)                                          tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                               <  1:495237/122‑1 (MQ=255)                                          tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                                                              >  1:992887/1‑139 (MQ=255)                                             tATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCa                                                                                           <  2:384885/139‑1 (MQ=255)                                                      tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGca                                                                                  >  2:1321330/1‑139 (MQ=255)                                                       aCAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcac                                                                                 <  2:1274234/139‑1 (MQ=255)                                                           aCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcaccact                                                                             <  2:840343/139‑2 (MQ=255)                                                                    ttAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                                    <  1:1321330/139‑1 (MQ=255)                                                                             aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                                                           >  1:593742/1‑139 (MQ=255)                                                                                    aTAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgc                                                    >  1:433058/1‑139 (MQ=255)                                                                                       aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                                                 <  1:805242/139‑1 (MQ=255)                                                                                         gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                               >  1:707706/1‑139 (MQ=255)                                                                                         gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                               <  2:707706/139‑1 (MQ=255)                                                                                         gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                               <  2:593742/139‑1 (MQ=255)                                                                                              tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            <  2:578965/57‑1 (MQ=255)                                                                                              tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            >  1:578965/1‑57 (MQ=255)                                                                                                cgccgcATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCa                                        <  2:992887/139‑1 (MQ=255)                                                                                                               tAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCaacaac                                                                                                           <  2:739575/57‑1 (MQ=255)                                                                                                               tAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCaacaac                                                                                                           >  1:739575/1‑57 (MQ=255)                                                                                                                  tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                      <  2:352001/75‑1 (MQ=255)                                                                                                                  tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                      >  1:352001/1‑75 (MQ=255)                                                                                                                  tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                           >  2:1301647/1‑86 (MQ=255)                                                                                                                  tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                           <  1:1301647/86‑1 (MQ=255)                                                                                                                     cAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGcc                                                                                            <  2:148224/66‑1 (MQ=255)                                                                                                                     cAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGcc                                                                                            >  1:148224/1‑66 (MQ=255)                                                                                                                     cAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGc                   >  2:598122/1‑139 (MQ=255)                                                                                                                             aCGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                                                 <  1:1008881/101‑1 (MQ=255)                                                                                                                             aCGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGAACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                                                 >  2:1008881/1‑101 (MQ=255)                                                                                                                                  cATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGaa      >  1:681288/1‑139 (MQ=255)                                                                                                                                   aTTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGaaa     >  1:775047/1‑139 (MQ=255)                                                                                                                                     tCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAAt   <  1:598122/139‑1 (MQ=255)                                                                                                                                      cAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATc  <  2:681288/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        GTCCGCCAATGGAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATC  >  NC_000913/2511332‑2511603

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