Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,469	C→T	G633D (GGC→GAC)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,469	0	C	T	98.8%	272.2 / ‑6.3	86	G633D (GGC→GAC)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/42);  minor base G (0/1);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGT  >  NC_000913/4297222‑4297712                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCtt                                                                                                                                                                                                                                                  >  1:129495/1‑251 (MQ=255)    cAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTCCTCCGGCGTTTTCGTAATCGGGCTTACGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGGGATACTTTTTTCTTTAc                                                                                                                                                                                                                                               >  1:450967/1‑251 (MQ=255)     aGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTAcg                                                                                                                                                                                                                                              <  2:550010/251‑1 (MQ=255)       ggCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTAcgcg                                                                                                                                                                                                                                            <  2:252921/250‑1 (MQ=255)               cgcgCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCAc                                                                                                                                                                                                                                     <  1:718885/250‑1 (MQ=255)               cgcgCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAAGTTTGTCTTTACGCGAGTGCAc                                                                                                                                                                                                                                     <  2:725181/250‑1 (MQ=255)                     ggCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAAcc                                                                                                                                                                                                                              >  1:360922/1‑251 (MQ=255)                      gCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCa                                                                                                                                                                                                                             <  2:77610/251‑1 (MQ=255)                        gAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGCCGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAAt                                                                                                                                                                                                                           >  1:683120/1‑251 (MQ=255)                                 aCTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCCGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCt                                                                                                                                                                                                                  >  1:275685/1‑251 (MQ=255)                                    tggtgtacccggtaaTCTCGTACTGCTCGGCGGCGCGCTGGTCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGGACTCCGGCGTTTTCGGAATCGGGCTTAAGTTTTCGGTAACCAGCTCGGTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCGTTACGCGAGTGTACCCAAACCAATGCCTCATCtt                                                                                                                                                                                                                 <  1:564321/236‑1 (MQ=255)                                     gttgtACTCGTCAATCCCGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTGTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCa                                                                                                                                                                                                               <  2:123444/250‑1 (MQ=255)                                     gttgtACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCaa                                                                                                                                                                                                              <  1:184033/251‑1 (MQ=255)                                     gttgtACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCa                                                                                                                                                                                                               <  1:123148/250‑1 (MQ=255)                                          aCTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCTCCCAAACCAATGCCTCATCTTCAATAccc                                                                                                                                                                                                         <  2:419892/251‑1 (MQ=255)                                             cGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATAccc                                                                                                                                                                                                         >  1:39428/1‑248 (MQ=255)                                              gTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCCCCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCCATGCCTCATCTTCAATACCCAGAc                                                                                                                                                                                                     >  1:40474/1‑251 (MQ=255)                                                cAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGt                                                                                                                                                                                                   >  1:333906/1‑251 (MQ=255)                                                 aaTCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGtt                                                                                                                                                                                                  >  1:441659/1‑251 (MQ=255)                                                  aTCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGGTAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCCATACCCAGACGttt                                                                                                                                                                                                 >  1:157741/1‑251 (MQ=255)                                                   tCACGCACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGCCGCGAACAGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGCTTTCGGTACCCCGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGCGGGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTg                                                                                                                                                                                                >  2:432021/1‑251 (MQ=255)                                                     aCGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGc                                                                                                                                                                                              >  1:302862/1‑251 (MQ=255)                                                      cGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGGTTGGCg                                                                                                                                                                                             <  2:452841/251‑1 (MQ=255)                                                         ccTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTc                                                                                                                                                                                           <  2:267558/249‑1 (MQ=255)                                                              tCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGt                                                                                                                                                                                     >  1:150325/1‑251 (MQ=255)                                                               ctgcggcgCGCTGAGTGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAAGCGGGCTTAAGTTTTCGGTAACCAGCGCGTTACAGGCACCAATCCACCACTGGGAGGTCATGTAAATCGCACCGTTGTTCGGACGAGCGCTGACCTGCGCACGGGTGATAAGTGTGGCTTTACGCGAGTGCACACAAACCAATGCCTCATCTCCAATACCCAGACGTTTGGCgggtgcgg                                                                                                                                                                                      <  2:325131/247‑8 (MQ=255)                                                                     cgcgcTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGAtt                                                                                                                                                                              <  2:524522/251‑1 (MQ=255)                                                                     cgcgcTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTAATCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATACACCACTGGTAGGTCATGAAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGAGATAATTTTGACTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGAtt                                                                                                                                                                              >  1:352572/1‑251 (MQ=255)                                                                        gcTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTtgt                                                                                                                                                                           >  1:433788/1‑251 (MQ=255)                                                                         cTGATCGGCGATCGGCTCGAAGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTtgtg                                                                                                                                                                          >  2:81696/1‑251 (MQ=255)                                                                                     cggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGtt                                                                                                                                                              >  1:260952/1‑251 (MQ=255)                                                                                          cGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTtcatc                                                                                                                                                          <  1:250567/250‑1 (MQ=255)                                                                                            aCGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATcagc                                                                                                                                                       <  1:81696/251‑1 (MQ=255)                                                                                            aCGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGAAATCGGGCTTAAGTTATCGGTAACCTGCTCGTTACAGGCACAAATCCACAACAGGCAGGTCATGTAAATGGACCCTTTGTTCGGACGATAGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGGGCACCCAACCCAATGCCTCATCTTCAATACCCAGCAGTTAGGCGTCTTCGGTATTGATTTGTGCGACGCCAGGTTCATcagc                                                                                                                                                       >  2:263040/1‑251 (MQ=255)                                                                                                          gTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCGGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTGCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTg                                                                                                                                         <  2:360922/251‑1 (MQ=255)                                                                                                          gTATTTGTACTCCGGCGTTTTCGTAAACGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCGGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTGGGCGTCTTCGGTATTGATTTGTGCGTAGCCATGTTGATCAGCCAGCGCCGCCAGTg                                                                                                                                         <  2:157741/251‑1 (MQ=255)                                                                                                              ttttACTCCGGCGGTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTGGTCTTTACGCGAGTGCACCCAAACCAATTCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGGGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCg                                                                                                                                      <  2:433788/247‑1 (MQ=255)                                                                                                                 tACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGcac                                                                                                                                   >  1:628913/1‑250 (MQ=255)                                                                                                                               cGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCg                                                                                                                    <  2:129495/251‑1 (MQ=255)                                                                                                                               cGTAAGCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATc                                                                                                                     <  2:39428/250‑1 (MQ=255)                                                                                                                                   aTCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCCGTGCCGCACAGTTACCGGTCATCGAACg                                                                                                                >  2:495361/1‑251 (MQ=255)                                                                                                                                   aTCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTGCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACg                                                                                                                <  1:229892/251‑1 (MQ=255)                                                                                                                                    acGGGCTTAAGTTTTCGGTAACCAGTTCGTTACAGGCCCCAATCCACCACTGGTAGGTCATGTAAATCGTCCCTTTGTTCGGACGATCTCTGCCCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTTAAGACCCAGACGTTTGGCGTCTGCGGTATTGATTTGTGCGTAGCTAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACgg                                                                                                               <  2:120227/250‑1 (MQ=255)                                                                                                                                       ggCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCa                                                                                                             >  1:202019/1‑250 (MQ=255)                                                                                                                                        gCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCCCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAag                                                                                                           <  1:590665/251‑1 (MQ=255)                                                                                                                                        gCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGCTCGCTGCCCTGAGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGCATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACCGTTACCGGACATCGacaggca                                                                                                             >  1:606870/1‑243 (MQ=255)                                                                                                                                            aaGTTTTCGGTAACCAGCTCGGTACAGGCACCAAGCCACCACTTGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAGGGGTGAGAATTTGGTCTTGACGCGAGTGCATCCAAACCAAGTCCTCATCTTCAAAACCCAGACGTTGGGCGTCTTCTGTATGGATTTGTGCGGAGCCATGTTCATCAGCATGCGCCGCCAGGGCCGCACAGTTACCGGTCATCGAACGGCGAGagta                                                                                                       <  2:634516/251‑1 (MQ=255)                                                                                                                                                   cGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTCCGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGAcc                                                                                                 >  1:653842/1‑250 (MQ=255)                                                                                                                                                    ggTAACAAGCTTGTTACAGGCACCACGCCACCACTGTGAGGTCAGGTAAATCGCCCCTTGGTTCGGCCGGTCGCGGTCCTGCTCACGGGTGATAATTTTTTCTTTACGCGAGTGCACCCAAACCAAGGCCTCATCTTCAATCCCAAGACTTTTGCCGTCTTCGGTATTGATTTGGGCGGAGCCAGGTTCATCAGCAAGCGCCGCCAGTGCCGCACAGTGACCGGTCATCGAACGGCAAGAGTAGTGACCaa                                                                                               <  2:606870/251‑1 (MQ=255)                                                                                                                                                       aaCCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTa                                                                                                                                                                                                                                                >  2:563387/1‑103 (MQ=255)                                                                                                                                                       aaCCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTa                                                                                                                                                                                                                                                <  1:563387/103‑1 (MQ=255)                                                                                                                                                        aCCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTc                                                                                           >  2:209875/1‑251 (MQ=255)                                                                                                                                                           aGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACg                                                                                        <  1:390327/251‑1 (MQ=255)                                                                                                                                                           aGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACg                                                                                        >  2:161380/1‑251 (MQ=255)                                                                                                                                                                gTTACAGGCACAAATCCCCCACGGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTGTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCGGGTTCATCAGCCATCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCg                                                                                   <  2:366045/251‑1 (MQ=255)                                                                                                                                                                gTTACAGGCAACAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCg                                                                                   >  2:628190/1‑251 (MQ=255)                                                                                                                                                                     aGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTCCGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGAc                                                                              >  2:713026/1‑251 (MQ=255)                                                                                                                                                                       gCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGAc                                                                              >  1:287601/1‑249 (MQ=255)                                                                                                                                                                              cccaccaCTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTGGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGCCGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATc                                                                     <  1:345102/250‑1 (MQ=255)                                                                                                                                                                                  ccacTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGt                                                                 >  1:535974/1‑251 (MQ=255)                                                                                                                                                                                     cTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACt                                                              >  1:478622/1‑251 (MQ=255)                                                                                                                                                                                       ggTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcg                                                            <  1:593165/251‑1 (MQ=255)                                                                                                                                                                                        gTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCTCGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATCCCCAGACGTTTGGCTTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgt                                                           <  1:119763/251‑1 (MQ=255)                                                                                                                                                                                        gTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgt                                                           >  2:42681/1‑251 (MQ=255)                                                                                                                                                                                          aGGTCATGTAAATCGCCCTTTTGTTCGGACGATCGCTGACCTGCGCACGGGTTATAATTTTGTCTTTACGCGAGGGCACCCAAACCAATGCATCATCTTAAATACCCGGACGGTTGGCGGCTGCGGTAGTGATTTGGGCGTAGCAAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgtc                                                          <  2:420536/250‑1 (MQ=255)                                                                                                                                                                                          aGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgtcg                                                         >  2:580613/1‑251 (MQ=255)                                                                                                                                                                                              cATGTAAATCGCCCCTTTGTTCGGACGGTCGCTGACCGGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCCCCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTGCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                                                      <  2:683120/250‑1 (MQ=255)                                                                                                                                                                                               aTGTAACTCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                                                      >  2:657837/1‑249 (MQ=255)                                                                                                                                                                                                    aaTCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGt                                               <  1:580613/251‑1 (MQ=255)                                                                                                                                                                                                        gCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGCCAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAt                                           >  2:4405/1‑251 (MQ=255)                                                                                                                                                                                                                gTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCt                                    >  1:446983/1‑250 (MQ=255)                                                                                                                                                                                                                     gACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAg                              <  1:161380/251‑1 (MQ=255)                                                                                                                                                                                                                      acaaTCGCTGCCCTGCGCACGGGTGATAATTTTGTCTTTCCGCGAGTGCACCCAACCAAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGt                             <  2:69393/248‑1 (MQ=255)                                                                                                                                                                                                                       cGATCGCTGTCTTGCGCACGGGTTATAATTTTGTCTTTCCGCGAGGGCACCCAAACCAATGCCTCACCTTCAATGCCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTc                            <  2:542836/251‑1 (MQ=255)                                                                                                                                                                                                                        gATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCg                           >  2:532090/1‑251 (MQ=255)                                                                                                                                                                                                                              tGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTg                     >  1:155887/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAACCCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAc                                  <  1:586898/235‑1 (MQ=255)                                                                                                                                                                                                                                 ccTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAc                                  >  2:586898/1‑235 (MQ=255)                                                                                                                                                                                                                                      gcACGGGTGATAATTTTGTCTTTACGCTGGTGCCCCCAAACCAATGCCTCATCTTAAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGTTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTGTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTg             <  2:532452/251‑1 (MQ=255)                                                                                                                                                                                                                                          gtgTGATAATTTTGTCTTGACGGGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGGCGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCACCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGGGAAGAACTGCGcc         <  1:231402/249‑1 (MQ=255)                                                                                                                                                                                                                                           ggTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTCGGCGTCTTCGGTATTGATTTGAGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGTTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGcc         >  1:19030/1‑250 (MQ=255)                                                                                                                                                                                                                                            gTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAg       <  1:209875/251‑1 (MQ=255)                                                                                                                                                                                                                                              gATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGAc     >  1:147998/1‑251 (MQ=255)                                                                                                                                                                                                                                               aTAATTTTGTCGTTCCGCGGGTGCACCCAAACCAATGCCTCATCTTCAATCCCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCTCAGGTGAAGAACTGCGCCAGAcc    <  2:275685/251‑1 (MQ=255)                                                                                                                                                                                                                                                tAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCg   >  2:601587/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaTTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGGCATCGAACGGCAAGAGTAGGGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCgg  >  1:206231/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGT  >  NC_000913/4297222‑4297712

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