Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,012,224	C→T	R198Q (CGA→CAA)	metR ←	DNA‑binding transcriptional dual regulator MetR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,012,224	0	C	T	100.0%	315.3 / NA	92	R198Q (CGA→CAA)	metR	DNA‑binding transcriptional dual regulator MetR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/51);  total (41/51)

GCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTG  >  NC_000913/4011978‑4012468                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACt                                                                                                                                                                                                                                                  >  1:496501/1‑251 (MQ=255)     ccGNAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGAACAGTCGGCGCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCGCTACTACCCAATGCGGTAGCGCGGCAATACCCAGCCGCGCGTCAACCATCTGGATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGATTCCTGCCGGCTGAAGAAAGTGCTGCCTGACAGCCAGTTGTCTACGc                                                                                                                                                                                                                                              <  2:599037/251‑1 (MQ=255)      cGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCt                                                                                                                                                                                                                                             >  2:141836/1‑251 (MQ=255)       gTAACTGGCTGACGCTGCTCGCCTTCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGTCCAGGGTATTGGTCTCCACCAGACCCCTGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATGCCCAACCGCGCGGCAACCATCTGAATCACCAATAAGGTGTTATCGACGCTTTTCACTGACGGGCTTACTCCTGCCGGCTGAAGAAAATGCCGTCAGACATCCAGTTGACTACGCTg                                                                                                                                                                                                                                            <  1:305561/251‑1 (MQ=255)       gTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTg                                                                                                                                                                                                                                            <  1:611886/251‑1 (MQ=255)            tGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCg                                                                                                                                                                                                                                       >  1:97209/1‑251 (MQ=255)                tGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGAt                                                                                                                                                                                                                                    >  2:154662/1‑250 (MQ=255)                    gctgctCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAAtt                                                                                                                                                                                                                               >  1:660415/1‑251 (MQ=255)                    gctgctCGCCATCGCGCACGGCGGCGTACAGTCAGCGCCACAAGCCTTCGCCCAGGGTTTTGGTCGCCACCAGACCCTGGCGCTCAAAACGCTCTACCACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGCCGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACAGCCAGTTGACTAAGCTGCACCGGATAAAtt                                                                                                                                                                                                                               <  2:191994/251‑1 (MQ=255)                     ctgctCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAAtt                                                                                                                                                                                                                               >  1:548389/1‑250 (MQ=255)                      tgctCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACCTCCAGTTGACTACGCTGCACCGGATAAATTa                                                                                                                                                                                                                              >  2:415732/1‑250 (MQ=255)                       gctCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAAt                                                                                                                                                                                                                            <  1:159749/251‑1 (MQ=255)                          cGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGc                                                                                                                                                                                                                         <  1:454647/251‑1 (MQ=255)                           gCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCg                                                                                                                                                                                                                        >  1:159497/1‑251 (MQ=255)                              aTCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACACATCCGCGCGGCAACCATCTGAATCAACAATACGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCAGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGtct                                                                                                                                                                                                                     >  2:167804/1‑251 (MQ=255)                                   gcACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCTCCACCATACCCTGGCGCTCAAAACTCTCTACTACCCAATTCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGTTGTTATCGACGCTTTTCAGTGACTGTCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTg                                                                                                                                                                                                                <  2:50191/251‑1 (MQ=255)                                   gcACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCAGCCGCGCGTCAACCATCTGAATCATCAATAAGGTGTTATCGACGCTTTGCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATGAATAGCGTCTCGCTg                                                                                                                                                                                                                <  1:471056/251‑1 (MQ=255)                                         cggcgTACAGTCGGCTCCACAAGCCTTCGCCTAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATCCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTACCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCgaga                                                                                                                                                                                                          <  2:245281/251‑1 (MQ=255)                                              tACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTc                                                                                                                                                                                                     >  1:213982/1‑251 (MQ=255)                                                 aGTCGGCTCCACAAGCCTTCTCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCgg                                                                                                                                                                                                  <  1:293148/251‑1 (MQ=255)                                                  gTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCggg                                                                                                                                                                                                 >  1:506271/1‑250 (MQ=255)                                                        cccACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAAtt                                                                                                                                                                                           <  2:150968/250‑1 (MQ=255)                                                         ccACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTc                                                                                                                                                                                          <  1:34996/251‑1 (MQ=255)                                                          cacaAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTcg                                                                                                                                                                                         <  1:228560/251‑1 (MQ=255)                                                          cacaAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTcg                                                                                                                                                                                         <  1:320177/251‑1 (MQ=255)                                                             aaGCCTTCGGCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTAAAAGCTCTCTACTACCCAATGCGGTAGCGCGGCAATTCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGt                                                                                                                                                                                      <  2:417707/251‑1 (MQ=255)                                                               gCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGttt                                                                                                                                                                                    >  2:560565/1‑251 (MQ=255)                                                                        cAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCa                                                                                                                                                                           <  2:159497/251‑1 (MQ=255)                                                                         aGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCgccgcc                                                                                                                                                                            >  2:321947/1‑249 (MQ=255)                                                                          gggTTTTGGTCCCCACCCGACCCTGGCGCTAAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGt                                                                                                                                                                         <  2:472740/251‑1 (MQ=255)                                                                              tttGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGa                                                                                                                                                                      <  2:660415/250‑1 (MQ=255)                                                                                tggccaccgcCAGACCCTGGCGCTCAAACCTCTCTACTACCCAATGCGGTAGCGCGGCAATAGCCATCCGCGGGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATATCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCTCGTTTTCGCCGCCAGTGGATg                                                                                                                                                                    <  1:419993/246‑1 (MQ=255)                                                                                 ggTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGGGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGt                                                                                                                                                                  >  2:188184/1‑251 (MQ=255)                                                                                  gTCACCACCAGCCCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATTCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTc                                                                                                                                                                 <  1:169608/251‑1 (MQ=255)                                                                                    caccacCAGACCCTGGCACTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGAGTTTTCGCAGATTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGTTCTTCCGGTGTCACTCTCGTTTTCGCCGCCAGTGGCTGGNc                                                                                                                                                                 >  2:587478/1‑249 (MQ=255)                                                                                        accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGc                                                                                                                                                           >  2:294415/1‑251 (MQ=255)                                                                                        accaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGc                                                                                                                                                           <  1:343591/251‑1 (MQ=255)                                                                                         ccaGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCt                                                                                                                                                          <  2:349692/251‑1 (MQ=255)                                                                                             aCCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAac                                                                                                                                                       >  2:342348/1‑250 (MQ=255)                                                                                               ccTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACAcc                                                                                                                                                    <  1:79030/251‑1 (MQ=255)                                                                                                 tGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACAcc                                                                                                                                                    >  1:218944/1‑249 (MQ=255)                                                                                                    cgcTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACg                                                                                                                                               <  1:655198/251‑1 (MQ=255)                                                                                                       tCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGGCAGGCGCTAACACCAAACGCAc                                                                                                                                            >  2:6539/1‑251 (MQ=255)                                                                                                           aaCTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCa                                                                                                                                        <  1:392389/251‑1 (MQ=255)                                                                                                              tctctACTCCCCAATGCGGCAGCGAGGCAATACCCATCCGCGCGGAAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATACATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCAt                                                                                                                                       >  2:310987/1‑249 (MQ=255)                                                                                                              tctctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCAt                                                                                                                                       >  2:567007/1‑249 (MQ=255)                                                                                                                tctACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTc                                                                                                                                   >  2:94071/1‑251 (MQ=255)                                                                                                                  tactacCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCg                                                                                                                                  >  1:336555/1‑250 (MQ=255)                                                                                                                      cccaaaTGCGGTAGCGCGTCAATCCCCACCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGGCGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTCCTCGGCACCAGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCTCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGaa                                                                                                                                <  1:407020/244‑1 (MQ=255)                                                                                                                          aaTGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGc                                                                                                                         >  2:201404/1‑251 (MQ=255)                                                                                                                               ggTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACACTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGact                                                                                                                     >  2:235526/1‑248 (MQ=255)                                                                                                                                  aGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATg                                                                                                                 <  1:102356/251‑1 (MQ=255)                                                                                                                                  aGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATg                                                                                                                 <  1:196999/251‑1 (MQ=255)                                                                                                                                   gcgcGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGc                                                                                                                <  1:458687/251‑1 (MQ=255)                                                                                                                                         cAATACCAATCCGCGAGGCAACAATCTGAATAAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTACAGGCTGAAGAAAATGCCGACAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGAGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGAAGCCAGGGGATGGTCAGGAGCTAACACCAGACGCCCTACATAGTCGAACATCGGCGATTAATGAAGGcc                                                                                                           >  2:474520/1‑250 (MQ=255)                                                                                                                                           aTACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACt                                                                                                        <  1:52512/251‑1 (MQ=255)                                                                                                                                           aTACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGCGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGCAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGAAGGCCACt                                                                                                        >  2:63899/1‑251 (MQ=255)                                                                                                                                             aCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTgc                                                                                                      <  1:560565/251‑1 (MQ=255)                                                                                                                                             aCCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTgc                                                                                                      >  2:48675/1‑251 (MQ=255)                                                                                                                                                       gcgGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAat                                                                                            <  1:631883/251‑1 (MQ=255)                                                                                                                                                         ggCAAGCATCCGAATCCACAATCAGGCGCTATCGACACTTTTCACTGACGGGCGTACGCCTGCCGGCTGACGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCGGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACAACAGCCGCAATTCATAGTCGAACATCGGCGAATCATGCATGCCACTGCGCGCCAGActat                                                                                          >  2:116830/1‑251 (MQ=255)                                                                                                                                                          gCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATc                                                                                         <  2:97209/251‑1 (MQ=255)                                                                                                                                                             aCCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGa                                                                                      <  2:432909/251‑1 (MQ=255)                                                                                                                                                              ccATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAACTGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTAATAGTAGAACAACGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCNGAc                                                                                     >  2:609728/1‑251 (MQ=255)                                                                                                                                                                  cTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTc                                                                                  >  1:49543/1‑250 (MQ=255)                                                                                                                                                                          acaaTAAGGTGTGATCGACGCTTGTCAGTGACGGGCTGACGCCTGCCGGGTGAAGAAAATGCCTCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAAGAAGGCAGGCCACTGCGCGTCAGGATATCGGGCGTCATTGCCAGa                                                                         <  2:308860/251‑1 (MQ=255)                                                                                                                                                                             aTAAGGTGTTATCGACCCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATcc                                                                      <  2:402931/251‑1 (MQ=255)                                                                                                                                                                               aaGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCTGCAGAATATCGGACGTCATTACCAGATCCAg                                                                    <  1:448393/251‑1 (MQ=255)                                                                                                                                                                               aaGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCa                                                                     <  2:620054/250‑1 (MQ=255)                                                                                                                                                                                  gtgtTATCGACGCTTTTCAGCGACGGCCTGACGCCTGCCGGCTGAAGACAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGAGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGctc                                                                 >  2:159969/1‑251 (MQ=255)                                                                                                                                                                                   tgtTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGCAAAAGCAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGctct                                                                >  1:454041/1‑251 (MQ=255)                                                                                                                                                                                          gACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGtt                                                         >  1:578850/1‑251 (MQ=255)                                                                                                                                                                                           ccGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGAGCCAGCTCTCCCTGTTg                                                        <  2:218944/250‑1 (MQ=255)                                                                                                                                                                                                       gACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCt                                             <  1:294415/250‑1 (MQ=255)                                                                                                                                                                                                         cGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCg                                          <  2:500020/251‑1 (MQ=255)                                                                                                                                                                                                             cTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTc                                      >  2:442342/1‑251 (MQ=255)                                                                                                                                                                                                                  gCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATg                                 >  1:490646/1‑251 (MQ=255)                                                                                                                                                                                                                    cTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTc                               >  1:38139/1‑251 (MQ=255)                                                                                                                                                                                                                          gCTGAAGAAAAGGCCGCCCGACATCCAGTTGACTACGCTGCAACGGATAAATTAATAGCGTCGCGGTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCCGTGGATGGACAGGTGCTAACACCAGACGCGCTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTATCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTTGGGGTCAAATGTCACGccc                         <  2:225414/251‑1 (MQ=255)                                                                                                                                                                                                                             gaagaaAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGa                       <  1:181788/250‑1 (MQ=255)                                                                                                                                                                                                                             gaagaaAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGAt                      >  1:231210/1‑251 (MQ=255)                                                                                                                                                                                                                               agaaAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGAttt                    <  1:240089/251‑1 (MQ=255)                                                                                                                                                                                                                               agaaAATGCCGCCAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCCGGTCAAATGTCACGCCCGAttt                    <  1:516787/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaaaTGCCGCCAGACATCCAGTTGACTACTCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCATCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGAGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACAAGATCCAGCTCTCCCTGTCGCAAGGCGGGCTGCGGGTCCAATGTCACGCCCGATTTa                   >  2:198185/1‑250 (MQ=255)                                                                                                                                                                                                                                        cgccAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCAt           <  2:445032/251‑1 (MQ=255)                                                                                                                                                                                                                                         gccAGACATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCCCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCAtc          <  1:494922/251‑1 (MQ=255)                                                                                                                                                                                                                                          ccAGACATCCAGGTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCAtct         >  2:50122/1‑251 (MQ=255)                                                                                                                                                                                                                                              aCATCCATTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGAGGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTAc     <  1:589685/251‑1 (MQ=255)                                                                                                                                                                                                                                              aCATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAgg                                                   <  2:125364/205‑1 (MQ=255)                                                                                                                                                                                                                                              aCATCCAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAgg                                                   >  1:125364/1‑205 (MQ=255)                                                                                                                                                                                                                                                  ccAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTg  <  1:46659/250‑1 (MQ=255)                                                                                                                                                                                                                                                  ccAGTTGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTg  <  1:641172/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGACTACGCTGCACCGGATAAATTAATAGCGTCTCGCTGGCGAGATCTTCCGGTGTAATTCGCGTTTTCGCCGCCAGTGGATGGTCAGGTGCTAACACCAGACGCACTTCATAGTCGAACATCGGCGAATAATGCAGGCCACTGCGCGGCAGAATATCGGACGTCATTACCAGATCCAGCTCTCCCTGTTGCAAGGCGGGCTGCGGGTCAAATGTCACGCCCGATTTAAAATCCATCTCTACCTG  >  NC_000913/4011978‑4012468

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