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

GAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATC  >  NC_000913/2511343‑2511603                                                                                                                              |                                                                                                                                        gAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            <  2:594593/139‑1 (MQ=255)  aGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcgc                                                                                                                           >  1:677577/1‑139 (MQ=255)         tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTcccc                                                                                                                    >  2:1191760/1‑139 (MQ=255)         tGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTcccc                                                                                                                    >  2:1580523/1‑139 (MQ=255)            cTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 <  1:761684/139‑1 (MQ=255)                 gCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGcc                                                                                                                          >  1:1086872/1‑125 (MQ=255)                 gCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGcc                                                                                                                          <  2:1086872/125‑1 (MQ=255)                   gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            <  1:1078948/121‑1 (MQ=255)                   gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                            >  2:1078948/1‑121 (MQ=255)                   gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACCCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAAcc                                                                                                          >  2:125853/1‑139 (MQ=255)                    aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCg                                                                                                                         >  1:1127459/1‑123 (MQ=255)                    aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCg                                                                                                                         <  2:1127459/123‑1 (MQ=255)                    aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                             <  1:1670645/135‑1 (MQ=255)                    aCGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaaca                                                                                                             >  2:1670645/1‑135 (MQ=255)                     cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATcc                                                                                                                                 >  2:245575/1‑114 (MQ=255)                     cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATcc                                                                                                                                 <  1:245575/114‑1 (MQ=255)                     cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCcca                                                                                                                   <  1:755799/128‑1 (MQ=255)                     cGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCcca                                                                                                                   >  2:755799/1‑128 (MQ=255)                       cTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGat                                                                                                      <  2:1620247/139‑1 (MQ=255)                           tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                                      <  1:20439/103‑1 (MQ=255)                           tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTAc                                                                                                                                      >  2:20439/1‑103 (MQ=255)                           tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 >  2:1457644/1‑124 (MQ=255)                           tGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 <  1:1457644/124‑1 (MQ=255)                               tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                                                              >  2:1247567/1‑139 (MQ=255)                               tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                                                              >  1:1316783/1‑139 (MQ=255)                                     cAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGa                                                                                        >  2:9929/1‑139 (MQ=255)                                      aGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                                       <  1:1247567/139‑1 (MQ=255)                                        gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAATGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                      >  2:1416328/1‑138 (MQ=255)                                        gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAATGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAcc                                                                                      <  1:1416328/138‑1 (MQ=255)                                                aCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                                                             <  2:1316783/139‑1 (MQ=255)                                                aCCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                                                             <  1:125853/139‑1 (MQ=255)                                                  cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                               <  2:155239/103‑1 (MQ=255)                                                  cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                               >  1:155239/1‑103 (MQ=255)                                                  cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                           >  1:1352894/1‑139 (MQ=255)                                                                  aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                                                           >  2:911921/1‑139 (MQ=255)                                                                        gATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCaacaac                                                                                                           >  2:661849/1‑85 (MQ=255)                                                                        gATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCaacaac                                                                                                           <  1:661849/85‑1 (MQ=255)                                                                         aTAAGGCGTTTACGCCGCATCCGCCAATAGCGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTAcgc                                                    <  2:314051/139‑1 (MQ=255)                                                                              gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 >  1:231818/1‑73 (MQ=255)                                                                              gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 <  2:231818/73‑1 (MQ=255)                                                                              gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                                                       <  2:845382/115‑1 (MQ=255)                                                                              gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                                                       >  1:845382/1‑115 (MQ=255)                                                                              gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                               <  1:583299/139‑1 (MQ=255)                                                                               cGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGc                                              <  1:1191760/139‑1 (MQ=255)                                                                                gTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCt                                             <  1:1580523/139‑1 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 >  2:524577/1‑68 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 <  1:524577/68‑1 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccac                                                                                                                  >  1:112046/1‑67 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccac                                                                                                                  <  2:112046/67‑1 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                                                       >  2:207567/1‑110 (MQ=255)                                                                                   tACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                                                       <  1:207567/110‑1 (MQ=255)                                                                                                 cAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCa                            <  1:9929/139‑1 (MQ=255)                                                                                                       tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacca                                                                                                                <  2:1470841/49‑1 (MQ=255)                                                                                                       tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacca                                                                                                                >  1:1470841/1‑49 (MQ=255)                                                                                                             aTGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGt                >  1:1525630/1‑139 (MQ=255)                                                                                                              tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTc               >  1:1797535/1‑139 (MQ=255)                                                                                                              tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTc               >  2:1626627/1‑139 (MQ=255)                                                                                                                 gACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 <  2:885156/38‑1 (MQ=255)                                                                                                                 gACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                                 >  1:885156/1‑38 (MQ=255)                                                                                                                  aCGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTg           >  1:36183/1‑139 (MQ=255)                                                                                                                           cAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATc  <  1:911921/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                        GAGCACAATGCCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAAATC  >  NC_000913/2511343‑2511603

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