Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,852,356	T→A	100%	L74Q (CTG→CAG)	hypD →	hydrogenase maturation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,852,356	0	T	A	100.0%	299.2 / ‑6.7	89	L74Q (CTG→CAG)	hypD	hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (54/34);  minor base C (1/0);  total (55/34)

GCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852260‑2852452                                                                                                 |                                                                                                 gCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGcc                                                                                               <  1:694613/100‑1 (MQ=255) gCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGcc                                                                                               <  1:1997505/100‑1 (MQ=255)  cGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCg                                                                                              >  1:2209031/1‑100 (MQ=255)  cGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCg                                                                                              <  1:2397288/100‑1 (MQ=255)    gTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGAt                                                                                            <  1:489854/100‑1 (MQ=255)    gTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGAt                                                                                            <  2:1488358/100‑1 (MQ=255)    gTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACATTGAGTTTATACACGGTCCAGGGTGCCCGGTGTGCGTACAGcc                                                                                               >  2:898451/1‑97 (MQ=255)     tCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATg                                                                                           >  2:201392/1‑100 (MQ=255)     tCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATg                                                                                           >  2:2340647/1‑100 (MQ=255)      cATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATgg                                                                                          >  2:1369627/1‑100 (MQ=255)       aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  1:2197550/100‑1 (MQ=255)       aTACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         <  1:2260349/100‑1 (MQ=255)        tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATggg                                                                                         >  1:948532/1‑99 (MQ=255)        tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        >  2:2076457/1‑100 (MQ=255)        tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGt                                                                                        >  2:104159/1‑100 (MQ=255)        tACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCCGTGTGCGTACAGCCGATGGGt                                                                                        <  2:2108405/100‑1 (MQ=255)        tACCCACGCTATCGTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACCGTCCGGGGTGCCCAGTGTGCGTACAGCCGATGGGt                                                                                        >  2:108655/1‑100 (MQ=255)         aCCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTa                                                                                       >  1:1084954/1‑100 (MQ=255)            cACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGaa                                                                                    <  2:1700782/100‑1 (MQ=255)               gCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGTTTTGCGTACAGCCGATGGGTAGAATCg                                                                                 <  2:165944/100‑1 (MQ=255)                 tATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGac                                                                               <  1:1990052/100‑1 (MQ=255)                    cTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACAcc                                                                            >  1:369671/1‑100 (MQ=255)                    cTTTAAATTCGGCCTCGACCAGTTACTGCCGCAAAACGATGAGATTATCCACGGTCCGGGGTACCCCATGTGCGTACCGCCGATGGGGAGCATCGccac                                                                             >  2:1517689/1‑99 (MQ=255)                     tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGCATCGACACCt                                                                           >  2:968409/1‑100 (MQ=255)                     tttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCt                                                                           >  1:1172421/1‑100 (MQ=255)                      ttAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTg                                                                          >  1:2005774/1‑100 (MQ=255)                         aaTTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGt                                                                       >  1:880320/1‑100 (MQ=255)                          aTTCGGCCTCGACCAGTTACTGCCGAAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCTGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTg                                                                      <  1:1444612/100‑1 (MQ=255)                              ggCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGgaga                                                                  >  2:1319594/1‑100 (MQ=255)                               gCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAt                                                                 >  1:728258/1‑100 (MQ=255)                                ccTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGAt                                                                 >  2:2128133/1‑99 (MQ=255)                                     aCCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccag                                                           >  1:786800/1‑100 (MQ=255)                                        aGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcca                                                        >  1:236152/1‑100 (MQ=255)                                         gTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGTTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTgccagcca                                                        >  2:1816491/1‑99 (MQ=255)                                         gTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCAt                                                       >  1:1633273/1‑100 (MQ=255)                                             cTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCgg                                                   >  2:1773148/1‑100 (MQ=255)                                              tGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGa                                                  >  2:1711418/1‑100 (MQ=255)                                               gCCGGCAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCCGCCATCCGGaa                                                 >  2:1922292/1‑100 (MQ=255)                                               gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGaa                                                 >  2:1139351/1‑100 (MQ=255)                                               gCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGaa                                                 >  1:1968888/1‑100 (MQ=255)                                                ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                >  1:118964/1‑100 (MQ=255)                                                ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                <  2:1293786/100‑1 (MQ=255)                                                ccGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                <  2:703548/100‑1 (MQ=255)                                                 cGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                >  1:1152734/1‑99 (MQ=255)                                                 cGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAg                                                >  1:2197899/1‑99 (MQ=255)                                                  ggAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTc                                              >  1:978618/1‑100 (MQ=255)                                                   gAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCa                                             >  1:1597865/1‑100 (MQ=255)                                                   gAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCAGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCa                                             >  2:825304/1‑100 (MQ=255)                                                     aaaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtc                                           >  2:1878819/1‑100 (MQ=255)                                                      aaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtct                                          <  1:2188895/100‑1 (MQ=255)                                                      aaCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtct                                          <  1:1369641/100‑1 (MQ=255)                                                       aCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtctt                                         >  1:1814794/1‑100 (MQ=255)                                                       aCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtctt                                         >  2:1778354/1‑100 (MQ=255)                                                       aCGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCAAGCCATAAGGAAGTAAtatt                                         >  1:1807144/1‑100 (MQ=255)                                                        cGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCAtcttc                                        <  2:595348/100‑1 (MQ=255)                                                          ttGAGTTTATCCACGGTACGAGGTACCCGGGGTACGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCACTGCGGAAGTAAtctt                                         >  2:1811213/1‑97 (MQ=255)                                                             aGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAc                                   >  1:1901729/1‑100 (MQ=255)                                                              gTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTAcc                                  <  1:378407/100‑1 (MQ=255)                                                                ttATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCtt                                >  2:959601/1‑100 (MQ=255)                                                                ttATCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCtt                                >  2:509757/1‑100 (MQ=255)                                                                   tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCttttg                             >  2:1857876/1‑98 (MQ=255)                                                                   tCCACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTgg                             >  1:1714236/1‑100 (MQ=255)                                                                    ccAGGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            <  2:1076998/100‑1 (MQ=255)                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:427812/1‑100 (MQ=255)                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGc                            >  2:774424/1‑100 (MQ=255)                                                                    ccACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTATgg                             >  1:1827239/1‑99 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:2366326/1‑100 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  1:376230/100‑1 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  1:559204/100‑1 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           >  2:827975/1‑100 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  2:1691790/100‑1 (MQ=255)                                                                     cACGGTCCGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCg                           <  2:1601220/100‑1 (MQ=255)                                                                           ccGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCa                     <  1:631529/100‑1 (MQ=255)                                                                           ccGGGGTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCa                     <  1:1094729/100‑1 (MQ=255)                                                                             ggggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATg                   <  2:1427814/100‑1 (MQ=255)                                                                             ggggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATg                   <  1:755560/100‑1 (MQ=255)                                                                              gggTGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATgc                  <  1:1954953/100‑1 (MQ=255)                                                                                 tGCCCGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGcgcg               >  1:2203337/1‑100 (MQ=255)                                                                                   cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             <  1:729740/100‑1 (MQ=255)                                                                                   cccGGTGTGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTg             <  1:1778375/100‑1 (MQ=255)                                                                                      ggTGTGCGTACAGCCGATGGGTAGAATCGATACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCg          <  2:475211/100‑1 (MQ=255)                                                                                       gtgtGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgg         <  2:645446/100‑1 (MQ=255)                                                                                       gtgtGCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCg          >  2:2256981/1‑99 (MQ=255)                                                                                        tgtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCggg        >  2:1204652/1‑100 (MQ=255)                                                                                         gtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGCAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgggg       >  1:461942/1‑100 (MQ=255)                                                                                         gtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgggg       <  1:1319608/100‑1 (MQ=255)                                                                                         gtgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCgggg       <  2:2017436/100‑1 (MQ=255)                                                                                          tgCGTACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGa      >  1:591054/1‑100 (MQ=255)                                                                                              tACAGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACa  >  1:2306598/1‑100 (MQ=255)                                                                                                 |                                                                                                 GCGGTCATACCCACGCTATCTTTAAATTCGGCCTCGACCAGTTACTGCCGGAAAACGTTGAGTTTATCCACGGTCCGGGGTGCCCGGTGTGCGTACTGCCGATGGGTAGAATCGACACCTGCGTGGAGATTGCCAGCCATCCGGAAGTCATCTTCTGTACCTTTGGCGACGCCATGCGCGTGCCGGGGAAACA  >  NC_000913/2852260‑2852452

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