Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,847,704	A→T	L513Q (CTG→CAG)	hycC ←	formate hydrogenlyase subunit HycC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,847,704	0	A	T	100.0%	132.6 / NA	39	L513Q (CTG→CAG)	hycC	formate hydrogenlyase subunit HycC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/23);  total (16/23)

CGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGC  >  NC_000913/2847577‑2847836                                                                                                                                |                                                                                                                                     cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGggcagg                                                                                                                          <  1:621310/140‑1 (MQ=255)   gCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGggcaggca                                                                                                                        <  1:548518/140‑1 (MQ=255)        gCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATcag                                                                                                                   <  2:368881/140‑1 (MQ=255)                 tGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGt                                                                                                          >  1:242725/1‑140 (MQ=255)                  gAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGggcaggca                                                                                                                        <  1:210646/125‑1 (MQ=255)                  gAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGggcaggca                                                                                                                        >  2:210646/1‑125 (MQ=255)                     cGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGatc                                                                                                      <  2:257163/140‑1 (MQ=255)                       gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCTCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGatcat                                                                                                    <  2:709092/140‑1 (MQ=255)                       gTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGatcat                                                                                                    <  1:475294/140‑1 (MQ=255)                             aCCATTGATTTTTCGTGGTCGTAACCGCACCCCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCCGCAGCAACGTGATCATCGGTTg                                                                                              <  2:231085/140‑1 (MQ=255)                                 ttGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                          >  2:254915/1‑140 (MQ=255)                                  tGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaaa                                                                                         >  1:108981/1‑140 (MQ=255)                                    aTTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACg                                                                                       <  2:88080/140‑1 (MQ=255)                                    aTTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACg                                                                                       <  2:242725/140‑1 (MQ=255)                                       tttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtg                                                                                    <  1:817490/140‑1 (MQ=255)                                              gTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                             >  1:171761/1‑140 (MQ=255)                                                cGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCa                                                                           <  2:841011/140‑1 (MQ=255)                                                        cacaCCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAg                                                                   >  1:712761/1‑140 (MQ=255)                                                         acacCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGc                                                                  <  1:575194/140‑1 (MQ=255)                                                            cccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGc                                                               <  2:19427/140‑1 (MQ=255)                                                             ccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATGATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCa                                                              >  1:453792/1‑140 (MQ=255)                                                             ccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGc                                                               <  1:101462/139‑1 (MQ=255)                                                             ccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGc                                                               >  2:101462/1‑139 (MQ=255)                                                                     cACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAca                                                      >  1:858928/1‑140 (MQ=255)                                                                     cACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAca                                                      >  2:439116/1‑140 (MQ=255)                                                                         gcgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcag                                                  >  2:150168/1‑140 (MQ=255)                                                                          cgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcagc                                                 >  1:167746/1‑140 (MQ=255)                                                                                   gATGGCAAACTATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGCGATCATCGGTTGAGAAGCGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                        <  2:455387/140‑1 (MQ=255)                                                                                   gATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                        <  1:150168/140‑1 (MQ=255)                                                                                   gATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                        <  2:724278/140‑1 (MQ=255)                                                                                                 gCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                        <  2:149703/126‑1 (MQ=255)                                                                                                 gCCTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCAt                                        >  1:149703/1‑126 (MQ=255)                                                                                                  ccTTTGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACgg                         <  1:254915/140‑1 (MQ=255)                                                                                                      tGCAAATCGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgca                     <  1:247269/140‑1 (MQ=255)                                                                                                             cGCCATGATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac              <  2:453792/140‑1 (MQ=255)                                                                                                             cGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAacac              >  2:263336/1‑140 (MQ=255)                                                                                                             cGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAaca               <  2:440186/139‑1 (MQ=255)                                                                                                             cGCCATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAaca               >  1:440186/1‑139 (MQ=255)                                                                                                                cATAATGATGAATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGATAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGc           >  2:489099/1‑140 (MQ=255)                                                                                                                         gaATGGCTGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACgc  <  2:644548/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                     CGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGC  >  NC_000913/2847577‑2847836

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