Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,185	C→T	T156T (ACG→ACA)	folE ←	GTP cyclohydrolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,185	0	C	T	98.9%	337.6 / NA	88	T156T (ACG→ACA)	folE	GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (43/44);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGA  >  NC_000913/2242938‑2243424                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                cgcataaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTc                                                                                                                                                                                                                                              <  1:166533‑M2/244‑1 (MQ=255)      aaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACg                                                                                                                                                                                                                                         >  1:280732/1‑251 (MQ=255)       aTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCTCAGCGCGCAGAACCTCGTGGCGCGCATTCTGACTGGATTTGAACACTCCCCCAAGCGAGGTCGTTGTCGTGGCCCTGGTTGCATCGCGGATGCTCCGCGCCTTCACGCAGTAATGCACCGCGTCGTTCGAGGCAGCCACGTTATTGGTGCCCAGCCGCGTTTGTCGCGCAATAAGACTTTGCTGTGTCAGatgt                                                                                                                                                                                                                                        >  2:356492/1‑248 (MQ=255)         cGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGtt                                                                                                                                                                                                                                       <  1:637411/250‑1 (MQ=255)          gAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCTTATTCTGACTGGATTTGAACAATCCACCAAGAGAGTTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGGGACAGCCACGTTATTTGTGCCCAGCAGCGGTTGTAGCGCAATAAGAAGTGGCTGTGTCATACGTTcc                                                                                                                                                                                                                                     <  2:385567/251‑1 (MQ=255)                gACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCAcc                                                                                                                                                                                                                               >  1:358945/1‑251 (MQ=255)                   gTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGCCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGc                                                                                                                                                                                                                            >  1:70919/1‑251 (MQ=255)                       cgcTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGAc                                                                                                                                                                                                                        <  2:70919/251‑1 (MQ=255)                           ccATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCt                                                                                                                                                                                                                     <  1:593944/250‑1 (MQ=255)                             attgTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTggg                                                                                                                                                                                                                  <  1:578214/248‑1 (MQ=255)                             attgTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCATTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTggg                                                                                                                                                                                                                  >  2:184021/4‑251 (MQ=255)                                  tCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCaa                                                                                                                                                                                                               >  2:536744/1‑249 (MQ=255)                                  tCCTGCCTTTTAACCAGTTGAGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGACGTAATGCACCGCGTCGATCGAGACAGCAACGTTCTTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCCGCGGACGATGGGCaa                                                                                                                                                                                                               >  2:425550/1‑249 (MQ=255)                                   ccTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGAGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCCGGGCaa                                                                                                                                                                                                               >  1:389368/1‑248 (MQ=255)                                      gCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACt                                                                                                                                                                                                         >  2:323297/1‑251 (MQ=255)                                                cAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTGCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCg                                                                                                                                                                                               <  2:84384/251‑1 (MQ=255)                                                 aGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCgg                                                                                                                                                                                              <  1:417060/251‑1 (MQ=255)                                                  gTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGAAAAGAACTGCACAATGCGGt                                                                                                                                                                                             >  2:214609/1‑251 (MQ=255)                                                     gtgATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTa                                                                                                                                                                                           >  1:37494/1‑250 (MQ=255)                                                     gtgATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTa                                                                                                                                                                                           >  2:345885/1‑250 (MQ=255)                                                        atgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAAGCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAttt                                                                                                                                                                                       <  2:162020/251‑1 (MQ=255)                                                          gaCGCACAGCGCGCAGACACTCGTGGCACGTATTCTGACTGGATTTGAACAATCCACCAAGCGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGACGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAttttt                                                                                                                                                                                     >  2:449044/1‑251 (MQ=255)                                                           aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTg                                                                                                                                                                                    >  2:180639/1‑251 (MQ=255)                                                               acaGCGCGCAGAAACTCGTGGCGCGTAGTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGAGGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGGTTGTATAGTAATAAGAATTTGCTGTGTCAGACGTTCTTGCCCCTGCGGACGGTGGGCAAAGAACTGCACAATTCTGTTAATTTTTGACAg                                                                                                                                                                                <  1:643084/251‑1 (MQ=255)                                                                 aGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTcc                                                                                                                                                                                                                                     <  1:571028/196‑1 (MQ=255)                                                                 aGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTcc                                                                                                                                                                                                                                     >  2:571028/1‑196 (MQ=255)                                                                    gcgcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCg                                                                                                                                                                            >  2:136316/1‑250 (MQ=255)                                                                        cgAAACGCGGGGCGCGTATTCTGACTGGATGTGAACAACCCACCAAGAAAGGTCGTGGGCGTCGCACTGGTTGCATCGCTGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGGTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCAc                                                                                                                                                                       <  2:600891/250‑1 (MQ=255)                                                                        aGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAAACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCAc                                                                                                                                                                       <  2:566503/251‑1 (MQ=255)                                                                        aGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCAATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACGGCACAATGCGGTTAATTTTTGACAGACCGATCAc                                                                                                                                                                       >  2:560096/1‑251 (MQ=255)                                                                         gAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCAcc                                                                                                                                                                      >  1:234081/1‑251 (MQ=255)                                                                          aaaCTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCCATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCCGGGCAAAGATCTGCACAATGCGGTTAATTTTTGACAGACCGATCACCg                                                                                                                                                                     >  2:610229/1‑251 (MQ=255)                                                                          aaaCTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCg                                                                                                                                                                     <  1:432171/251‑1 (MQ=255)                                                                           aaCTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCg                                                                                                                                                                     >  2:8639/1‑250 (MQ=255)                                                                                gTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCtt                                                                                                                                                               >  1:281430/1‑251 (MQ=255)                                                                                  ggCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCTTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTc                                                                                                                                                             <  1:214609/251‑1 (MQ=255)                                                                                  ggCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTc                                                                                                                                                             >  2:284852/1‑251 (MQ=255)                                                                                   gcgcGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCg                                                                                                                                                            <  1:128206/251‑1 (MQ=255)                                                                                    cgcgTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCgg                                                                                                                                                           <  1:248776/251‑1 (MQ=255)                                                                                          ttCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGTATAGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGatat                                                                                                                                                      >  2:506837/1‑250 (MQ=255)                                                                                              gACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGAGCACCGAATCTTTCGGGATATAgg                                                                                                                                                   <  2:281430/249‑1 (MQ=255)                                                                                                   gagttGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACGGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAAGAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACAGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGt                                                                                                                                            <  1:356492/248‑1 (MQ=255)                                                                                                        gAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCt                                                                                                                                        <  1:536744/250‑1 (MQ=255)                                                                                                          aCAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTg                                                                                                                                     >  2:642080/1‑251 (MQ=255)                                                                                                             aTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                   >  1:369477/1‑250 (MQ=255)                                                                                                                    aaGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGACATAGGCCACCGTCGCTTTGCCATCGATgg                                                                                                                           <  1:136316/251‑1 (MQ=255)                                                                                                                         agGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAAc                                                                                                                       >  2:209463/1‑250 (MQ=255)                                                                                                                           gTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACaaa                                                                                                                    <  2:280732/251‑1 (MQ=255)                                                                                                                              gTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATg                                                                                                                 <  2:37494/251‑1 (MQ=255)                                                                                                                                tGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGgt                                                                                                               >  2:164134/1‑251 (MQ=255)                                                                                                                                   cGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGgtgt                                                                                                             <  1:284852/250‑1 (MQ=255)                                                                                                                                    gTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTc                                                                                                           <  2:107263/251‑1 (MQ=255)                                                                                                                                      ggCACTGGTTGCATCGCGGACGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTcac                                                                                                         <  1:449044/251‑1 (MQ=255)                                                                                                                                            ggTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTg                                                                                                    <  1:323297/250‑1 (MQ=255)                                                                                                                                                gCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGt                                                                                               >  2:595133/1‑251 (MQ=255)                                                                                                                                                gCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGt                                                                                               <  1:168320/251‑1 (MQ=255)                                                                                                                                                gCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGt                                                                                               <  1:425550/251‑1 (MQ=255)                                                                                                                                                  aTCGCGGATGCTACGCGCCTTCACGCCGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACACTGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTACCAAAATGGTGTTCACAGGTGCTGGTc                                                                                              >  1:637786/1‑250 (MQ=255)                                                                                                                                                          tGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGat                                                                                      >  1:164791/1‑250 (MQ=255)                                                                                                                                                            ctaCGCGCATTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTCGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACCAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATc                                                                                   >  2:247970/3‑251 (MQ=255)                                                                                                                                                            ccACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATc                                                                                   <  1:180639/251‑1 (MQ=255)                                                                                                                                                                gcgcCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCgcgc                                                                               >  1:472128/1‑251 (MQ=255)                                                                                                                                                                gcgcCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCATTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCgcgc                                                                               <  1:184021/251‑1 (MQ=255)                                                                                                                                                                    cTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACgg                                                                           >  2:187278/1‑251 (MQ=255)                                                                                                                                                                                  gCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATc                                                              <  2:369477/250‑1 (MQ=255)                                                                                                                                                                                    cccGCGTCGTTCGAGACAGCCACTGTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACATTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGAc                                                           <  2:421367/250‑1 (MQ=255)                                                                                                                                                                                      cgcgTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCAGTTCATCGACCt                                                         <  1:209463/251‑1 (MQ=255)                                                                                                                                                                                          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAt                                                     >  2:189541/1‑251 (MQ=255)                                                                                                                                                                                               cgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGt                                                >  1:182641/1‑251 (MQ=255)                                                                                                                                                                                                  gaCAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAGCAAAATGGTGTTCACAGGTGCTGTTCAGAGTGATAGCGCGCACGGTGACCATTTCATCGACCTTCATTTTGtttt                                             <  1:86920/251‑1 (MQ=255)                                                                                                                                                                                                     aGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGACTTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGCCCTTCATTTTGTTTTCa                                           >  1:34910/1‑250 (MQ=255)                                                                                                                                                                                                     aGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCa                                           >  1:414030/1‑250 (MQ=255)                                                                                                                                                                                                     aGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAAGCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCaa                                          <  1:635821/251‑1 (MQ=255)                                                                                                                                                                                                       ccACGTTATTGGTGCCCAGCATCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCGCGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATg                                        <  1:506837/251‑1 (MQ=255)                                                                                                                                                                                                        cACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGa                                       <  1:642080/251‑1 (MQ=255)                                                                                                                                                                                                           gTTATTGGTGCTCAGCAGCGTTTGTAGCGAAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCATTGAggg                                    >  2:369436/1‑251 (MQ=255)                                                                                                                                                                                                                tGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGAtt                               >  1:545359/1‑251 (MQ=255)                                                                                                                                                                                                                 ggTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCGTCAGTTTGTTTTCAATGAGGGTGAtt                               <  1:637694/250‑1 (MQ=255)                                                                                                                                                                                                                     cccAGAAGCGTTTGAAGGCCAATCAGATTTTGCTGTGTAAGCAGTTACTGCACCTGCGCACGCTGGGAAAATAACTGCACAATCCGGTTAATTTTTGACAGACAGATAAACGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCTATAGCTGGTCACACAATGGTGTTCACAGGTGCTGGTCCGCGTGATATCGCGCACGGTGACCCTTTCATCGACATTCCATTCGTTCTCACTGAGGGTGGTTTTCgg                          >  2:670409/1‑251 (MQ=255)                                                                                                                                                                                                                      ccAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCggg                         >  2:591619/1‑251 (MQ=255)                                                                                                                                                                                                                       cagcagCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGa                        <  1:164134/251‑1 (MQ=255)                                                                                                                                                                                                                          cagcGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAAt                     <  1:554997/251‑1 (MQ=255)                                                                                                                                                                                                                          cagcGTTTGTAGCGCAAGAAGAATTCGCTGTGGCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCCCCGCATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGATGGTGATTTTCGGTaaag                     <  2:210079/251‑2 (MQ=255)                                                                                                                                                                                                                                   tAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAAt            >  1:284001/1‑251 (MQ=255)                                                                                                                                                                                                                                       gcAATCAGAATTTGCTGCGACAGACGTTCCTGAACCTGCGGTCGCTGGGCAACGATTTGCACAATGCGGTTAATTTTTGCCAGACCGATCACCGAATCTTTCGGGCTATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATcc          >  2:367717/1‑249 (MQ=255)                                                                                                                                                                                                                                             aGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGga  <  1:189541/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaTTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGa                                 >  2:505443/1‑218 (MQ=255)                                                                                                                                                                                                                                               aaTTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGa                                 <  1:505443/218‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                AACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGA  >  NC_000913/2242938‑2243424

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