Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,602,017	(C)7→8	coding (200/2019 nt)	hyfB →	hydrogenase 4, membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,602,010	1	.	C	100.0%	185.9 / NA	54	A65P (GCC→CCC)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (22/32);  total (22/32)

GCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATAT  >  NC_000913/2601875‑2602125                                                                                                                                         |                                                                                                                    gCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑ccc                                                                                                                  >  2:686708/1‑139 (MQ=255)      gCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑cccccccc                                                                                                             >  2:466880/1‑138 (MQ=255)      gCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑cccccccc                                                                                                             >  2:1245891/1‑138 (MQ=255)       cTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑cccccccct                                                                                                            >  1:1181378/1‑137 (MQ=255)                    gTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATAACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCg                                                                                                >  2:931621/1‑139 (MQ=255)                         gACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTa                                                                                           <  1:931621/139‑1 (MQ=255)                              aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       >  2:135179/1‑138 (MQ=255)                              aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       <  1:135179/138‑1 (MQ=255)                              aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCt                                                                                      >  2:540819/1‑139 (MQ=255)                              aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCt                                                                                      >  1:583935/1‑139 (MQ=255)                                 ggCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                                                   <  1:1075834/139‑1 (MQ=255)                                  gCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGa                                                                                  <  2:228196/139‑1 (MQ=255)                                   cTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGaa                                                                                 >  1:1198097/1‑139 (MQ=255)                                        aaGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                                                            <  2:1041160/139‑1 (MQ=255)                                          gCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAg                                                                          <  2:1084653/139‑1 (MQ=255)                                                cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                    <  2:583935/139‑1 (MQ=255)                                                cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                    >  1:818242/1‑139 (MQ=255)                                                cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                                                    >  2:441547/1‑139 (MQ=255)                                                   cATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTATGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                 <  2:373269/139‑1 (MQ=255)                                                            gCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGt                                                        <  1:44090/139‑1 (MQ=255)                                                             cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                       >  1:505140/1‑139 (MQ=255)                                                              tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGc                                                      <  2:1198097/139‑1 (MQ=255)                                                              tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGc                                                      <  2:505140/139‑1 (MQ=255)                                                                  gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGt                                                  <  2:1181378/139‑1 (MQ=255)                                                                     cggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGt                                                  >  1:684805/1‑136 (MQ=255)                                                                     cggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGt                                                  <  2:684805/136‑1 (MQ=255)                                                                     cggcgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAt                                               <  2:192088/139‑1 (MQ=255)                                                                        cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTAtctc                                            >  1:1222894/1‑139 (MQ=255)                                                                          tGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCt                                          >  1:1072311/1‑139 (MQ=255)                                                                            aTTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTg                                        <  1:441547/139‑1 (MQ=255)                                                                             ttGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGc                                       <  2:1072311/139‑1 (MQ=255)                                                                                  atcatcAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                  <  1:725361/139‑1 (MQ=255)                                                                                     atcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTa                                                            >  1:326862/1‑110 (MQ=255)                                                                                     atcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTa                                                            <  2:326862/110‑1 (MQ=255)                                                                                      tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtgg                              >  1:1088206/1‑139 (MQ=255)                                                                                       caGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTCAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                             <  2:929957/139‑1 (MQ=255)                                                                                              aTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCctc                                                                      >  1:520135/1‑91 (MQ=255)                                                                                              aTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCctc                                                                      <  2:520135/91‑1 (MQ=255)                                                                                                        tACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCCATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt            <  2:144537/139‑1 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCctct                                                                     <  2:757353/73‑1 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCctct                                                                     >  1:757353/1‑73 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc                       >  2:659017/1‑119 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc                       <  1:659017/119‑1 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa             >  1:491363/1‑129 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCa             <  2:491363/129‑1 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata   <  1:540819/139‑1 (MQ=255)                                                                                                                 gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCCGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata   <  2:1088206/139‑1 (MQ=255)                                                                                                                  tgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       <  1:565077/54‑1 (MQ=255)                                                                                                                  tgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       >  2:565077/1‑54 (MQ=255)                                                                                                                  tgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTatat  <  2:818242/139‑1 (MQ=255)                                                                                                                    tAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       >  1:1279967/1‑52 (MQ=255)                                                                                                                    tAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAAccc                                                                                       <  2:1279967/52‑1 (MQ=255)                                                                                                                           aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                             >  1:765256/1‑103 (MQ=255)                                                                                                                           aGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                             <  2:765256/103‑1 (MQ=255)                                                                                                                                cccGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGc                                                                                    >  2:16198/1‑43 (MQ=255)                                                                                                                                cccGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGc                                                                                    <  1:16198/43‑1 (MQ=255)                                                                                                                                cccGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                       >  1:413850/1‑72 (MQ=255)                                                                                                                                cccGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                       <  2:413850/72‑1 (MQ=255)                                                                                                                                         |                                                                                                                    GCTTCGCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑CCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATAT  >  NC_000913/2601875‑2602125

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