Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	621,713	0	C	T	76.2%	146.0 / 39.7	82	A159A (GCG→GCA)	fepD	ferric enterobactin ABC transporter membrane subunit FebD
Reads supporting (aligned to +/- strand):  ref base C (8/11);  new base T (32/29);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGTTTGCATTAACGCCCCG  >  NC_000913/621468‑621959                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       acAATCAACTGTGTGCGAGCCAGGCGACTGCCCAGCGCCGGCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGTGCAGTTGCTCCGGCGAGCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGAGTGCGAATAGCCAGCGAACCGGCTTGCCAGAAACGCAACTGAGCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAg                                                                                                                                                                                                                                                   <  2:175850/250‑1 (MQ=255)      cAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGcc                                                                                                                                                                                                                                               <  1:315972/250‑1 (MQ=255)      aaaCTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCGCAAACTGTTCAGCGCGCGACTCAGGAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCa                                                                                                                                                                                                                                              <  2:254151/250‑1 (MQ=255)        aCTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCAcg                                                                                                                                                                                                                                            >  1:254390/1‑251 (MQ=255)                    ccACGCGACGGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCCCCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGtg                                                                                                                                                                                                                                <  2:97152/251‑2 (MQ=255)                    ccACGCGACGGCCCAGCGCCGTCGCGGTGGCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGtt                                                                                                                                                                                                                                <  2:354831/251‑1 (MQ=255)                                  aGCGCCGTCGCGGGTTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGTCTCAGTAATAGCGCAGTTGCGCCGGCGATCAGCACCGGGATCAGCACCTCTTTTAAGGTATGTAGAGTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGTCGGTGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCAAGCCCGCCAGGGTTAAACGCACCGGACt                                                                                                                                                                                                                  <  2:280913/251‑1 (MQ=255)                                   gcgcCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACtt                                                                                                                                                                                                                 <  1:80197/251‑1 (MQ=255)                                              gtggcGCTGCCGAGGCTCAAACTGTGCGACGGGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTGTTAAGGTATGTAGATTGCGAAGATCCAGCGAACCGGCTTGCAAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCTCCTCGCCCGCCAGGGTTAAACGCACCGGACGTAACTGCccgcc                                                                                                                                                                                                      <  2:421082/247‑1 (MQ=255)                                               tgtCGCTGCCGCGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCccgccg                                                                                                                                                                                                     >  1:98387/1‑251 (MQ=255)                                                gtCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCccgccg                                                                                                                                                                                                     <  1:221951/250‑1 (MQ=255)                                                gtCGCTACCGAGGCTCAAACTGTTCAGCGCGCGACTCAGAAATCGCGCAGTAGCTCCGGCGCTCAGCACCGGGATCAGACCCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACAGGCTTGCCAGAAACGCACCTGATCGTAGACGTCAGGATTAAGCAGGGCGCTGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGCCTTACCTGCCagccgc                                                                                                                                                                                                    >  2:385415/1‑251 (MQ=255)                                                   gCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGAGGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCt                                                                                                                                                                                                 <  1:385415/251‑1 (MQ=255)                                                      gCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCCGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGAc                                                                                                                                                                                              >  1:178730/1‑251 (MQ=255)                                                          aGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTGCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGcc                                                                                                                                                                                          <  1:187414/251‑1 (MQ=255)                                                            gCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCgg                                                                                                                                                                                        <  1:311436/251‑1 (MQ=255)                                                              tCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGAACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGt                                                                                                                                                                                       >  1:256343/1‑250 (MQ=255)                                                                 aaCTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAg                                                                                                                                                                                   <  1:144924/251‑1 (MQ=255)                                                                  aCTGTTCAGCGCGCGACTCAGAAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCCGGGTTACACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTaacgg                                                                                                                                                                                  >  2:22182/1‑248 (MQ=255)                                                                     gTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTACGGAATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGATAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGca                                                                                                                                                                                >  1:260726/1‑250 (MQ=255)                                                                            gggcgACTCAGTAAAAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCCCTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCACCTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTCCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAAcaatcaa                                                                                                                                                                        <  2:341711/249‑1 (MQ=255)                                                                                    cAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTACGGTATGTAGCTTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCCGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGcc                                                                                                                                                                 >  1:387333/1‑250 (MQ=255)                                                                                       tAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGccacc                                                                                                                                                              >  1:115509/1‑250 (MQ=255)                                                                                         aTAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAgc                                                                                                                                                           >  1:282876/1‑251 (MQ=255)                                                                                           aGCGCAGTTGCTCCGGCGACCAGCACCGGGATCCGCACCACTTTTACGGTATGTAGCTAGCGCATATCCAGCGCACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTACGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTACACGCACCGGACTTCACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAgcgc                                                                                                                                                         >  2:74201/1‑251 (MQ=255)                                                                                                    gctccgactagcAGCACCGGTATCAGCACCCCGGTTAAGGTAGGTAGATTGCGAATAGCCAGCGAACCGGTTTGCCGGAACCGCAACGGATCGAAGACGCAAGGATTAAGCGGGGCGAT‑CCGGCGGTCCGGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTGAAGCGCACCGGGCTTAACGGCCCGCCGTCCTGACTGCCGGTAAAGGCAACATTCAATGATGCCACCAGCGCCCCGGCGaa                                                                                                                                                <  2:194421/240‑1 (MQ=255)                                                                                                           cGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCTAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGGCATGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATc                                                                                                                                         <  1:172573/251‑1 (MQ=255)                                                                                                             aTCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGc                                                                                                                                       >  1:9130/1‑251 (MQ=255)                                                                                                             aTCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACGGCCCGCCGCCCTGACTGCCGGTAAAGGCAACCATCAATGAGGCCACCCGCGCCCCGGCGAAGGCCACCGc                                                                                                                                       >  2:132512/1‑251 (MQ=255)                                                                                                               cagcacctggaaacgCACCACTTATAAGGTATGTACATTGCGAATATGCAGCGAACCGGCTTGCCAAAAACGCAACTGATCGTAGACGTCAGAATGAAGCAGTGCGATGCCGCTGGTC‑AGTCCTTCCACCACCGCTGCCAGCGCCACTCCCGCCAGGGTTAAACGCACCGGTCTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGTCCATCGCCa                                                                                                                                     <  2:105753/237‑1 (MQ=255)                                                                                                               cAGCACCGGGATCAGCACCGCTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAACCGCAACTGATCGTAGACGTCAGGATGAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCTGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCa                                                                                                                                     <  1:292167/251‑1 (MQ=255)                                                                                                                aGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAg                                                                                                                                    <  2:115509/251‑1 (MQ=255)                                                                                                                 gCACGGGGATCAGCCCCACGTGTAAGGTATGGAGATTGCGATTATCCAGCGACCCGACTGGCCAGAAACGCATCTGATCGAAGACGTCAGGATTAAGCAGGGCAACTCCGCTGGTC‑ATTCCATCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTACCTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAAGCAATGAGGCCACCAGCGCTCCGGCGAAGTCCATCGCCAgt                                                                                                                                   <  2:40624/251‑1 (MQ=255)                                                                                                                   aCCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAgttg                                                                                                                                 >  1:74616/1‑251 (MQ=255)                                                                                                                    ccGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCACCAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAgttgt                                                                                                                                >  2:148427/1‑251 (MQ=255)                                                                                                                     cGGGATCAGCACCTCTTTTAAGGTATGTAGATGGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACCTCGGGATTAGGCAGGGCGATGCCGCGGGTC‑ATTCCTTCCCGCACCGCTGCCAGCGCCACGCCCGCCAGGGGTAAACGCACCGGGCGTAACTGCCCGCCGCCCTGACTGCCTGTAAAGGCAACAATCAATTAGGCCACCAGCGCCCCGGCTAAGGCGAGCGCCAgttgtt                                                                                                                               <  2:384634/251‑1 (MQ=255)                                                                                                                               accacTTTTAAGGTATGTAGATTGAGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATCCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCCCGCCCGCAAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGAGAAGGCCATCGCAAGTTGTTCCTGCGCgg                                                                                                                      >  2:417118/1‑250 (MQ=255)                                                                                                                                  acTTTTAAGGTATGTAGATTGCGCATATCCAGCGAACCCACTTGCCAGAAACGCAACTGATCGTCGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCACGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAagg                                                                                                                  >  1:153501/1‑250 (MQ=255)                                                                                                                                          ggTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACCATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGCAGAGTAACCa                                                                                                           >  2:353179/1‑250 (MQ=255)                                                                                                                                          ggTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCaa                                                                                                          <  1:147901/251‑1 (MQ=255)                                                                                                                                               gTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAgc                                                                                                     >  1:18120/1‑251 (MQ=255)                                                                                                                                                 aGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGCGGCCACCAGCGCCCCGGCGAAGGCCATCGCCCGTTGTTCCTGCGCGGCAGAGCAACCAAACAgagc                                                                                                   >  2:937/1‑248 (MQ=255)                                                                                                                                                  gATTGCGAATATCCAGCGAACCGGCTTGCCCGAAACGCAACAGATAGTAGACGCCCGGATTAACCAGGGCGATCCCGCTGGTC‑CGTCCTTCCAGCACCGCTGCCAGCCCCACGCCCGCCCGGGTTAAACCCACCGGACTTAACTGCCCGCCCCCCTGACTGCCGGGAAAGGCAACAATCAATGAGGCCACCCGCGCCACGGCGAAGCCCAACGCCAGTTGTTCCTGCGCGGAAGGGTAACCaaccactgcc                                                                                                  >  2:376978/1‑243 (MQ=255)                                                                                                                                                  gATTGCGAATATCCAGCGAACCGGCTTGCCAGAAGCGCAACTGATCGTAGACGCCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCACCGCCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCGCCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCGGTTGTTCCTGCGCGGAAGAGTAACCAAACAGcgcc                                                                                                  <  2:76357/251‑1 (MQ=255)                                                                                                                                                  gATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACGGATCGGAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCGATCGCCAGTTGTTCCTGCGCGTAAGAGTAACCAAACAGcgcc                                                                                                  <  1:164535/251‑1 (MQ=255)                                                                                                                                                   aTTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCAc                                                                                                                                                                                                                                             >  1:135837/1‑111 (MQ=255)                                                                                                                                                   aTTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCAc                                                                                                                                                                                                                                             <  2:135837/111‑1 (MQ=255)                                                                                                                                                       cGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTCACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCCATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGGGTAACCAACCAGCGCCGCAcc                                                                                             >  2:298142/1‑251 (MQ=255)                                                                                                                                                       cGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAcc                                                                                             <  1:306825/251‑1 (MQ=255)                                                                                                                                                       cGAATATACAGCGAACCGGCTTGCCAGAAACGCAACTGATCGAAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCCGTTGTTCCTGCGCGGACGAGTAACCAAACAGCGCCGCAcc                                                                                             >  2:173188/1‑251 (MQ=255)                                                                                                                                                        gAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                                                            >  1:78342/1‑251 (MQ=255)                                                                                                                                                        gAATATCCAACGAACCGGCTTGCCAGAACCGCAACGGAGCGTAGACGTCAGGATTAAGCAGGGCGATGCCTCTGGTC‑AGTCCTTCCAGCTCCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGTCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                                                            <  1:287885/251‑1 (MQ=255)                                                                                                                                                        aaaTAGCCAGCAAACGGCCTTGCCAAAAACGCACCTGCACGGACACGCCGGGATTAAGCGGGGCGATACCGCTGGCC‑AGCCCTTCCACCCCCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCCCCGGACTTAACTGCCCGCCGCCCGGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                                                            <  2:229644/250‑1 (MQ=255)                                                                                                                                                          atatCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTGCTTCCAGCACCGCTGCCAGCGCCCCGCCCGCAAGGGTTAAACGCACCGGACTTCACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAACGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAccc                                                                                            >  1:402800/1‑249 (MQ=255)                                                                                                                                                                aGCGAACAGGCTCGCCAGAAACGCAACCGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGCC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAAt                                                                                    >  2:33899/1‑251 (MQ=255)                                                                                                                                                                        ggCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGACGCTGGTC‑AGTCATTCCCGCAACACTGCCAGCGCCACGCCCGCCAGGGTTAAAAGCACCGGACTTAACTCCCCGCCGACCAGACTGCCGATAAATGCAAAAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGAAACCAAACAGCGCCGCACCCAGCACAAGGGCAAAGc                                                                            >  1:194155/1‑251 (MQ=255)                                                                                                                                                                          aacctcAGAAACGCAACTGCGCGTAGCCGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGCCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTg                                                                          <  2:66730/246‑1 (MQ=255)                                                                                                                                                                             gCCAGAAACGCAACTGAGCGTAGACGTCAGGATTAAGCAGGGGGATGCCGCTGTTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCTGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAGTAGTAATCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGc                                                                        <  1:429052/250‑1 (MQ=255)                                                                                                                                                                             gCCAGAAACGCAACGGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCCCCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCt                                                                       <  2:408640/251‑1 (MQ=255)                                                                                                                                                                              ccAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTc                                                                      >  2:194285/1‑251 (MQ=255)                                                                                                                                                                               cAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTcc                                                                     <  1:325104/251‑1 (MQ=255)                                                                                                                                                                               cAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTg                                                                          >  1:66730/1‑246 (MQ=255)                                                                                                                                                                                  aaaCGCAACTGATCGTAGACTTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCCCGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGc                                                                  <  2:9130/251‑1 (MQ=255)                                                                                                                                                                                         aCTGATCGTACACGTCAGGACTAAGCAGGGCGATGCAGCTCGCC‑AGTCATTCCAGCCCCGATCCCAGCGCCCCGCCCGACAGGGTTACCCGCACCGGACTCAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCACCAATCAATGAGGCCACAAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGCAGCGTACCCACACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTcccaccttccc                                                            >  2:367725/1‑241 (MQ=255)                                                                                                                                                                                            gATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCCCCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCTCCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCca                                                        <  1:17172/251‑1 (MQ=255)                                                                                                                                                                                            gATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGcc                                                         >  1:219437/1‑250 (MQ=255)                                                                                                                                                                                                   gactccagtaTTAAGCAGGGCGATGCCGCTGGCC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCTCCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCTCCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCTCCGCATCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCaagcaagc                                                 <  1:417118/242‑1 (MQ=255)                                                                                                                                                                                                    ccGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑CGCCCTTCCATCCCCGTGGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCTCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCTGAAGAGTATCCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTGCGGCGTTCACGCCAAGCAAGcc                                                <  2:179970/250‑1 (MQ=255)                                                                                                                                                                                                       tCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCggg                                             >  1:271553/1‑251 (MQ=255)                                                                                                                                                                                                       gcgggATTAAGCAGGGCGATGCCGCTGGTC‑ATTCCTCCCGGCACCCCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCggg                                             <  2:98387/248‑1 (MQ=255)                                                                                                                                                                                                           gATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCGCCTGCGTTCACGCCAAGCAAGCCggggacg                                         >  2:111126/1‑248 (MQ=255)                                                                                                                                                                                                               aaGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCACCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCa                                      >  1:164682/1‑250 (MQ=255)                                                                                                                                                                                                                   aGGGCGATGCCGCTGGTC‑AGGCCTTCCAGCACCGCGGCCAGCGCCTCGCCCGCCAGGGTTAAACGCACCGGCTTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCATCAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCAGCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTgg                                 <  1:74201/251‑1 (MQ=255)                                                                                                                                                                                                                           gCCGCTGGTC‑AGTCCGTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGCCTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCATACAGCTCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCGAGCCGGGGTCGGCAAGGGGGTTTCggg                         <  2:387333/251‑1 (MQ=255)                                                                                                                                                                                                                                     cgtccctcacaacACCGCCGCCGCCGCCCCGCCCGCCAGGGTTAAACGCACCGGACTTATATGCCCGCCGCCCGGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCAACCAGCACAATGGCAAAGCTGGCTCCGGCGTTCCCGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGttt                <  2:351454/239‑1 (MQ=255)                                                                                                                                                                                                                                     cagcccTTCCAGCCCCTCCCCCAGCACCCAGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCAGCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGGACAATGGCAAAGGTGGCTCCGGCGTTCTCGCCAAGCAAGCCGGGGTCGGCAGGTGAGTTTCGGGCGATGGttt                <  2:344375/247‑1 (MQ=255)                                                                                                                                                                                                                                      agcccTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTc                                          <  2:170430/221‑1 (MQ=255)                                                                                                                                                                                                                                      aGTCCTTCCAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTc                                          >  1:170430/1‑224 (MQ=255)                                                                                                                                                                                                                                      aGTCCTTCCAGCACCGCCGCCAGCCCCACGCCCGCCAGGGTTACATGCACCGGACTTAACTGCCCGCAGCCCTGACTGCCGGTCAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCACCGCCAGTTGTTCCTGCGCGGCAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCGGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGGTTCGGGGGAGGGttgg               >  2:32194/1‑249 (MQ=255)                                                                                                                                                                                                                                          cTTCCAGCACCGCTGCCAGCGCCAAGCACGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGAAACCAAACAGCGCCGCAACCAGCACAATGGCAAAGCTGGCTCAGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGATTGCAtt           >  2:213274/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccAGCACCGCTGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGTTTGCATTAAc        >  1:315248/1‑251 (MQ=255)                                                                                                                                                                                                                                                gCACCGCTGCCAGCGCAACGCCCGCCAGGGTTAAACGCAACGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAACGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGTTTGCATTAACGcc     >  1:161261/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ccgctgCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGTGTGAGGGTTTGCAGTAACGCCCCg  <  1:343503/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTC‑AGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAATCAATGAGGCCACCAGCGCCCCGGCGAAGGCCATCGCCAGTTGTTCCTGCGCGGAAGAGTAACCAAACAGCGCCGCACCCAGCACAATGGCAAAGCTGGCTCCGGCGTTCACGCCAAGCAAGCCGGGGTCGGCAAGTGGGTTTCGGGTGAGGGTTTGCATTAACGCCCCG  >  NC_000913/621468‑621959

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