Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,529,185	C→T	94.9%	pseudogene (1284/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,529,185	0	C→T	94.9%	3.4 0.0 ‑0.4	97	pseudogene (1284/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  new base (39/51):  ref base (3/3):  total (42/55)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.91e-01

CCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATG  >  NC_000913/1529005‑1529482                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                             ccACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCAc                                                                                                                                                                                     <  2:550855/301‑1 (MQ=11) ccACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCAc                                                                                                                                                                                     <  2:307876/301‑1 (MQ=11) ccACGGCGTGGTGGTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGGGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCTGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCAc                                                                                                                                                                                     <  1:25857/301‑1 (MQ=2)           ggggtgcgggCTGGACGGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAACCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTAc                                                                                                                                                                            <  1:662653/292‑1 (MQ=11)            tgtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAACCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                                                                                                                                          <  2:277885/301‑1 (MQ=9)            ggtTCCCGGCTGAACTGGGGCGGCTGCGGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTAcc                                                                                                                                                                           <  1:548212/299‑1 (MQ=255)             gtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                                                                                                                                          >  2:744622/1‑300 (MQ=14)               tCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGAGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCTAAGACGGCAATACGGCGGGGCGCGGGGAGGATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGCACCCGCAGCCCCTGCCCCAGCGGTCCCGt                                                                                                                                                                         >  1:642410/1‑299 (MQ=11)                        aaCGGGGGCGGCGGCGGGAAAGGCGGGAGGAAGAAATCCGGGCAGACGGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGTAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgc                                                                                                                                                               <  1:1008929/300‑1 (MQ=11)                            ggTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGGCGGCAATACGGCGGGGCGCGGGGAGTATGATGAAT‑GGGGCCACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCCTCTGCCCGGgccgcag                                                                                                                                                          >  1:769078/1‑301 (MQ=11)                             ggtgggCTGGGGGACGGGCGTGAGGGAAGAAGCCGTGCAGGCGGCGTGAGCAGAGGAGGGCGGGCGGGGCTGGGTAAGAGCGGGCGTAAGGGGGGGGAACGGGACAGAAAGGGGGAGGGGGAATACAGAGCGGGGGGAAAGGATCATCGATATGACTGCGACCGACGGGGAGGGCCGCGGGCGCTTAGCAGCGAAGGCGGAAAAACGGCGGGGCGCGGGAAGTATGATGAAT‑GGGGCAACCAGCTGAATGAGGAGAGCGCGGATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcag                                                                                                                                                          <  1:163083/296‑1 (MQ=255)                             gTGCGGCTGCTGGCCAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGTAGCAACTGGCCCGACAGGTTGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTACCACTCCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAACACGGCGTCGCGCGGGGACTAGGATCGAT‑GGGGCAACCACCTTACTGAGGAGACCCCGCAGCCACTCCACCACCCGTCCCgccgccctggcacccag                                                                                                                                                          >  1:1006084/1‑284 (MQ=11)                                   ctgctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg                                                                                                                                                   >  2:476168/1‑301 (MQ=14)                                     gctgGGAAGGCGGGAGGAAGAAATCAGGGCAGACCGGGTGAGCATGGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCGGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                                                                                                                  <  1:301405/300‑1 (MQ=11)                                      ctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgagg                                                                                                                                                <  1:789646/301‑1 (MQ=14)                                              gCTGGAGGAAGAGATCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg                                                                                                                                         <  1:118256/300‑1 (MQ=255)                                               cTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAAAACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                       <  2:386562/301‑1 (MQ=255)                                                     gaagaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGGGTATGCTGAAT‑GGGGCAACCAGCTGAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGACCGGGCACCAGCATGCTGCGGAGTCGGTGCTGTAc                                                                                                                                  >  1:805984/1‑300 (MQ=11)                                                     gaagaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTa                                                                                                                                   >  2:140559/1‑299 (MQ=18)                                                     gaagaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc                                                                                                                                  <  2:1052642/300‑1 (MQ=18)                                                     gaagaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATAAGCGAAGACGTCAATACGGCGTGGGCCGGGGATTATGATGAAT‑GGGGCAACCAGCTTAATGCGGACAACCCGCATCACCTGCACACGCCGTCCCGTCTGCCAGGGCAGCAGCATGATGAGGGTTCGGGGCTGtg                                                                                                                                   >  1:1087234/1‑298 (MQ=11)                                                     gaaaaaATCCGGGAAGACCGCGGGAGCAGGGAAAGCGGGGCGTGGCTTGGGCAGGGCGGGCTGAGGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc                                                                                                                                  <  1:3402/300‑1 (MQ=11)                                                         aaaTCCGGGCAGACCGCGTGGGCAGGGAAAGCCGGGCGGGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATaa                                                                                                                             <  1:744600/301‑1 (MQ=255)                                                          aaaccGGGCAAACGGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTAGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAc                                                                                                                            <  2:257285/298‑1 (MQ=255)                                                          aaTCCGGGGAGACCGCGGTCGCAGTGAAAGCCGTGCGTGGATTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAACACCCACCGGATCGAACAGCTCATCTTTATCACTGCGGCCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAc                                                                                                                            <  1:309180/301‑1 (MQ=11)                                                            tCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCg                                                                                                                          >  1:721498/1‑301 (MQ=18)                                                             ccGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc                                                                                                                           <  1:140557/299‑1 (MQ=16)                                                                  cAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTGTGATGAAT‑GGGGCAACCAGATTAATGAGGAGAACCGGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCCGCAGCCTGCTGAGGAGTCGGGGCTGTACTATCACCGTCAcc                                                                                                                     >  1:841258/1‑300 (MQ=11)                                                                    gACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGGACTATAACCGTCACCg‑g                                                                                                                  >  2:422300/1‑300 (MQ=16)                                                                         cgTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTa                                                                                                             <  2:973483/300‑1 (MQ=25)                                                                         cgTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTa                                                                                                             >  2:484607/1‑300 (MQ=25)                                                                           ggAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTC‑GGGCTGTACTATAACCGTCACCG‑GCACTACGa                                                                                                          <  2:642421/299‑1 (MQ=255)                                                                            gAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGaa                                                                                                         >  2:140128/1‑300 (MQ=19)                                                                            gAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGa                                                                                                          <  2:1144948/300‑1 (MQ=25)                                                                               cAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGGGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATcc                                                                                                       <  1:725665/300‑1 (MQ=255)                                                                               cAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATcc                                                                                                       <  2:56093/300‑1 (MQ=25)                                                                                 gTGAAAGCCGGGCGTGGCTTTGGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATcc                                                                                                       <  1:609347/298‑1 (MQ=255)                                                                                 gTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGCGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGCCCGCAGCATGATGAGGAGTCGGGGCTGTCCTATAACCGTCACCG‑GCCCTACGATCCGtt                                                                                                    >  1:1066535/1‑301 (MQ=255)                                                                                 gTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATcc                                                                                                       >  2:609357/1‑298 (MQ=25)                                                                                  ttaaaGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCg                                                                                                      <  2:277346/296‑1 (MQ=255)                                                                                     aaGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑g                                                                                               <  2:140317/301‑1 (MQ=21)                                                                                           ggCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTCGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACAAACGGGGCCTGCCGCTGGCGCTTATAAGCGCAGTCGGCAATACGGAGCGGCGCGGGACGGATGATGAAT‑GGGGCAACCAGATTCATGAGGAGAACCAGCAGCACCTGCACCAGACGTACCGTCTGCCAGGGCCGCAGCATGAGTAGGAGTCGGGGCGCGCCTATAACCGTCACAGCGCA‑TTCGGTCCGTTGCA‑GGGGCgg                                                                                         >  1:475133/1‑301 (MQ=11)                                                                                                ggCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGCGTATGATGAATGGGGGCAACCAGCTTAATGAGGATAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTACCCG‑GCACTAGGAACCGTTGCATGGGGCgt                                                                                         >  1:559782/1‑297 (MQ=11)                                                                                                ggCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTatat                                                                                    >  2:877938/1‑301 (MQ=21)                                                                                                   ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAc                                                                                 <  2:294742/301‑1 (MQ=255)                                                                                                   ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAc                                                                                 <  2:808670/301‑1 (MQ=255)                                                                                                   ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAc                                                                                 <  2:379323/301‑1 (MQ=255)                                                                                                      cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAccc                                                                               <  1:748107/300‑1 (MQ=255)                                                                                                      cgcAGTGCGGGCGGACGGTGGAGCAACTGGCCAGACAGTGTGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGTGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCATGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAGCCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAccc                                                                               <  2:420566/300‑1 (MQ=11)                                                                                                          gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATGCGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCCACAGCCGTACCGTCTGCCCGGGCAGCAGCATGAGGAGGAGTCGGGGCTGGACTAAAACCGTCACCG‑GGATTACGAACCGTTGCG‑GGGGGCGTGTACCACCCCGGa                                                                           >  1:13416/1‑300 (MQ=11)                                                                                                          gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGG‑GCCGGAATCCACACCGGCGCGGAAAGCTCATCTTTATCACTGGCACCAACGGGGGCTGGCCCCGGCGCCTCTCCGCCCAAACGGCAAAACCGCGGGGCCGGGGGAGGATGGTGAAT‑GGGGCAACCAGCTTAATGAGGAAAAACCGCAACACCTGCACCAGCCGTACCGTCTCCCAAGGCCGCCGGCTGATAGGGAATCGGGGGTGTACCATGACCTGCAACG‑GCACTACGGACCGTTGCA‑GGGGCGGGATTTACCACCGGAc                                                                          >  1:470875/1‑300 (MQ=2)                                                                                                            gCGGGCTGACGGTGGAGCAACTGGCCAGAAAGGTGGAGCCGGAATCCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGAcc                                                                         <  1:167014/300‑1 (MQ=255)                                                                                                                    ccggggggccAACTGCCCAGACAGGGGGACCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGt                                                                <  1:435106/292‑1 (MQ=14)                                                                                                                       gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGCACTATCACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTgg                                                             >  2:235972/1‑300 (MQ=255)                                                                                                                       gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCGGTTGCC‑GGGGCGGGATATCACCCCGGACCCGAGTGGGtggg                                                             >  1:64864/1‑298 (MQ=255)                                                                                                                       gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTgg                                                             >  2:87639/1‑300 (MQ=255)                                                                                                                          gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTgaga                                                           <  2:241394/300‑1 (MQ=255)                                                                                                                            gCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAgg                                                         <  1:726572/300‑1 (MQ=255)                                                                                                                            gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAgg                                                         <  2:960178/300‑1 (MQ=255)                                                                                                                              aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgg                                                      <  2:25984/301‑1 (MQ=255)                                                                                                                                  ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGAACCGCTTGGGTTGGGAGGgggatgg                                                  >  1:832920/1‑301 (MQ=255)                                                                                                                                     cggACAGGTGGGGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAG‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAat                                               <  1:502999/299‑1 (MQ=255)                                                                                                                                       gACAGGTGGAGCCGGAATACACACCGACGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAACCGTACCGGCTGCCAGGGCAGCAGCATGAGGAGGAGTCGGGGCTGCACTATAACCGTCACAG‑GCACTACGATCAGTTGCA‑GGGGCGGACTAGCAACCAGGGCCCGGTTGGGTTGAGAGGTGGCTGGAatat                                             >  1:212027/1‑301 (MQ=14)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGCTGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGGTCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGGGGGGGAAGGAATTTGTAt                                         >  1:1126552/1‑301 (MQ=14)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATaa                                                                                                                             <  2:909024/219‑1 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑gggg                                                                                            <  2:680798/250‑1 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑gggg                                                                                            >  1:680785/1‑250 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                                          <  2:222075/300‑1 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCAGTTGCA‑GGGGCGGTATATCACCCAGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTAt                                         >  2:669669/1‑301 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAGGATGAGGAGTCGGGGCTGTACTATaa                                                                                                                             >  1:908996/1‑219 (MQ=255)                                                                                                                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCCGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCC‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGGATATGt                                           >  1:530448/1‑299 (MQ=255)                                                                                                                                           ggTGGACGCGGAACACACACCGGGGGGGAAATCGCATCTTGAGCAGCGCGACCACCGGGGAGTGCCGCTGGCGGTTAGCAGCGAAGACGGCAAGACGGCGTGGCGCGGGGAGGATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGGACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                                          <  1:138115/300‑1 (MQ=14)                                                                                                                                            gTGGAGCCGGAATACCAACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                                          <  2:1066576/299‑1 (MQ=255)                                                                                                                                            gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATGTACCCCCGGACCCGGATGGGTTGGGGGGGGGGAGGAATAGGTa                                          >  1:222073/1‑299 (MQ=14)                                                                                                                                            gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTa                                          <  2:1144651/299‑1 (MQ=255)                                                                                                                                            gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTAt                                         <  2:769101/300‑1 (MQ=255)                                                                                                                                                aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGtt                                    >  2:754674/1‑300 (MQ=255)                                                                                                                                                aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTa                                    >  2:79465/1‑301 (MQ=255)                                                                                                                                                  ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGCACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTa                                    <  2:920122/299‑1 (MQ=255)                                                                                                                                                  ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATc                                  <  2:1126598/301‑1 (MQ=255)                                                                                                                                                  ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATc                                  <  2:848432/301‑1 (MQ=255)                                                                                                                                                    ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGCTTGGGTTGGGAGGTGGATGGAATATGTATCAGTATCCg                                >  1:498821/1‑301 (MQ=255)                                                                                                                                                       aTACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGCTTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTg                             >  2:376307/1‑301 (MQ=255)                                                                                                                                                       aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTCAGATCCGTTGCC‑GGGGCGGTATATCACCCCGGACCCGAGTGGGTTGAGAGGTGGGTGGACTATGTATCAGTCACCGTTg                             >  1:572273/1‑301 (MQ=14)                                                                                                                                                       aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAg                                      <  2:898955/292‑1 (MQ=255)                                                                                                                                                       aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGGTGAGAGGGGGATGGAATATGTATCAg                                      >  1:898928/1‑292 (MQ=255)                                                                                                                                                       aTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTAGGGGCTGTACTATAACCGTCACCG‑GAACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTAGCAGTATCCGTTg                             >  2:438147/1‑301 (MQ=255)                                                                                                                                                            cacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATcc                        <  2:721512/301‑1 (MQ=255)                                                                                                                                                                ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGGATATCACCCCGGACCCGATTGGGTTGGGAGGTGGATGGAATATGTAACCGTATCCGTTGAATCCCATa                    >  2:1105557/1‑301 (MQ=255)                                                                                                                                                                  cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                   <  1:393450/300‑1 (MQ=255)                                                                                                                                                                    cgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTGTCAGTATCCGTTGAATCCCATACAAg                <  1:204501/301‑1 (MQ=255)                                                                                                                                                                    cgAAAAGCTCAGCTTTATCCCTTCGACCACCGGGGACTGCCGCTGGCGCTTATCAGGGAAGACGGCAAGACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAg                <  1:335822/301‑1 (MQ=255)                                                                                                                                                                         agctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGGTGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTTATa           >  2:864258/4‑301 (MQ=255)                                                                                                                                                                            cgagctttATCACTGCGAGCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAc        <  1:134499/296‑1 (MQ=255)                                                                                                                                                                              atctttACCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAccc      <  1:422294/298‑1 (MQ=255)                                                                                                                                                                                ctttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATGGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCaa    >  2:283833/2‑301 (MQ=255)                                                                                                                                                                                  ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGCGGTGGATGGAATATGTTACAGTATCCGTTGAATCCCATACATGTGATAGACCCAATg  >  2:911070/1‑301 (MQ=255)                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                             CCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATG  >  NC_000913/1529005‑1529482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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.