Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	913,349	G→A	L156L (CTG→TTG)	hcp ←	hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	913,349	0	G	A	100.0%	168.1 / NA	54	L156L (CTG→TTG)	hcp	hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/24);  total (30/24)

GCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCA  >  NC_000913/913235‑913467                                                                                                                   |                                                                                                                       gCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGc                                                                                                                      <  1:946026/117‑1 (MQ=255) gCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGc                                                                                                                      >  2:946026/1‑117 (MQ=255)        ccccAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCg                                                                                         >  1:1172339/1‑139 (MQ=255)           cAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCa                                                                                      >  1:1149324/1‑139 (MQ=255)                  gCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                               >  2:1134658/1‑139 (MQ=255)                    cATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                             <  2:1172339/139‑1 (MQ=255)                      tGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                            <  2:208896/106‑1 (MQ=255)                      tGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACg                                                                                                            >  1:208896/1‑106 (MQ=255)                               ggTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                  <  2:1149324/139‑1 (MQ=255)                                  aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                                                                              >  1:163262/1‑124 (MQ=255)                                  aCTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                                                                              <  2:163262/124‑1 (MQ=255)                                   cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAggg                                                              >  2:1026339/1‑139 (MQ=255)                                   cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAggg                                                              >  1:86572/1‑139 (MQ=255)                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGTCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                           >  2:1143438/1‑139 (MQ=255)                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                           >  2:1022989/1‑139 (MQ=255)                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                           >  2:313851/1‑139 (MQ=255)                                      ggCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTa                                                           >  2:1062156/1‑139 (MQ=255)                                        cATAAATATCGTTGTCGTATTGACCGAGAAAATGCGCGTGTTCCATATAGGACGCCGCACCTTAAAGGCCATAAAAGCACAGAAGACGCAGGCCGAGAATGTTTTCGACAATCGCCGATTTATCTTTGTAAGGGGTaaa                                                         >  1:1052357/1‑139 (MQ=255)                                                tCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                               <  1:602160/109‑1 (MQ=255)                                                tCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                               >  2:602160/1‑109 (MQ=255)                                                 cGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                  >  2:743139/1‑121 (MQ=255)                                                 cGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTCTATCTTTGtt                                                                  <  1:743139/121‑1 (MQ=255)                                                      tCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                 <  2:1014242/117‑1 (MQ=255)                                                      tCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                 >  1:1014242/1‑117 (MQ=255)                                                           ttGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACCGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCATTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGc                                      >  1:1089371/1‑139 (MQ=255)                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTa                                                                          >  2:1044095/1‑98 (MQ=255)                                                                cGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTa                                                                          <  1:1044095/98‑1 (MQ=255)                                                                         gcgcGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                    <  1:715801/127‑1 (MQ=255)                                                                         gcgcGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                    >  2:715801/1‑127 (MQ=255)                                                                          cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGc                                      >  2:287925/1‑124 (MQ=255)                                                                          cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGc                                      <  1:287925/124‑1 (MQ=255)                                                                          cgcgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAtcg                       <  1:1134658/139‑1 (MQ=255)                                                                            cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                 >  1:960582/1‑95 (MQ=255)                                                                            cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                 <  2:960582/95‑1 (MQ=255)                                                                                 tCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCa                <  1:1026339/139‑1 (MQ=255)                                                                                   cATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAcc              <  2:142664/139‑1 (MQ=255)                                                                                     tataGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                  >  2:343041/1‑69 (MQ=255)                                                                                     tataGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                  <  1:343041/69‑1 (MQ=255)                                                                                         ggCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGc        <  1:1062156/139‑1 (MQ=255)                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                   >  1:973692/1‑63 (MQ=255)                                                                                          gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                   <  2:973692/63‑1 (MQ=255)                                                                                           ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAAt                                                                                    <  1:973466/61‑1 (MQ=255)                                                                                           ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAAt                                                                                    >  2:973466/1‑61 (MQ=255)                                                                                           ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc               <  1:833995/130‑1 (MQ=255)                                                                                           ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAc               >  2:833995/1‑130 (MQ=255)                                                                                            cgccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                               >  2:132574/1‑65 (MQ=255)                                                                                            cgccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccg                                                                               <  1:132574/65‑1 (MQ=255)                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt    <  2:1089371/139‑1 (MQ=255)                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt    >  1:738701/1‑139 (MQ=255)                                                                                             gccgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGt    >  2:73582/1‑139 (MQ=255)                                                                                               cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCa  >  2:1071586/1‑139 (MQ=255)                                                                                                                aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                             <  2:628875/47‑1 (MQ=255)                                                                                                                aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                             >  1:628875/1‑47 (MQ=255)                                                                                                                aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGc                               <  1:6078/93‑1 (MQ=255)                                                                                                                aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGc                               >  2:6078/1‑93 (MQ=255)                                                                                                                   |                                                                                                                       GCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCA  >  NC_000913/913235‑913467

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