Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	913,520	Δ1 bp	coding (295/1653 nt)	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,518	0	T	.	100.0%	206.5 / NA	43	coding (297/1653 nt)	hcp	hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/25);  total (18/25)

GAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGCTTTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCACGCGCT  >  NC_000913/913368‑913662                                                                                                                                                       |                                                                                                                                                 gagaATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGCt                                                                                                                                                  >  2:142637/1‑151 (MQ=255)             gCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCCCCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGCttcag                                                                                                                                              <  2:540285/143‑4 (MQ=255)             gCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTcgc                                                                                                                                     >  2:143692/1‑151 (MQ=255)                   cgccgcTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAggg                                                                                                                               >  2:605145/1‑151 (MQ=255)                               tttGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCAc                                                                                                                   >  1:44022/1‑151 (MQ=255)                                 tGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACg                                                                                                                  >  1:570362/1‑150 (MQ=255)                                      ggggTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTa                                                                                                            <  2:503706/151‑1 (MQ=255)                                                 tgctgcCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATAc                                                                                                 >  2:127293/1‑151 (MQ=255)                                                         gACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGgaga                                                                                          <  2:603703/150‑1 (MQ=255)                                                         gACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGgaga                                                                                          <  1:229169/150‑1 (MQ=255)                                                                cAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCg                                                                                      >  2:657090/1‑147 (MQ=255)                                                                cAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCg                                                                                      <  1:657079/147‑1 (MQ=255)                                                                      gCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGt                                                                            <  2:75099/151‑1 (MQ=255)                                                                       ccGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGtt                                                                           <  2:403043/151‑1 (MQ=255)                                                                        cGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTg                                                                          <  1:569061/151‑1 (MQ=255)                                                                        cGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTg                                                                          <  1:576290/151‑1 (MQ=255)                                                                         gagaTCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTgg                                                                         >  1:184515/1‑151 (MQ=255)                                                                         gagaTCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTgg                                                                         <  1:157225/151‑1 (MQ=255)                                                                          agagCGTCGCCTACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCTTTTGCAGCTACAGCCAGGCATTGCGC‑TTCAGCGCATCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTg                                                                          <  1:611672/149‑1 (MQ=255)                                                                              cgtcgCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAg                                                                     <  1:49724/150‑1 (MQ=255)                                                                               gtcgCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAggg                                                                   >  2:343560/1‑151 (MQ=255)                                                                                     cACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTga                                                               <  1:662504/149‑1 (MQ=255)                                                                                     cACCAGTGGCAGGTCAGCCATCGGGTTATCGTCGCGCGCGTTGGCAGCTACAGCCTGGCATTGCGC‑TTCATCGCCTCGCGCAGAGCACTCTCTCCACGAGCGTAGCCGACAAGACGCGGAGAATCGAAGTTAACGTTGGTCAGTGTTgaga                                                             <  1:332722/151‑1 (MQ=255)                                                                                        cAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGaaa                                                           >  1:237918/1‑150 (MQ=255)                                                                                             gCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCa                                                     <  2:531907/151‑1 (MQ=255)                                                                                                ggTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGt                                                  <  1:127290/151‑1 (MQ=255)                                                                                                 gTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTg                                                 >  2:104034/1‑151 (MQ=255)                                                                                                        aTCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGaa                                          <  1:605139/151‑1 (MQ=255)                                                                                                           gggTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCt                                       >  1:91598/1‑151 (MQ=255)                                                                                                            ggTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTg                                      <  1:293089/151‑1 (MQ=255)                                                                                                              ttATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTc                                    >  2:171808/1‑151 (MQ=255)                                                                                                               tATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTc                                    <  2:300085/150‑1 (MQ=255)                                                                                                               tATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTc                                    >  1:131076/1‑150 (MQ=255)                                                                                                                   gACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACAt                               >  2:267640/1‑151 (MQ=255)                                                                                                                       cgcgcgTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTg                           <  1:104031/151‑1 (MQ=255)                                                                                                                        gcgcgTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTgg                          >  2:517073/1‑151 (MQ=255)                                                                                                                           cgTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTtg                       <  1:648369/151‑1 (MQ=255)                                                                                                                           cgTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTtg                       <  2:97307/151‑1 (MQ=255)                                                                                                                            gTTGGCATCTACAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTtg                       >  2:23625/1‑150 (MQ=255)                                                                                                                                      aCAGCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGcc               <  1:143688/148‑1 (MQ=255)                                                                                                                                         gCCAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTc         <  1:62538/151‑1 (MQ=255)                                                                                                                                          ccAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTc         >  1:207340/1‑150 (MQ=255)                                                                                                                                           cAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCAc       <  1:707/151‑1 (MQ=255)                                                                                                                                           cAGGCATTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCAc       >  1:531109/1‑151 (MQ=255)                                                                                                                                                aTTGCGC‑TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCACGCGCt  <  1:171806/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTATCGACGCGCGCGTTGGCATCTACAGCCAGGCATTGCGCTTTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCACGCGCT  >  NC_000913/913368‑913662

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