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

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,529,175	0	T→C	100.0%	3.3	92	pseudogene (1274/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  new base (37/50):  ref base (2/2):  total (39/53)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.34e-01

TATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAG  >  NC_000913/1528876‑1529468                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         tATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAACAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGGAAAGtt                                                                                                                                                                                                                                                                                                         >  2:548220/1‑301 (MQ=11)                                                                                                                                  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:379323/301‑1 (MQ=255)                                                                                                                                                                                                                                    ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAc                                                                   <  2:294742/301‑1 (MQ=255)                                                                                                                                                                                                                                    ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCAc                                                                   <  2:808670/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                                                                              >  1:680785/1‑250 (MQ=255)                                                                                                                                                                                                                                                                            ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑gggg                                                                              <  2:680798/250‑1 (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:848432/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                   ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATc                    <  2:1126598/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)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         TATCCAGACGGTATACGAGCCGGGAAGCTTCACGCCGCTCATCCGGGTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAT‑GGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCG‑GCACTACGATCCGTTGCA‑GGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAG  >  NC_000913/1528876‑1529468

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.