Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,707	A→G	G166G (GGT→GGC)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,707	0	A	G	99.0%	315.1 / NA	97	G166G (GGT→GGC)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (43/53);  total (43/54)
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.74e-01

AGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC  >  NC_000913/4193470‑4193950                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    aGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGccc                                                                                                                                                                                                                                        >  2:483063/1‑250 (MQ=255)  gACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCTGCTTGAGTAATGTCGCGGACGACGGCGATACTGCCGACTCCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCAc                                                                                                                                                                                                                                       <  1:34850/251‑1 (MQ=255)   aCCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTAGCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACg                                                                                                                                                                                                                                      <  1:518631/251‑1 (MQ=255)     cGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATAATGCCGACAACCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACggg                                                                                                                                                                                                                                    >  1:433756/1‑250 (MQ=255)     cGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGt                                                                                                                                                                                                                                    >  1:423043/1‑251 (MQ=255)       aTAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTgg                                                                                                                                                                                                                                  <  2:4910/251‑1 (MQ=255)       aTAATCAGCACCTGGCTGTCGAGCAGTTTTTGCCGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTgg                                                                                                                                                                                                                                  <  1:379743/251‑1 (MQ=255)          aTCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGAt                                                                                                                                                                                                                               <  1:155925/251‑1 (MQ=255)                 ccTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATGTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGc                                                                                                                                                                                                                        <  1:96063/251‑1 (MQ=255)                  cTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGcc                                                                                                                                                                                                                       >  1:14926/1‑251 (MQ=255)                   tGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGTAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCAAATCGCCACGGGGGGATAATCCGcc                                                                                                                                                                                                                       >  1:47389/1‑250 (MQ=255)                   tGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGAGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGccc                                                                                                                                                                                                                      >  1:461933/1‑250 (MQ=255)                    ggCTGTCTAGCTGTTTTTGCAACCCGTCGAGAGCGCTATCGTCCAGCAGGATTTGGCGGCTCTAACGCATAACGTCC‑CGGTCATTCATCACCAACTCATCCAATTTCCCGCCACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGTGTAATGGCGCTGACTACGGCGATACTGCCGACACCCGT‑TGCTCACACCGCAGGCGCGCGGGCCAGAATGATGTCGCCACTCGCCACGGTGGGGTAATCCGcc                                                                                                                                                                                                                       >  2:285770/1‑249 (MQ=255)                    ggCTGTCGAGCAGTTTTGGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGcc                                                                                                                                                                                                                       <  2:546694/249‑1 (MQ=255)                    ggCTGTCGAGCAGTGTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAg                                                                                                                                                                                                                     <  1:284122/251‑1 (MQ=255)                         gcGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCt                                                                                                                                                                                                                <  1:210133/250‑1 (MQ=255)                          cGAGCAGTTTTTGCTGCCCGTCCAGCGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCCGACTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAACCGCCACGGTGGGATAATCCGCCAGTCGcac                                                                                                                                                                                                               >  1:296657/1‑249 (MQ=255)                                gTTTTTGCTGCCCGTCCAGCGAGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATg                                                                                                                                                                                                         >  1:431125/1‑251 (MQ=255)                                gTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCCGGCGCGCGCGCCAGACTGATGCCGCCAATCGCCACGGCGGGATAATCCGCCAGTCGCTCAACATg                                                                                                                                                                                                         >  1:105202/1‑251 (MQ=255)                                       ctgcCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGcc                                                                                                                                                                                                  >  1:140949/1‑251 (MQ=255)                                         gcCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCa                                                                                                                                                                                                 <  1:529443/250‑1 (MQ=255)                                             gTCCAGAGCGATATCGTCGGTCAGGTTTTGGCGGCTAGAACGCATAAAGTCA‑CGGTCTTTCAACGCCAACTCCTGCAATTTCCAGCAACTGGGCCGTTGCCAAACGCCAGTCTGCGGCTTTAGTAATGGCGCTGACGACGGCGATACTGCCGACCCCCGT‑TGCTATCACCGCAGGCCCGCGTGCCAGACTGATGCCTCCAAGCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGt                                                                                                                                                                                            <  2:390759/251‑1 (MQ=255)                                                     cGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGccc                                                                                                                                                                                    <  1:133273/251‑1 (MQ=255)                                                         atCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGGTCCAGCCCCTGc                                                                                                                                                                                >  2:368248/1‑251 (MQ=255)                                                           cgtcgAGCAGGATTTGGCGGCTATAACTCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCGGACGACGGCGATTCTGCCGACACCCGT‑GGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCgg                                                                                                                                                                              <  1:569015/251‑1 (MQ=255)                                                             ccgAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTg                                                                                                                                                                            <  1:79734/250‑1 (MQ=255)                                                                           gcggcTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgt                                                                                                                                                              <  1:84571/251‑1 (MQ=255)                                                                             ggcTATAACGCATAAAGGCAGTGGT‑ATCCATCGCCAACTCCTGCAATTGCCAGCAGCGGGGCCGTGGCCAAACGCCAGTCTGCGGCTTGAGTAAGGGCGCTGACGACGGGGATATTGCCGACACCCGT‑GGCTATCACGGCAGGCGCGCGTGCAATACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGATCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgttt                                                                                                                                                            <  2:448857/251‑1 (MQ=255)                                                                             ggcTATAACGCAAAAAGTCA‑CGGTCATTTATCGCCACCTCCTGCAATTACCAGCAACTGTGCCGTTGCTAAAGGCCGGTCTGCGGCTGGAGTAATGGCGCTGACGACGGCGATACTGCCGAGACCCGT‑TGCTAGCAGCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCATCATGCCg                                                                                                                                                                                                      <  1:302021/209‑1 (MQ=255)                                                                              gcTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCCGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGACGGCATCTgtttg                                                                                                                                                           >  2:360057/1‑251 (MQ=255)                                                                                 ataACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAAGCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttggtt                                                                                                                                                        <  1:487133/251‑1 (MQ=255)                                                                                     cGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCg                                                                                                                                                    <  1:160309/251‑1 (MQ=255)                                                                                     cGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGAAGCCGGCGATACTGACGCAACCCGT‑TGCCATCACCGAAGACGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGATGGCATCTgtttggattg                                                                                                                                                      >  1:526017/1‑249 (MQ=255)                                                                                          aaaGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAGCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCggg                                                                                                                                               <  2:443153/251‑1 (MQ=255)                                                                                          aaaGTCA‑CGGTCATTCATCGCCAACTCCTGCAATGTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGTCTTGAGTAATGGCGCTTACGACGGCGATACTGCCGACACCCGT‑TGCTAGCACCGCAGGCGCGCGTGCGAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCGGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCggg                                                                                                                                               <  1:364341/251‑1 (MQ=255)                                                                                               cA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAacac                                                                                                                                          >  1:461339/1‑251 (MQ=255)                                                                                                 acggtacttaatCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGAGTGCCAGACTGAGGCCGCCAATCGCCACGGAGGGGTAATCCGCCAGTCGCTCAACAGGCCGTGCCAGCTGTTCCAGCCCCTTCGGGGCAGAAGGCCTCTGTTTGGTTTGCGTCGGGAacac                                                                                                                                          >  2:433054/11‑250 (MQ=255)                                                                                                   ggTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATg                                                                                                                                       >  1:147336/1‑251 (MQ=255)                                                                                                   ggTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATg                                                                                                                                       <  1:25865/251‑1 (MQ=255)                                                                                                      cattcatCGCCAACCCCTGCATTTTCCAGCATCTGTGCCGTTGCTAAACGCCAGTCTGCGGCTGAAGTAATGGCGCTGACGACGGCGATACTGCCGCCTCCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTct                                                                                                                                    <  2:451066/251‑2 (MQ=255)                                                                                                       attcatCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTctc                                                                                                                                   <  2:150361/251‑3 (MQ=255)                                                                                                         tcatCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACCCCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTctca                                                                                                                                  <  2:432069/250‑4 (MQ=255)                                                                                                            ccGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTTCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGGTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAgcgc                                                                                                                              <  2:217485/250‑1 (MQ=255)                                                                                                              gccaaccacTGCAATTTCCAGCAACTGAGCCGTTGCCAAACGCAAGTCCGCGGCTTGAGTAATGGCGCTGAAGACGGCGATACTGCCGACACCCGT‑TGCTATGACCGCAGGCGCGCGTGACAGACTGATGCCGCCAATCGCCACGGGGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCAGTTCCAGCCCCTGCGGTGAAGACGGCCTCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGcgcg                                                                                                                             >  2:53092/9‑250 (MQ=255)                                                                                                                 aaCTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGatat                                                                                                                         >  1:576890/1‑251 (MQ=255)                                                                                                                     ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCCATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGCACACATGTCTCAGCGCGAtata                                                                                                                        >  1:100145/1‑248 (MQ=255)                                                                                                                     ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAag                                                                                                                      <  1:472114/250‑1 (MQ=255)                                                                                                                        gCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgg                                                                                                                   <  1:291805/250‑1 (MQ=255)                                                                                                                         cAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAgggg                                                                                                                 <  1:243068/251‑1 (MQ=255)                                                                                                                          aaTTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGc                                                                                                                >  1:473500/1‑251 (MQ=255)                                                                                                                           aTTTCCAGCAACTGTGCCGTTGCCAATCGCAAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCg                                                                                                               <  2:505905/251‑1 (MQ=255)                                                                                                                                cAGCAACTGTGCCGTTGCCCAACCCCAGTCTGCGGCTTGAGTCATGGCGCTGACGCCGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCCCGGTGGGATCATCCGCCCGTCGCTCAACATGCCGTGCCCGCTGTTCCAGCCCCTGCGGTGCAGAATGCATCTGTTTGGTTTGCGTCGGCAACACCTGTCTCAgcgctatctaaccgggtctttcct                                                                                                          >  2:135261/1‑231 (MQ=255)                                                                                                                                cAGCAACTGTGCCGTTGCCAAACGCCAGTCGGCGGCTTGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgctg                                                                                                          <  2:338682/251‑1 (MQ=255)                                                                                                                                   cAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGccc                                                                                                       >  2:129082/1‑250 (MQ=255)                                                                                                                                     aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAgc                                                                                                     <  1:453009/251‑1 (MQ=255)                                                                                                                                        gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGGGTAATGGCGCTGACGACGGCGCTCCGGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTGCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCTCGATATAAGAGGGGCGTGCTGCCAgcgc                                                                                                   <  2:589360/250‑1 (MQ=255)                                                                                                                                         tgCCGTTGCCAAACGCCAGTCTGCGGGTTGCGTATTGGCACTGACGCCGGCGATACTGCAGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTGCGCCAGTCGCTCAACATGCTGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCATCGCGCTATACGAGGGGCGTGCTGCCAGcgcgc                                                                                                 >  2:64253/1‑250 (MQ=255)                                                                                                                                                gCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAt                                                                                          <  2:16044/251‑1 (MQ=255)                                                                                                                                                 ccaagctccAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                         <  1:190041/244‑1 (MQ=255)                                                                                                                                                 ccAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAtt                                                                                         <  1:460487/251‑1 (MQ=255)                                                                                                                                                 ccAAACGACACTCTGCGGCTTGAGTAATGGCGCTGACGACGCCGATACTGCCGCAACCAGT‑TGCTATCACCGCAGGCGAGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATACTCCGCCCGCCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGCAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCCGCGCGGTATAAGAGGGGCGTGCTGCCCGCGCGACGCCGAtt                                                                                         >  2:378499/1‑251 (MQ=255)                                                                                                                                                  caaccGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                        <  2:481917/247‑1 (MQ=255)                                                                                                                                                  cAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTACCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                        >  2:263377/1‑251 (MQ=255)                                                                                                                                                   aaaCGCCAGGCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                       >  1:425398/1‑251 (MQ=255)                                                                                                                                                      cGCCAGTCAGCGGGTTGAGTAAGGGCGCTGACGACGGCTATTCTGCCGACACCCGT‑TGCTAACACCGCAGGCGCGTGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGAGAATCCGCCAGGCGCTAAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCAACTGTTAGGTTTGCGTCGGGAACACATGTCTCAGTGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCat                                                                                    <  1:44708/251‑1 (MQ=255)                                                                                                                                                       gCCAGTTTGCGGCTTAAGTAATGGTGCTGACGTTGGCTATACGGTCGCCACCCGGGTG‑TATCACCGCCGGCGCCGGTTCCAGACTGATGCCGCCAATCTCCACGGTTGGATAATCCTCCAGTCGCTCAACATGCCGCGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGTATCTGTGTGGTTTGCGTCTGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGGCGAtgaccata                                                                                   <  2:541795/251‑8 (MQ=255)                                                                                                                                                           gTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGCCACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGt                                                                               >  1:210875/1‑251 (MQ=255)                                                                                                                                                            tCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAAAGGGGCGAGCTGCCAGCGCGACGTCGATTTCCCTATCGTc                                                                              >  1:401953/1‑251 (MQ=255)                                                                                                                                                             cTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTc                                                                              >  2:569788/1‑250 (MQ=255)                                                                                                                                                              ggCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAt                                                                            <  1:401030/250‑1 (MQ=255)                                                                                                                                                               gCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAtg                                                                           <  1:368248/251‑1 (MQ=255)                                                                                                                                                                cGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAtgt                                                                          >  1:213780/1‑251 (MQ=255)                                                                                                                                                                    ttGAGTAAGGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTc                                                                       <  2:440716/250‑1 (MQ=255)                                                                                                                                                                                gctgaagccgGCGATACTGCCGACACCCAT‑TGCTATCACCGCAGGCGCGCGAGCCAGCCTGCTGCCGCCAATCGCCACGGGGGGACAATACGCCGGTCGCCCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGCCGTCTAGTCTCATATCGTCATATGTCGAAACGCCCCGcc                                                          >  2:570979/9‑251 (MQ=255)                                                                                                                                                                                        cgGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGGCATGTGTCGAAACGCCCAGCCGCAGGCCt                                                  <  1:483063/251‑1 (MQ=255)                                                                                                                                                                                         ggCGATACTGCCGTCACCCGT‑GGCGATCACCGCGGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGGGTCGAAACGCCCAGCCGCATGCCTg                                                 <  2:583866/251‑1 (MQ=255)                                                                                                                                                                                           cGATCCTGCCGACCCCCGT‑TGCGATCCCCGCAGGCGCGCGTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGcc                                               <  2:165818/251‑1 (MQ=255)                                                                                                                                                                                            aaTACTGCCGACACCCGT‑TGATATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGGCGCTCAACATGCCGTGCCAGATGTTCCAGCCCCTGCGGTGCAGGAGCCATCGGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGTGGCGTGCTGCCAGCGCGACGTCGCTTTCCCTATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGccc                                              >  2:481798/2‑250 (MQ=255)                                                                                                                                                                                                cTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGAGCCAGCTGTTCCAGACCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCcgcgcg                                          >  1:212143/1‑251 (MQ=255)                                                                                                                                                                                                     gACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAACAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATCGTCATGTGTCGAACCGCCCAGCCGCAGGCATGCCCCGCGGCTgg                                     >  2:486090/1‑251 (MQ=255)                                                                                                                                                                                                          ccGT‑TGCGATCACCGCAGGCGGGCGTGCCAGACTGATGCCGACAAGCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTg                                <  1:331313/251‑1 (MQ=255)                                                                                                                                                                                                            gT‑TGCGATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgag                              <  2:461339/251‑1 (MQ=255)                                                                                                                                                                                                                gCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                           <  1:430389/251‑1 (MQ=255)                                                                                                                                                                                                                gCTATCACCGCAGGCGCTCTTGCCAGACGGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCCTAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGGCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCTTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                           <  2:187380/251‑1 (MQ=255)                                                                                                                                                                                                                gCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATc                           <  1:243289/251‑1 (MQ=255)                                                                                                                                                                                                                    tCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTg                       >  2:137907/1‑251 (MQ=255)                                                                                                                                                                                                                      cccGCAGGCGCGCGTGCCAGACTGTTGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                     <  1:559215/250‑1 (MQ=255)                                                                                                                                                                                                                      aCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAACCTCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGc                     >  1:71286/1‑251 (MQ=255)                                                                                                                                                                                                                         gCAGGCGCGCGTGCCAGACTGATGCAGCCAATCCCCACGGTGGGATAATCCACCAGTCGCTCAACATGCCGTGCCAGCTGTTCCCGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGACCACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGCCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGCGATCGGTGGCTTg                  >  2:574211/1‑251 (MQ=255)                                                                                                                                                                                                                          cAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                 >  1:556188/1‑251 (MQ=255)                                                                                                                                                                                                                           aGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGCTAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGc                 >  1:451203/1‑250 (MQ=255)                                                                                                                                                                                                                               gctcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc            >  1:461831/4‑251 (MQ=255)                                                                                                                                                                                                                               gcgcgTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATc            <  2:140949/251‑1 (MQ=255)                                                                                                                                                                                                                               gcgcgTGCCAGACTGATGCCGACAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGCAATGGGTCGAAACGCCCAGCCGCAGGCCTGCCGAGCGGGTGGCATTGAGATCGGTGGCTTGCaactc            >  2:193453/1‑248 (MQ=255)                                                                                                                                                                                                                                     gCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTg      <  1:111914/251‑1 (MQ=255)                                                                                                                                                                                                                                       cAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGcc    <  1:137907/251‑1 (MQ=255)                                                                                                                                                                                                                                         gACTGATGCCGACAATCGCCGCGGGGGGATAATCCGCCAGGCGCTCAAAATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGAAATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGACCAGCCGCAGGCCTGCAGAGCGGATGGCATTGCGATCGGTGGCTTGCAAATCTTCCTGcccc  >  2:46950/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    AGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCA‑CGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGT‑TGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC  >  NC_000913/4193470‑4193950

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