Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,111,727	2 bp→TT	coding (865‑866/2544 nt)	opgH →	osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,111,727	0	C	T	95.6%	279.0 / 2.8	92	P289M (CCG→TTG)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (49/38);  total (53/39)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01
*	NC_000913	1,111,728	0	C	T	95.6%	282.4 / 0.6	92	P289M (CCG→TTG)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (50/37);  total (53/39)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ATTCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCG  >  NC_000913/1111486‑1111970                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                    agtCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGGTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATcg                                                                                                                                                                                                                                             <  2:63354/249‑2 (MQ=255)    cTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCg                                                                                                                                                                                                                                          >  1:685358/1‑251 (MQ=255)      ggACGGCGTTATTGGGCTTCCTGCAACTTCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTc                                                                                                                                                                                                                                        <  2:555991/251‑1 (MQ=255)      ggACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTTGATATCTGCGTc                                                                                                                                                                                                                                        <  1:503347/251‑1 (MQ=255)        aCGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGCGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGc                                                                                                                                                                                                                                      >  1:554632/1‑251 (MQ=255)         cGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAACGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTCTAACTTGGATATCTGCGTCGc                                                                                                                                                                                                                                      >  2:688252/1‑250 (MQ=255)         cGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCa                                                                                                                                                                                                                                     <  2:582149/251‑1 (MQ=255)          ggCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGc                                                                                                                                                                                                                                      >  1:530519/1‑249 (MQ=255)                   gggCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCaaaa                                                                                                                                                                                                                             >  1:171018/1‑249 (MQ=255)                           tGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCtggatg                                                                                                                                                                                                                   >  1:411184/1‑251 (MQ=255)                                 tGCTGATTGGTCGCGATAAATACAGTATATCTTCGTAAACAGTTGGCGGTGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAAAGAAGCCGTGAACCGTGTTTTTGCTGGACTGCGTGCAACGTGGGAAGCAGAAAAAGCCACAAGGAACGCCAAACACTTTGATGTCTACCTTCTTAGTGACAGTTCTAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCtt                                                                                                                                                                                                             >  2:439461/1‑251 (MQ=255)                                           tCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGa                                                                                                                                                                                                     >  1:453586/1‑249 (MQ=255)                                            cgcgATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAg                                                                                                                                                                                                                           <  2:638992/226‑1 (MQ=255)                                            cgcgATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAg                                                                                                                                                                                                                           >  1:638992/1‑226 (MQ=255)                                            cgcgATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGt                                                                                                                                                                                                  <  1:183291/251‑1 (MQ=255)                                                 tAAATACAGTATATCTGCGTCAACAGTTGGCGATGATCCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTTGGAATCAGTAAAAGCCAACGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGACGTCGGTg                                                                                                                                                                                             >  2:329795/1‑251 (MQ=255)                                                  aaaTACAGTATATCTGCGTCAACAGTTTGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATTCCTATCTGTAACGAAGACGTGAACCGTGGTTTTGCTGGCCTGCGCGCAACGGTGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTGTGATGTCTACACTCTTAGGGACAGTTATAACTTGGATAGCTGCGTCGCAGAGCAAAAATCCTCGATGTAGCTTAGCGCTTAAGTCGGTg                                                                                                                                                                                             <  2:304543/250‑1 (MQ=255)                                                  aaaTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGCTGGAGCTTATCGCTGAAGTCGGTgg                                                                                                                                                                                            >  2:141192/1‑251 (MQ=255)                                                  aaaTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGGTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTgg                                                                                                                                                                                            <  1:321867/251‑1 (MQ=255)                                                    aTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTATGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGAGGGAGCTTATCGCTGAGGTCGGTGGCg                                                                                                                                                                                          >  2:582887/1‑251 (MQ=255)                                                      aCAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGTAGCCGGTGGCgc                                                                                                                                                                                         >  2:357903/1‑249 (MQ=255)                                                          tataTCTGCGTCAACAGTTGGCGATGAACCATTAAACCAGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGACACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGACTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAgga                                                                                                                                                                                     >  1:12669/1‑249 (MQ=255)                                                                   gTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCt                                                                                                                                                                           >  1:558487/1‑251 (MQ=255)                                                                        cAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATcgc                                                                                                                                                                      <  1:439461/251‑1 (MQ=255)                                                                         cgggggctaTGATCCATAAAAACCGGACCATCGCACGGCGTGGCAT‑ATGCCTATCTGTAACGAAGACGGGAACCGGGTTTTTGCTGGACTGCGTGCAACGGTGGAAGCAGTAAATGCCACCGGGAATTCGAAACACTTTGATGTCTACATTCTGAGTGACAGTTATAACTGGGAGATCTGCGGCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGGGGCGAAGAACAGATTTTCTATcgcc                                                                                                                                                                     <  2:675261/243‑1 (MQ=255)                                                                                 aTGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAACCACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTTTCGCCGCCGCCGt                                                                                                                                                             >  2:69250/1‑251 (MQ=255)                                                                                 aTGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTATTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGCGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGt                                                                                                                                                             >  2:335888/1‑251 (MQ=255)                                                                                  tGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCAACGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAACAGCCTGGATGGAGCTTATCGCTGAAGTCGGGGGCGAAGGTCAGATTTTCTATCGccgccgccggc                                                                                                                                                            >  2:311898/1‑249 (MQ=255)                                                                                        aTTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCcgcg                                                                                                                                                      >  1:577540/1‑251 (MQ=255)                                                                                               ccggaagaTCGCACGGCGTTG‑ATCATACCTATCTGTAACGGAGGCGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAACCACCGGGCATGCCAAACACTTTGATGTCTACCTTCTTATGGACAGTTATAACTTAGATATCTGCGTCGCAGAGCAAAAAACCTGGATGGAGCTAAGCGCTGAAGTCGGTGGCGAGGGTCAGATTTTCTAAcgccgccgcagtacccgctggcagcg                                                                                                                                               >  2:293054/8‑234 (MQ=255)                                                                                                       tCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGCATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCATCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCg                                                                                                                                       >  1:322470/1‑251 (MQ=255)                                                                                                       tCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCg                                                                                                                                       >  2:240482/1‑251 (MQ=255)                                                                                                         gCACGGCGTTG‑ATCATGCCTATCTTTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                     <  2:500330/251‑1 (MQ=255)                                                                                                         gCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                     <  1:181311/251‑1 (MQ=255)                                                                                                          cACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                     >  1:499822/1‑250 (MQ=255)                                                                                                               gctgttaTCATGCCTATCTGTAACGAAGTCGTGAACCGTGCTTTTGCTGGCCTGCTTGCAACGGTGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTTATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAatat                                                                                                                                <  2:710559/245‑1 (MQ=255)                                                                                                                gTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCg                                                                                                                              <  1:335888/251‑1 (MQ=255)                                                                                                                 ttG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTATTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCACTGAAGTCGGTGTCGAAGGTCAGATTTTCTATCGCCGCCGTCGTCGCCGGGTGAAGCGTAAAAGCGGTCATATCg                                                                                                                              >  2:190707/1‑250 (MQ=255)                                                                                                                   g‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCCGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATg                                                                                                                           >  1:393541/1‑251 (MQ=255)                                                                                                                     aTCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGa                                                                                                                          <  1:183217/251‑1 (MQ=255)                                                                                                                       cATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTTGGAATCAGTAAAAGCCAACGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACt                                                                                                                        <  1:329795/251‑1 (MQ=255)                                                                                                                        aTGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGc                                                                                                                                                                                           <  2:25248/183‑1 (MQ=255)                                                                                                                        aTGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGc                                                                                                                                                                                           >  1:25248/1‑183 (MQ=255)                                                                                                                          gCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCt                                                                                                                     >  2:403099/1‑251 (MQ=255)                                                                                                                          gCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCt                                                                                                                     >  1:342863/1‑251 (MQ=255)                                                                                                                             tATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGcc                                                                                                                  >  1:550539/1‑251 (MQ=255)                                                                                                                                     aCGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTAGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTgg                                                                                                          >  1:56241/1‑251 (MQ=255)                                                                                                                                        aaGACGTGAACCGTGTTTGTGCTGGCCTGCGTGCAACGTGGGAAGCAATAAAAGCCACGGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGGCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGGCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGTCTTCTGCCGTCGCTGGGGc                                                                                                       <  2:442745/251‑1 (MQ=255)                                                                                                                                         agacggtaaCCGTGTTTTGGGTGGCCGACGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTTCATTCTTAGGGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGACAGATTTTCTATCGCCGCCGCCGTCGCCGCGGGAAGCGTAAAAGCGGTAATAACGATGACTTCTTCCGTTGCTGGGGCa                                                                                                      <  2:487829/244‑1 (MQ=255)                                                                                                                                          gACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGccg                                                                                                     >  2:4034/1‑249 (MQ=255)                                                                                                                                              tGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGcc                                                                                                   >  1:627961/1‑249 (MQ=255)                                                                                                                                                 aCCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTAc                                                                                              <  2:411184/251‑1 (MQ=255)                                                                                                                                                   cGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAg                                                                                            <  1:688252/251‑1 (MQ=255)                                                                                                                                                            gCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGcagatagg                                                                                    >  2:229178/1‑243 (MQ=255)                                                                                                                                                            gCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGc                                                                                           <  1:229178/243‑1 (MQ=255)                                                                                                                                                             cTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggtg                                                                                  >  1:264988/1‑251 (MQ=255)                                                                                                                                                                 ccTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGc                                                                              >  1:37133/1‑251 (MQ=255)                                                                                                                                                                         cAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCt                                                                      >  1:323943/1‑251 (MQ=255)                                                                                                                                                                           aCGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGa                                                                    <  2:171018/251‑1 (MQ=255)                                                                                                                                                                            cGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGAc                                                                   <  1:141192/251‑1 (MQ=255)                                                                                                                                                                             gTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACt                                                                  <  2:627961/251‑1 (MQ=255)                                                                                                                                                                              tGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTg                                                                                                                    <  1:230960/200‑1 (MQ=255)                                                                                                                                                                              tGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTgcacc                                                                 >  2:341859/1‑247 (MQ=255)                                                                                                                                                                              tGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGACAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTg                                                                                                                    >  2:230960/1‑200 (MQ=255)                                                                                                                                                                               gggAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTAGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCg                                                                <  1:582887/251‑1 (MQ=255)                                                                                                                                                                                      aGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATc                                                                                                                               >  1:555319/1‑181 (MQ=255)                                                                                                                                                                                      aGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATc                                                                                                                               <  2:555319/181‑1 (MQ=255)                                                                                                                                                                                       gTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGAc                                                        <  1:403099/251‑1 (MQ=255)                                                                                                                                                                                           aaGCCCCCGGGAATTCCAAACACTTGTATGTCGTCATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCGGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAAGATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCGGGATTCTGACTCGGTAATTACCGGt                                                    <  2:453586/251‑1 (MQ=255)                                                                                                                                                                                                   gggAATGCCAAACACTTTGATGTCTACCTTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCg                                                                                                                              >  2:367417/1‑169 (MQ=255)                                                                                                                                                                                                   gggAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCg                                                                                                                              <  1:367417/169‑1 (MQ=255)                                                                                                                                                                                                    ggAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAttgtttg                                           <  1:6714/251‑1 (MQ=255)                                                                                                                                                                                                      aaTGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCATCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGTTGCTGACTCGGTAATGACCGGGGCATGTTtgtg                                         >  2:578014/1‑251 (MQ=255)                                                                                                                                                                                                      aaTGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGGCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCt                                                                                                                     <  2:249954/175‑1 (MQ=255)                                                                                                                                                                                                         gCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑Cg                                      <  1:4034/250‑1 (MQ=255)                                                                                                                                                                                                         gCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGAAATGACCGGTGATTGTTTGTG‑Cgg                                     >  1:600543/1‑251 (MQ=255)                                                                                                                                                                                                              acacTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTgg                                >  1:72878/1‑251 (MQ=255)                                                                                                                                                                                                                acTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGATTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTg                              >  1:663069/1‑251 (MQ=255)                                                                                                                                                                                                                  tttGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTgcg                            >  1:606072/1‑251 (MQ=255)                                                                                                                                                                                                                   ttGATGTCTACATTCTTCGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCg                                                                >  1:573935/1‑215 (MQ=255)                                                                                                                                                                                                                    tGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGCTGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCCGTACAGCTACATGGTGGGGCTGGATGCTGACTCGGTAATGACCGGTGATTGATTGTG‑CGGGCTGGTGCGcc                          >  2:240707/1‑251 (MQ=255)                                                                                                                                                                                                                     gATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATCTCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACa                                                                                       >  2:206645/1‑190 (MQ=255)                                                                                                                                                                                                                     gATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCt                         <  2:264988/251‑1 (MQ=255)                                                                                                                                                                                                                     gATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATAGCGAGGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACa                                                                                       <  1:206645/190‑1 (MQ=255)                                                                                                                                                                                                                             cATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGGTGGa                   >  2:619105/1‑249 (MQ=255)                                                                                                                                                                                                                              catCTTAGAGACAGTTATCACGTGGATATATGCGTCGCCGAGCAAAAAGCCTGGATGGAACTTATAGCTGACGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGACGCCGCGTGGAGCGTAAAAGAGATAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAAGGTGGCGCGTGATGCTGACTAGGTAATGACCGGTGTTTGTT‑GTGGCGGGCAGGGGCGCCTGATGAAAGc                >  2:205651/3‑251 (MQ=255)                                                                                                                                                                                                                                   tAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAAc            <  2:165661/250‑1 (MQ=255)                                                                                                                                                                                                                                       gACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGaa       <  2:37133/251‑1 (MQ=255)                                                                                                                                                                                                                                       gACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGa        >  2:377776/1‑250 (MQ=255)                                                                                                                                                                                                                                         aaGTTATAACTTGGATATCTGCGTCGAAGAGAAAAAAGCCGGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACg     <  2:196239/250‑1 (MQ=255)                                                                                                                                                                                                                                            ttATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGTGCCTGATGGAAGCCAACCCGAccgccg  >  2:474750/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCg  >  1:84481/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                    ATTCTGGACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTG‑ATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTG‑CGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCG  >  NC_000913/1111486‑1111970

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