Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,938,768	IS2 (+) +5 bp	coding (542‑546/993 nt)	rbsR →	DNA‑binding transcriptional dual regulator RbsR

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3187426	NA (NA)	97 (1.010)	83/494	0.1	96.9%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3938768 =	3 (0.030)					coding (542/993 nt)	rbsR	DNA‑binding transcriptional dual regulator RbsR
*	?	NC_000913	= 3938772	3 (0.030)	106 (1.100)	89/494	0.1	97.2%	coding (546/993 nt)	rbsR	DNA‑binding transcriptional dual regulator RbsR
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GGTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATT  >  NC_000913/3187178‑3187674                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           ggTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑t                                                                                                                                                                                                                                                         >  1:178196/1‑251 (MQ=11)   ggCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTa                                                                                                                                                                                                                                                       <  2:724461/250‑1 (MQ=14)    gCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTat                                                                                                                                                                                                                                                      <  1:338702/251‑1 (MQ=14)     gggAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTatc                                                                                                                                                                                                                                                     <  2:966310/250‑1 (MQ=14)      ggAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTatca                                                                                                                                                                                                                                                    <  2:978296/251‑1 (MQ=32)           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTTATAACAG‑ATGTCTGGAAATATAGGGGCAAACCCAG‑TTATCATCAtt                                                                                                                                                                                                                                               <  1:753540/251‑1 (MQ=21)           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCAtt                                                                                                                                                                                                                                               >  1:183858/1‑251 (MQ=34)           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCAtt                                                                                                                                                                                                                                               <  1:505392/251‑1 (MQ=34)           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGAATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCAtt                                                                                                                                                                                                                                               >  1:394183/1‑251 (MQ=21)                    aaCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTCCATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCAAG‑TTAACATCATTGGTTGTGCg                                                                                                                                                                                                                                      >  2:975473/1‑251 (MQ=255)                     aCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGt                                                                                                                                                                                                                                     >  1:891817/1‑251 (MQ=39)                          aTAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCt                                                                                                                                                                                                                                <  1:615884/251‑1 (MQ=255)                          aTAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGCGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGTGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAATTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCt                                                                                                                                                                                                                                <  2:98679/251‑1 (MQ=255)                            aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTa                                                                                                                                                                                                                              <  1:608132/251‑1 (MQ=255)                                gagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGcc                                                                                                                                                                                                                          <  1:602572/251‑1 (MQ=255)                                   agCTTCGTGAAAACGAGAAAGCGTGAACACATCAGTATCATGCCCTAACCAGGTGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGGTCGAGGATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATAGCTTCGGCAGCGGGCTTGAAATGGGTTAAGGGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTAAGCGTATGCCGt                                                                                                                                                                                                                        <  2:924220/250‑1 (MQ=14)                                     cTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTg                                                                                                                                                                                                                      >  1:228379/1‑250 (MQ=255)                                      ttCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGa                                                                                                                                                                                                                     >  2:797626/1‑250 (MQ=255)                                         gTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGGGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTc                                                                                                                                                                                                                 <  1:921264/251‑1 (MQ=255)                                                 gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGa                                                                                                                                                                                                         <  1:24782/251‑1 (MQ=255)                                                 aaTAAAGCGGTACTACATCAGTATCATGCCCAAACAAGAAGGTTTAACGGCAGCAAAGAACCTTGTAGAGTCGGTCGAGCACTATAACGAATGGCATCCGCATAGTGCGCTGGTTTATCTCTCGCCACGGGAATATCTTCGGCAGCGGGTTTGTAATTGGTTAATTGATAACAG‑ATGTCTTGAAATATAGTGGCAAATCCAG‑TTATCATCATGGGTTTTGCGACATCTGATGTCGTTGAATTCAACAAAGa                                                                                                                                                                                                         <  2:885333/250‑1 (MQ=14)                                                  aTAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAT‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGAt                                                                                                                                                                                                        <  1:811293/251‑1 (MQ=255)                                                  aTAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGAt                                                                                                                                                                                                        <  2:804949/251‑1 (MQ=255)                                                        cGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATc                                                                                                                                                                                                  >  1:202924/1‑251 (MQ=255)                                                         gTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                                                                                                                                                 <  1:470725/251‑1 (MQ=255)                                                          tGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGCATTCAACAAAGATTTAATCg                                                                                                                                                                                                 >  2:791653/1‑250 (MQ=255)                                                             cTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGc                                                                                                                                                                                             >  2:136122/1‑251 (MQ=255)                                                                cATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCg                                                                                                                                                                                           >  1:886624/1‑250 (MQ=255)                                                                      ggtCATGCGCAAACCCGGCGGGTTAACGGCAGTAAAGAACCTTGCAGAGGCGGTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGTCATGGGAATATCTGAGAAAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGACAGCTTATGAGGTTGAATTGAACAAAGATGTAATCGAAGCCGAAGATct                                                                                                                                                                                    <  2:200177/249‑2 (MQ=17)                                                                       aTCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGTGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCg                                                                                                                                                                                    <  1:128889/250‑1 (MQ=255)                                                                       aTCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGa                                                                                                                                                                                          >  2:303381/1‑244 (MQ=255)                                                                        tCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTg                                                                                                                                                                                  >  1:148353/1‑251 (MQ=255)                                                                            gCCCAAACCAGACGTGTTAACGGCAGCAAAGAACCTTGCAGAGGCGGTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAAGATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGaa‑a                                                                                                                                                                              <  2:492153/250‑1 (MQ=255)                                                                            gCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGCAGATCGTGaa‑a                                                                                                                                                                              >  1:63863/1‑250 (MQ=255)                                                                                 aaCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTa                                                                                                                                                                        <  1:77943/251‑1 (MQ=255)                                                                                 aaCCAGACGGGGTAACGGCAGCCAAGACCCTTGCAGAGGCGTTCGAGCATTATAACGAAGTGCATCCGCATAGTGCGCTGGGTTATCGCTAGCCACGGGAATATCTGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑ATATCATCATTGGTGGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGACGCCGAAGATCGTGAA‑TACGTTa                                                                                                                                                                        <  2:672301/251‑1 (MQ=18)                                                                                          ggTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTccc                                                                                                                                                               <  2:376192/251‑1 (MQ=255)                                                                                          ggTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCAGAAGATCGTGAA‑AACGTTAACCTTTccc                                                                                                                                                               >  1:612467/1‑251 (MQ=255)                                                                                             tAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCaa                                                                                                                                                             <  2:934401/250‑1 (MQ=255)                                                                                              aaaggCAGCAAAGACCCTTGCAGGAGCGTGCGATCATTCTAACTAATTGTCGCCACATTGTGCGCTGGGTTATCGCTCTCCTCGTGAATATCTGAGGAAGCGGGCATGGAATGGGTTAAGTTATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAAt                                                                                                                                                            <  2:18088/247‑1 (MQ=21)                                                                                               aCGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAAttt                                                                                                                                                          <  1:73054/251‑1 (MQ=255)                                                                                                cGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTg                                                                                                                                                         <  1:775611/251‑1 (MQ=255)                                                                                                   cagcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCCACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGaa                                                                                                                                                       >  1:532002/1‑250 (MQ=255)                                                                                                   cagcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGaaa                                                                                                                                                      <  1:887230/251‑1 (MQ=255)                                                                                                    agcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGAGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAAc                                                                                                                                                     >  1:616788/1‑251 (MQ=255)                                                                                                     gcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGGGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCAGTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACt                                                                                                                                                    <  2:578814/251‑1 (MQ=255)                                                                                                                   agagGCGATCGAGTATTATAACGAATGGCATCCGCAGAGTGCGCTGGGTTATCGCTCGACACGGGAATATCTGTGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGGCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGACAATTAcc                                                                                                                                      <  2:921671/251‑1 (MQ=25)                                                                                                                    gagGCGTTCCAGCATTATAACGAATGGAATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AGCGTTAACCTTTCCCAATTTGAAACTGATGGGCAATTAccc                                                                                                                                     <  2:833897/251‑1 (MQ=255)                                                                                                                    gagGCGTGCGAGCAATATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCGCGCCACGGGAACATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTAGCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTAccc                                                                                                                                     <  2:95068/251‑1 (MQ=255)                                                                                                                            cGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCa                                                                                                                              >  1:655173/1‑250 (MQ=255)                                                                                                                             gAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCaa                                                                                                                             >  1:655023/1‑250 (MQ=255)                                                                                                                               gCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAat                                                                                                                           >  2:47318/1‑250 (MQ=255)                                                                                                                               gCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAat                                                                                                                           >  1:408081/1‑250 (MQ=255)                                                                                                                                   tataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTCAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTa                                                                                                                                                                                                                              >  2:443014/1‑148 (MQ=255)                                                                                                                                   tataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTa                                                                                                                                                                                                                              <  1:443014/148‑1 (MQ=255)                                                                                                                                        cGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTa                                                                                                                   >  2:621773/1‑249 (MQ=255)                                                                                                                                         gAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGc                                                                                                                <  1:791653/251‑1 (MQ=255)                                                                                                                                          aatggaagccgaaTAGTGCGCTGGGATATCGCCCGCCACGGGAATATCGGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCa                                                                                                               <  2:228379/239‑1 (MQ=38)                                                                                                                                                   cGCATATTGCGCGGGGTTATCGCTCGTCAAGGGAAAATCTGCGGTAGAGGGCGTGTAAGGGGTTAAGTGATAACAG‑ATGTCTGGAAAGATAGGGGCAAATCCAG‑TTATCATTATTGGTTGTGCGTCAGCTTATGCCGTTGAAATCAACAAATATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTTAAACTGATGGCCAATTCCCCGTCGGCAAATATTCACTAAGCACTCTGagga                                                                                                      <  2:356797/251‑4 (MQ=255)                                                                                                                                                      aTAGTGCGCTGGGTTATCGCTCGCCACGGGCATTTCGGCGCCGGCGG‑CTTGGAAAGGGTTAAATGATAACAGGATGCCTGGAAA‑ATTGGGGCAAATTCAGGTTATCA‑CCTTTGTTGTGCGTTCGCTTAAGGCGTTGAATTTCACCAAGGTTTAATCGGAGCCCGAGGTCCTGAACACCGTTAACCTTTCC‑AATTTTGAACTGATTGTCACTTACCCGGTGGCCAATATTCACTAAGGCCTCTGaataatt                                                                                                   >  1:302450/1‑249 (MQ=255)                                                                                                                                                        aGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGTGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGTCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATtaataat                                                                                                 <  2:721459/251‑1 (MQ=255)                                                                                                                                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAaga                                                                                             <  1:840354/251‑1 (MQ=255)                                                                                                                                                             gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGACGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATAACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAag                                                                                              >  2:760282/1‑249 (MQ=255)                                                                                                                                                                    tATCGCTCGCCACGGGAATATCTGCGACAGCGGGCTTGTAATTGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGcc                                                                                      <  2:672371/250‑1 (MQ=255)                                                                                                                                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGCTAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCAcc                                                                              >  2:476968/1‑251 (MQ=255)                                                                                                                                                                            gcCACGGGAATATCTGCGGCAGCGGGCATGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTGGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCt                                                                             <  2:655023/251‑1 (MQ=255)                                                                                                                                                                              cACGGGTATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACACGTCGGCAAATATTAACTAAGCACTCTGATTCATAATAAGAGGACGCCAATCCACCTTg                                                                           >  2:431994/1‑251 (MQ=255)                                                                                                                                                                               aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGa                                                                          <  1:542961/251‑1 (MQ=255)                                                                                                                                                                                         cTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATgg                                                                >  1:925409/1‑251 (MQ=255)                                                                                                                                                                                            cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATggga                                                              >  2:4256/1‑249 (MQ=255)                                                                                                                                                                                                 gCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGCCGCCAATCCACCTTGACCTCCAATGGGTATTAAt                                                        >  1:685403/1‑251 (MQ=255)                                                                                                                                                                                                         gTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcc                                                <  1:883380/251‑1 (MQ=255)                                                                                                                                                                                                         gTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAc                                                 <  2:408081/250‑1 (MQ=255)                                                                                                                                                                                                           aaTGTGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGACGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcct                                               <  2:183858/250‑2 (MQ=255)                                                                                                                                                                                                              gggTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACt                                           <  2:886624/251‑1 (MQ=255)                                                                                                                                                                                                                 ttAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGc                                         >  1:890029/1‑250 (MQ=255)                                                                                                                                                                                                                  tAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGc                                         >  2:276076/1‑249 (MQ=255)                                                                                                                                                                                                                    aGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                     <  1:880207/251‑1 (MQ=255)                                                                                                                                                                                                                     gTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGttt                                    <  1:175134/251‑1 (MQ=255)                                                                                                                                                                                                                     gTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGttt                                    <  1:708580/251‑1 (MQ=255)                                                                                                                                                                                                                      tGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACAAAACTGCAGTTTg                                   >  1:674557/1‑251 (MQ=255)                                                                                                                                                                                                                                aTGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGCCAcc                          >  2:725787/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACAcc                          >  1:221505/1‑249 (MQ=255)                                                                                                                                                                                                                                      ggAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGc                     <  1:276076/250‑1 (MQ=255)                                                                                                                                                                                                                                      ggAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCa                    <  1:841803/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATAGTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCTCCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGc                     <  1:922475/250‑1 (MQ=255)                                                                                                                                                                                                                                        aaaTATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGCAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGCCACCAGAGCAAt                  >  1:358657/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaTATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAAtt                 <  1:520705/251‑1 (MQ=255)                                                                                                                                                                                                                                            attggggCAAATCCAT‑TTATCATCATTGGTTGGGCGTCAGCTCATGCCGTTGATTTCAACAAAGATTTAAACGAAGCCGAAGAGCGTCAA‑ACCGTTACCCTGTCCCAATTTGAAACTGATGGCATATTGCCCGTCGGCAAATATGCACTAAGCACTCGGATTAATAATAAGAGGATGCCAATCCACCTTGACCTCCAATGGGTATTAGGTGATAACCAAACTGCAGTTTGCGTGACACCAGTGCAATTaa               <  2:448016/247‑1 (MQ=255)                                                                                                                                                                                                                                                 tgCAAATCAAG‑ATATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGCATTCCACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGAcc                                                                        >  2:84672/2‑188 (MQ=255)                                                                                                                                                                                                                                                 ggCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGAcc                                                                        <  1:84672/188‑1 (MQ=255)                                                                                                                                                                                                                                                   cAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTACCattatt        >  2:452856/1‑250 (MQ=255)                                                                                                                                                                                                                                                    aaaTCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCAGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattattcg      >  2:44114/1‑249 (MQ=255)                                                                                                                                                                                                                                                     aaTCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAgg     >  1:851328/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ccAG‑TTATCATCATTGGTTTTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCATAGCAATTAACATTATTAGGAtt  <  1:822205/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           GGTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAG‑ATGTCTGGAAATATAGGGGCAAATCCAG‑TTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAA‑AACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATT  >  NC_000913/3187178‑3187674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.