Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,601,372	C→T	P58S (CCT→TCT)	hyfA →	hydrogenase 4 component A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,601,372	0	C	T	100.0%	124.7 / NA	37	P58S (CCT→TCT)	hyfA	hydrogenase 4 component A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/13);  total (24/13)

GTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGT  >  NC_000913/2601133‑2601614                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   gTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGc                                                                                                                                                                                                                                           >  1:124857/1‑249 (MQ=255)     aTCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCt                                                                                                                                                                                                                                     <  2:62445/251‑1 (MQ=255)       cAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGc                                                                                                                                                                                                                                   <  1:73378/251‑1 (MQ=255)           aTGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCgg                                                                                                                                                                                                                               >  1:49447/1‑251 (MQ=255)                aTCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACATCAACAATAACTGCTCCTGTCGTGGGTCATCACTGCGAGGAAGCCTCTTGCCTGCAGGTCTGTCCGGTCAAt                                                                                                                                                                                                                          >  1:276573/1‑251 (MQ=255)                    tttATGCAATATGAAATGCAATGTTTCATCTCATTTTAAAGGAGCCGAAATGCACCGATCTGTGGTGGCAGAACAACTGTGAAGTACAGGATGTAATACCTGTCTCGCTACATGTTCGGGAGTGCATAAAAAGCAAGGTTAACCGCAACACACGCGCCTGAACCTGGCGAAGACGTAAACAATCAATGCCCCTGTCGTGTGTCATAACTGTGAGGAAGCCTCTTGCCTGCGGGTCTGCCAGGTCAATGcc                                                                                                                                                                                                                       >  2:21470/1‑250 (MQ=255)                        tGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGCGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGGTGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCATGTCTGACCGGTCAATGCCAtctc                                                                                                                                                                                                                  >  1:90573/1‑251 (MQ=255)                                 aaaTGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAggg                                                                                                                                                                                                          >  1:69361/1‑250 (MQ=255)                                     gCAAGGTTGCATAGCATTGTCAAGGAGCCGACATGAACCGCTTTGTGGGGGCCGAACCCCTGTGGTGTACAGGATGTAATCCCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGACCTGGCGAAGACGTCAACAATCACTGCCCCTGCCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCACGGTCAATGCCATCTCTCAGAGGgatgat                                                                                                                                                                                                     <  2:129211/251‑1 (MQ=255)                                                               gCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTg                                                                                                                                                                           >  2:24804/1‑251 (MQ=255)                                                                                    atggCCGAACCACTGTGGTGTACATGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGACTCTTGCCTGCATGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATCATGCGATCAAACTCAACGAAAGCCTCTGGATTGGCTGCAAGCTTTgcgc                                                                                                                                                      >  2:207511/2‑251 (MQ=255)                                                                                      ggCCGAACCACTGCGGTTTACAGGATGTAATACCTGTCTCGCTGAATGTTCGCACGTGCATAAAACGCAAGGTTTACAGCACCACCCGCGCCTGGCCCTGGAGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGCGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGCTGATGCGATCCAACTCAACGAAAGACTCTGTATTTGCTGCAAGCTTTGCGCCg                                                                                                                                                    >  2:273614/1‑251 (MQ=255)                                                                                                  gtggtgTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCAt                                                                                                                                        >  1:111492/1‑251 (MQ=255)                                                                                                                       cTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTc                                                                                                                     >  1:280888/1‑249 (MQ=255)                                                                                                                                 gCCTGTGCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTc                                                                                                          <  2:180623/250‑1 (MQ=255)                                                                                                                                  ccTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCg                                                                                                        <  1:24804/251‑1 (MQ=255)                                                                                                                                      atCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTg                                                                                                     >  2:70642/2‑250 (MQ=255)                                                                                                                                                       cGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAatat                                                                                   >  1:164472/1‑251 (MQ=255)                                                                                                                                                         cAAGGTGTACATCAACCCCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCGGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGt                                                                                 <  2:234180/251‑1 (MQ=255)                                                                                                                                                                  cctcaccacCCGCGGCTGGCCCTGTCGAAGACGTCATCAATCGCTGTCCCTGTCGTGTGTTATCCCTGTGAGTACGCCTCTTGCCTGCCGGTCTGCCCGGTCAATGCCATCTCTCAGAGCGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAATCTTTGCGCCGTGGTCTGCTCATTTGGCGCAATCAGCGCTTCAGGAAGCCTTCCGGTGAATGCCCCTGCGCAATATGTTTTTCAGGc                                                                        <  2:204193/245‑1 (MQ=255)                                                                                                                                                                   aGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCt                                                                       >  2:283675/1‑251 (MQ=255)                                                                                                                                                                   aGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCt                                                                       <  2:280888/251‑1 (MQ=255)                                                                                                                                                                      aaCACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGaa                                                                    <  1:64148/251‑1 (MQ=255)                                                                                                                                                                       acacCCGCGCATGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGACTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTNa                                                                     >  2:253920/1‑249 (MQ=255)                                                                                                                                                                          cccGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGct                                                                <  2:111492/251‑1 (MQ=255)                                                                                                                                                                                     ccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGAc                                                     <  1:267158/251‑1 (MQ=255)                                                                                                                                                                                                           cACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGCCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacac                               >  1:190793/1‑251 (MQ=255)                                                                                                                                                                                                              tGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacacaac                            >  1:37923/1‑251 (MQ=255)                                                                                                                                                                                                              tGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacacaac                            >  1:106148/1‑251 (MQ=255)                                                                                                                                                                                                                  ccggTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGc                        <  2:124857/248‑1 (MQ=255)                                                                                                                                                                                                                  ccTGTCGTGTGTCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGc                        >  1:144462/1‑251 (MQ=255)                                                                                                                                                                                                                   ctgtcgcgcgacataaccgtcacgAAGCCTCTTGCCTGAACGTCCGCCCGGTCACTCCCATCTCTCAGAGCGATGATGCGATCCACCCCAACCAAAACCTCTGTATTCGCTGCAAGCTTTGCGCCCTGGTCTGCCCCTTTAGCGCAATCAGCGCTTCAGGACGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCCGCCTGAAGGCTCACTCAACGACGGCGAAGAAAACGCGCCATCACAACATGcc                       >  2:125018/24‑250 (MQ=255)                                                                                                                                                                                                                          gtgtCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTtgctgc                >  1:73396/1‑251 (MQ=255)                                                                                                                                                                                                                           tgtCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTgcg               >  1:2882/1‑251 (MQ=255)                                                                                                                                                                                                                            gtCATCACTGTGAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTgcgc              >  1:189797/1‑251 (MQ=255)                                                                                                                                                                                                                                    tgtgAGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAAcc      >  1:191890/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGGAAGCCTCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGgt  <  1:283675/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   GTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCCGAACCACTGTGGTGTACAGGATGTAATACCTGTCTCGCTGCCTGTTCGGACGTGCATAAAACGCAAGGTTTACAGCAACACCCGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCTGCGCTGGGAACCTGGT  >  NC_000913/2601133‑2601614

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