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

CCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAA  >  NC_000913/2511353‑2511600                                                                                                                    |                                                                                                                                     ccTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccac                                                                                                               >  1:289837/1‑139 (MQ=255)  cTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCccacc                                                                                                              <  2:436094/139‑1 (MQ=255)    gATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                            <  2:296564/139‑1 (MQ=255)         gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                            <  1:232199/134‑1 (MQ=255)         gACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                            >  2:232199/1‑134 (MQ=255)                     tCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAgcg                                                                                           >  2:728005/1‑139 (MQ=255)                              gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                         <  2:311058/100‑1 (MQ=255)                              gCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAgcg                                                                                                                         >  1:311058/1‑100 (MQ=255)                                 tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                                    <  2:89761/134‑1 (MQ=255)                                 tACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGAc                                                                                    >  1:89761/1‑134 (MQ=255)                                   cAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATGGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACg                                                                                      <  1:753719/130‑1 (MQ=255)                                   cAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATGGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACg                                                                                      >  2:753719/1‑130 (MQ=255)                                        cGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAAt                                                                        >  2:372181/1‑139 (MQ=255)                                         ggCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATc                                                                       >  2:138133/1‑139 (MQ=255)                                          gCTCTTTAACGCGTAGGGCGGGTAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCATCAACCAGATATTCAGCGCCACGGCCAGCACCACAATCa                                                                      <  2:657352/139‑1 (MQ=255)                                           ctctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAc                                                                     >  2:622598/1‑139 (MQ=255)                                            tctTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCAcc                                                                    <  1:138133/139‑1 (MQ=255)                                                    gcgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgt                                                            <  2:289837/139‑1 (MQ=255)                                                     cgTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtc                                                           >  1:46735/1‑139 (MQ=255)                                                        aGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgtctgt                                                        >  1:308765/1‑139 (MQ=255)                                                          ggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgt                                                            <  2:675477/133‑1 (MQ=255)                                                          ggCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCctgt                                                            >  1:675477/1‑133 (MQ=255)                                                             ggATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                                                          >  2:771116/1‑116 (MQ=255)                                                             ggATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACcaca                                                                          <  1:771116/116‑1 (MQ=255)                                                                  aGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGtt                                              <  2:152449/139‑1 (MQ=255)                                                                    gCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTg                                            <  1:338382/139‑1 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                  <  1:39976/44‑1 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTccc                                                                                                                  >  2:39976/1‑44 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                            <  1:681215/50‑1 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACcaa                                                                                                            >  2:681215/1‑50 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                                 >  1:139456/1‑93 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAg                                                                 <  2:139456/93‑1 (MQ=255)                                                                                             tGCCAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACt                   >  2:649186/1‑139 (MQ=255)                                                                                               ccAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcac                                                                              <  1:535232/78‑1 (MQ=255)                                                                                               ccAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGcac                                                                              >  2:535232/1‑78 (MQ=255)                                                                                               ccAGATGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTg                 >  2:325222/1‑139 (MQ=255)                                                                                                    tGCGACGCTCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTc            >  2:290811/1‑139 (MQ=255)                                                                                                            tCATTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGa    >  2:64004/1‑139 (MQ=255)                                                                                                              aTTCAATTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGaaa  >  1:135346/1‑139 (MQ=255)                                                                                                                    |                                                                                                                                     CCTGATGCGACGCTAGTGCGTCTTATCAGGCCTACAAACCGGCTCTTTAACGCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGCCAATAGTGCCAGATGCGACGCTCATTCAACTACAATCCCAGCGCCGTCCCCACCAACAACCAGATATTCAGCGCCACGACCAGCACCACAATCACCCAGCCTGTCTGTTTTACGCGTTTGCTGTTCACCAGATCGCCCATCAACTTGCTGTCACTGGTGAAA  >  NC_000913/2511353‑2511600

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: ‑