Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,924	C→T	Q605* (CAG→TAG)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,924	0	C	T	94.6%	199.6 / 0.7	76	Q605* (CAG→TAG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (27/43);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

CAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGAC  >  NC_000913/4445686‑4446165                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  cAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGAACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCgggg                                                                                                                                                                                                                                       >  1:146490/1‑250 (MQ=255)  aGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGa                                                                                                                                                                                                                                       >  1:459834/1‑250 (MQ=255)     gCCGAAGGGAAAGGCGAGCTGGGGGTAAAAGAGCGCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGAGCGCCGGGAAACt                                                                                                                                                                                                                                   <  2:28798/251‑1 (MQ=255)      cgtagtttaatggCGAGCTGGGGATAAAAGAATTCAATCTTGAGCCCTCCAGCAACGCGCCGAGGCTGAATACCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCCTGTGAA‑AGTACGCTGCACGGTGTAAGCGACACAACTACTGGCAGTTAGCCCCGCGCGCGTCCTTCGCTTGCAGGAACTTTCCTT‑GGCGCAGGTTCGCGTCGAAGGCGATATCAA‑ATCCTCCGATTAGATCGCCGGTAAAcgc                                                                                                                                                                                                                                  <  2:294913/240‑3 (MQ=255)       cGAAGGGAAAGGCTATCGGGGGGTAAAAGATCTCAATCTGGTCGCCACCACCCACGCCCCGGGGCTGGATATCGCGCTGCCGGGGCTTGGCGGGACAGCGAAAGGGCTTGTTAA‑AGGACGCGGCACGGTGGATGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCCCAGGCTCGCGTGGAAGGCGCTATCAA‑ATCCTCCGATTAGATCGCCGGGAAAcgcg                                                                                                                                                                                                                                 <  2:323675/251‑4 (MQ=255)         aaGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGAGCAGGGTCGCGGGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAACCTCgcc                                                                                                                                                                                                                               >  2:361122/1‑249 (MQ=255)          aGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGGTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACg                                                                                                                                                                                                                              <  2:215878/251‑1 (MQ=255)            tGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGt                                                                                                                                                                                                                             >  2:154358/1‑250 (MQ=255)             gAAAGGGGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGTTGGAAGCGCCACAACTACCGGCCGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACGTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTAc                                                                                                                                                                                                                           <  2:75888/251‑1 (MQ=255)                  gCGAGCGGTGGGTAAAAGATCGCAAGCGTGATGCCACCATCAACGCGCCGGGGCTGGAGAACGCGCTGCCGGGGCGTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCATGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATAAGACCTCCGGGAAACTCGACGTACGCGtt                                                                                                                                                                                                                      <  2:292584/251‑1 (MQ=255)                       cTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGGGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGGAACTCGACGTACGCGTTGAGc                                                                                                                                                                                                                  >  2:195047/1‑250 (MQ=255)                             gTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGAGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGAGATCAA‑ATCCACCGATCAGATCGCCGGGGAACTCGACGTACGCGTTGAGCAAATTTc                                                                                                                                                                                                           >  2:38230/1‑251 (MQ=255)                                        caatcttgagtcccccactaccgcgTCGGGGCTGGATAACGCGCGGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGG‑GAACAGTACCCGGCACGGTGGAAGCTCCACAACTACTGGCTGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTGCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCTGGGAAACTCGACTTACGCGTTGAGCAAATTTCGCAGCCGGATg                                                                                                                                                                                                <  2:478238/230‑1 (MQ=255)                                         aaTCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGt                                                                                                                                                                                               <  1:185630/251‑1 (MQ=255)                                         aaTCTTGATGCCACCATCAACGCGCCAGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGt                                                                                                                                                                                               <  1:56039/251‑1 (MQ=255)                                          aTCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGt                                                                                                                                                                                               >  1:62642/1‑250 (MQ=255)                                              tgaggCCACCATCAACGCGGCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTa                                                                                                                                                                                              <  2:428883/243‑1 (MQ=255)                                              tGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCCGGAACTTTCCGT‑GGCGCCGGTTCGCGTGGAAGGCGATATCAA‑AACCCCCGATTAGATCGCCGGGAAACTCGCCGTACGCGTTTCGCAAATTTAGACGCCGGGTGTAAat                                                                                                                                                                                           >  2:24562/1‑250 (MQ=255)                                                aTGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATc                                                                                                                                                                                        <  1:228139/251‑1 (MQ=255)                                                     accaTCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAACGGGCTGGTGAA‑AGTACGCGGCACGGTGGCAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCCGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCCA‑ATACACCGATCCGATCTCCGGGTAACTCGACGTACGCGTTGCGCAAATTTct                                                                                                                                                                                                          >  2:145819/1‑227 (MQ=255)                                                     accaTCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCg                                                                                                                                                                                                          <  1:145819/228‑1 (MQ=255)                                                        aTCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGt                                                                                                                                                                                <  1:400073/251‑1 (MQ=255)                                                          cAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTc                                                                                                                                                                               <  1:450150/250‑1 (MQ=255)                                                           aaCGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAc                                                                                                                                                                             >  1:71965/1‑251 (MQ=255)                                                             cgcgCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGAGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAccc                                                                                                                                                                           <  1:38230/251‑1 (MQ=255)                                                              gcgcCGGGGCTGGATAACGCGCGGGCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCt                                                                                                                                                                          <  1:3714/251‑1 (MQ=255)                                                                 ccGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCATCACCCTGa                                                                                                                                                                        >  1:419720/1‑250 (MQ=255)                                                                        tGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCCACGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGGTCGCGTGGAAGGCGGTATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGAACGCGTTGAGCAAATTTCGCAGCCGGAAGAAAATATCAACCTCGTCACCCTGAATGCCac                                                                                                                                                                 >  2:84789/1‑249 (MQ=255)                                                                         ggATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGGTATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGCAGGCCaa                                                                                                                                                                 >  2:483852/1‑249 (MQ=255)                                                                         ggATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaa                                                                                                                                                                 <  1:297412/249‑1 (MQ=255)                                                                             aaCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGTGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATAGCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCCAAGGCAg                                                                                                                                                           <  2:416805/251‑1 (MQ=255)                                                                                  gcTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGaaa                                                                                                                                                      <  2:108726/251‑1 (MQ=255)                                                                                  gcTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGaaa                                                                                                                                                      <  1:195047/251‑1 (MQ=255)                                                                                      ccGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagca                                                                                                                                                  <  1:251874/251‑1 (MQ=255)                                                                                      ccGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagc                                                                                                                                                   >  1:297767/1‑250 (MQ=255)                                                                                        ggggCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGAGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGAGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGAAGTACGGGTTGAGCAAATTTCGCAGCCGGATGTAAATATCCACCTCGTTACCCGAAATGCCAAAGGCAGCGaaaaac                                                                                                                                                   >  2:283666/1‑246 (MQ=255)                                                                                              tGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAg                                                                                                                                           <  1:184166/250‑1 (MQ=255)                                                                                               ggCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTTGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCTTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCt                                                                                                                                         <  1:200270/251‑1 (MQ=255)                                                                                                     aCAGGGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAATCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGtt                                                                                                                                   <  2:459834/251‑1 (MQ=255)                                                                                                             aGGGCTGGTGAA‑AGTACGCGGCTCGGTGGAAGCGCCACAACTACTGGCAGATATCCCCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCTCCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTc                                                                                                                           <  1:505097/251‑1 (MQ=255)                                                                                                                      gAA‑AGTAAGCGACACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGAAACTTTCCGT‑GGCGCAGGTTCGCGTGGCAGGCGATATCAA‑ATCCACCGCTTAGATCGCCGGGTAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGACGCACGAGCTCCATTTGCGGATTCAGGGGGAGc                                                                                                                  >  1:434188/1‑251 (MQ=255)                                                                                                                        a‑aGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCt                                                                                                                >  1:460537/1‑251 (MQ=255)                                                                                                                               gcgGCAGGGTGTAAGCGCCGCAACTACTGGCAGCTATCACCGCGCGCGTCCTGCGCGGGCGGGAACTTCCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCTCCGATTAGATCTCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCTAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAgggggcgaccgtctcc                                                                                                          <  2:438156/251‑12 (MQ=255)                                                                                                                                cgGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCg                                                                                                         <  1:154358/251‑1 (MQ=255)                                                                                                                                 ggCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAatat                                                                                                                                                                                         >  2:166923/1‑170 (MQ=255)                                                                                                                                 ggCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAatat                                                                                                                                                                                         <  1:166923/170‑1 (MQ=255)                                                                                                                                        tGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGACGGCGATATCAA‑ATCCACCGATTAAATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTAACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGCTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCtt                                                                                                 >  1:475845/1‑251 (MQ=255)                                                                                                                                               gcCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCTGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTgg                                                                                          <  2:484189/251‑1 (MQ=255)                                                                                                                                                 cacaACTACTGGCAGATATCACCGCGCGCTGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTGCGCGTGGAAGGCGATAGCAA‑ATCCACCGATTAGATCGCCGGGTAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGTCAGCGAAAAGCATCTCTATCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCa                                                                                        <  1:329689/251‑1 (MQ=255)                                                                                                                                                 cacaACTACTGGCAGATATCACCGCGCGCGGCCTGCGATGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCTTGGAAGGCGATATCAA‑ATCCACCGATTAGAGCGCCGGGAAACTCGACGTACGCGTGGAGCAAATTTCGCAGACGGATGTAAATATCAACCTCGTCTCCCTGAATGCCAAAGGCAGCGAATAGCAGCACGAGCTACTGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCa                                                                                        <  2:112367/251‑1 (MQ=255)                                                                                                                                                  acaACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGc                                                                                         >  2:295415/1‑249 (MQ=255)                                                                                                                                                   caACTACTAGCAGACATCAACGAGCGCGGCCTGCGCTGGCAGGCACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGCTATCAA‑ATCCCCCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGCATGCCAAAGGCAGCGAAAAGCAGCACGCGCTACAGTTGCGGATTCCGGGGGAGCCTGTCTCCGGGCCGCTTACTCTGGCAgg                                                                                      >  2:506289/1‑251 (MQ=255)                                                                                                                                                        actgtgaGATATCACCGCGCGCGGACTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATAGCAA‑AGCCACCGAGTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCATCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGGCTGTTTCCGGGCAGCTTAATCTGGCAGGAAGt                                                                                  <  2:41370/244‑1 (MQ=255)                                                                                                                                                           ggCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGGCTCCGGGCAGCTTAATCTGGCAGGAAGtttt                                                                               <  1:114283/250‑1 (MQ=255)                                                                                                                                                           gccAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCTGGCAGCTTAATCTGGCAGGAAGtttt                                                                               <  1:36226/248‑1 (MQ=255)                                                                                                                                                            gCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTg                                                                              >  1:15067/1‑250 (MQ=255)                                                                                                                                                             ccgATACCACCGCGCGCGGTCTGCGGTGGCAGGAACTT‑CTGTTGGCTCAGGTGCGCGTGGAAGGCGATA‑CAACATCCACCTAATAGAGCGCCGGGAAACGCGACGTACGCGTTGAGCAAATTTCGCTGCCTGAGGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCGCGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTggaatgaattgttgat                                                                            <  2:36636/249‑15 (MQ=255)                                                                                                                                                               gATAGCACCGCGCGCGGCCTGCGCTGGCAGGAACGTTCCGG‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCg                                                                          <  2:463252/251‑1 (MQ=255)                                                                                                                                                                atatCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGCTCGCGTGCCAGGCGATATCAA‑ATCCACCGATTACATCGCCGGGAAACTCGACGTACGCGTTGAGCAAGTTTCGCAGCCGGATGTAAATATCAACCTCGTCCCCATGAATGCCAAAGGCAGCGACAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCGGCTTAATCTGGCAGGAGGTTTTGATccc                                                                         >  2:162770/1‑249 (MQ=255)                                                                                                                                                                 tatCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGATTAATCTGGCCGGAAGTTTTGATCGc                                                                         >  1:452017/1‑250 (MQ=255)                                                                                                                                                                  atCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGCAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGAACGAGCTACAGTTGCGGCTTCAGGGGGAGCCAGTCTCAGGGCAGCTTAATATGGCAGGAAGTTTTGATCGCa                                                                        >  2:466506/1‑250 (MQ=255)                                                                                                                                                                          gcgcgGCCGTCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAAc                                                               <  2:155147/251‑1 (MQ=255)                                                                                                                                                                           cgcgGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACg                                                              >  1:417808/1‑251 (MQ=255)                                                                                                                                                                            gcgGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGc                                                             <  1:361122/251‑1 (MQ=255)                                                                                                                                                                                        gggaGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCGGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAAct                                                 <  2:389840/248‑1 (MQ=255)                                                                                                                                                                                                  ttCCGT‑GGAGCAGGTTCGTGTGGAAGGCGATATCAA‑ATCCTCCGATTAGATCGCCGGGAAACTCGACGGACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCGCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCGGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCaa                                         <  2:373053/249‑1 (MQ=255)                                                                                                                                                                                                       ggggCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCATCGATTAGATCGCCGGGAAACGCGACGTGCGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCg                                   <  2:454859/249‑1 (MQ=255)                                                                                                                                                                                                         ggCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCt                                 >  1:281363/1‑251 (MQ=255)                                                                                                                                                                                                         ggCGCAGGTGCGCGTGGAAGGCGATATCAA‑ATCCTCCGATTAGATCGCCGGGAAACTCGACGGTCGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTATCAATACCCGCt                                 <  2:107367/251‑1 (MQ=255)                                                                                                                                                                                                          gcgcAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCtt                                <  1:461014/251‑1 (MQ=255)                                                                                                                                                                                                                gTTCGCGTGGAAGGAGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCGGGAAGGGAACTCTTAGCAATACCCGCTTCCAGAc                          <  1:135163/251‑1 (MQ=255)                                                                                                                                                                                                                       tGGAAGGCGATATCAA‑ATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAGAGAAGAACGCTGGACAGGAACTCTTAGCACGACCCGCTTCAAGCCGCGGGtt                   >  2:135007/1‑251 (MQ=255)                                                                                                                                                                                                                                  atCAA‑ATCCACAGATAAGATCGCCGGGAAACTCGACGTACGCGATGAGCAAATTTCGCAGCCGGATGTAAATATCAACCACGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCCGCTTAATCTGGCAGGAAGTTTTGATCGCAACGACGAACGCTGGAAGGGAACTATTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTc        >  1:446785/1‑251 (MQ=255)                                                                                                                                                                                                                                        acccACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTAGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGa   <  2:15067/249‑1 (MQ=255)                                                                                                                                                                                                                                        aTCCACCGATTAGATCGCCGGGAAACTAGACGCACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGCTTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGGGGTCGcggc   >  2:495240/1‑248 (MQ=255)                                                                                                                                                                                                                                          ccACCGATTTGAACGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAAGATCAACCTCGTCCCCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCGACAGTTGCGGATTCAGGGGGAGCCTGTCGCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGAc  <  2:379626/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  CAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAA‑AGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGT‑GGCGCAGGTTCGCGTGGAAGGCGATATCAA‑ATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGAC  >  NC_000913/4445686‑4446165

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: ‑