Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,523	G→C	7.6%	intergenic (+130/+585)	eco → / ← mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	2,304,523	0	G→C	7.6%	5.0 ‑0.8 0.0	71	intergenic (+130/+585)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  new base (5/2):  ref base (34/30):  total (39/32)
Fisher's exact test strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTCCGAGCCGGT  >  NC_000913/2304255‑2304808                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                               caacgcgGATCTGGGCGATGCTGGAATGTGGCGTTAAAACAGGAAGCTGCCGATCGTGGTGTGTACGCCAGACAATGTATGTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAAGTGACAACGGGGTAGTTCGCTAAACTGCCGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAAAGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAcgccg                                                                                                                                                                                                                                                                <  2:138461/298‑1 (MQ=25)   ccGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACC‑CCGC‑TGTGAAGGGCGGCACCCCGTAggctggctcttatacccatctgacg                                                                                                                                                                                                                                                              >  2:260369/1‑278 (MQ=255)   ccGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAgg                                                                                                                                                                                                                                                                                     <  1:260371/278‑1 (MQ=255)                  aTGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGc                                                                                                                                                                                                                                               <  1:874499/301‑1 (MQ=255)                   tGCTGGAATGCTGCGTTACAACAGCAAACTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCCTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGGCTCTAGCCCGTATTATCAG‑CCCTACA‑CCAC‑CGTGAAGTGCTGCACCCCAGAGGTCAGACAAgtcagccaacacccacccgccctcccacgc                                                                                                                                                                                                                                               >  2:27908/1‑271 (MQ=25)                     cTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCg                                                                                                                                                                                                                                            >  2:317340/1‑301 (MQ=255)                                                 gCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCTACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCTGCACCCCGTAGTTCAGACAACGCTGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGGCGCTGgcc                                                                                                                                                                                                                 >  2:457589/1‑299 (MQ=255)                                                     aTCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTTCGGCACCCCGTAGGTCAGACCAGGCGGTCACCCCCCATCCGACATCCAACGCCCGCGCCGGTTGCCTGCTGCGACGCTGgccccgc                                                                                                                                                                                                             >  2:409521/1‑295 (MQ=255)                                                      tCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAG‑GCCTACACCCAC‑TG‑GAAGTGCGGCACCCCGCCGGGCAGACAAGGCGTCAACGCCGCACACGACATCCCACGCCCGAGCAGGTTCCCTGATGCGACGCTgtccggcatt                                                                                                                                                                                                           >  2:566377/1‑292 (MQ=21)                                                             gtATACGCCAGACAATGTAGATGTGATGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGCCACCCCGTAGGTCAGACACGGCGGTCACGCCGCATCCGCGATCCAACGCCCGACCCGGTTTCCTGATGCAACCCTGGCGCCTCTTATCAgg                                                                                                                                                                                                     >  2:175648/1‑300 (MQ=255)                                                             gtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGACGTGCTGCACCCCGTAGGTCAGACAAGGCGGCCACGCCGCATCCGAAATCCAACGCCCCAGCCGGTTGCCTGATGCGACGCTGTCCCGTCttcacaag                                                                                                                                                                                                     >  2:721880/1‑294 (MQ=255)                                                             gtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCCgg                                                                                                                                                                                                     >  1:605386/1‑300 (MQ=255)                                                             gtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                                                                                                                                                                                     >  1:592040/1‑300 (MQ=255)                                                                   cccAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                              <  1:317343/300‑1 (MQ=255)                                                                    ccAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                              <  1:44912/300‑1 (MQ=255)                                                                      aGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCg                                                                                                                                                                                            >  1:603860/1‑300 (MQ=255)                                                                       gACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCTCATCCGACATCCAACGCCCGAGCCGGTTGCCTGGTGCGGCGCTGGCGAGTCTGTTCAGGCCTACACCGc                                                                                                                                                                                           >  2:751358/1‑300 (MQ=255)                                                                        aCAATGTAGATGGGAATTACCGCGTCTGGAAGACGGAAGGGAAAATTGACAACGCGGTATTTCGCTAAACCGCCGTGAAGGGCGGCACCCCTTAGGGCAGACAAGGCGGTCACGCCGCATCCGACATCCGACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGc                                                                                                                                                                                           <  2:592023/299‑1 (MQ=255)                                                                              tAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                                                                                                                                    <  1:625903/300‑1 (MQ=255)                                                                                      aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACATCGCGGTAGTTCGCTAAACTGCCGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGCCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAc                                                                                                                                                                            <  1:630224/300‑1 (MQ=255)                                                                                               gtagggAGGGCGAAAGAGACAATTGACAACGCGGCGGGGTGCGAAACGGCGGTGAAGGTCGGCAGCCGGTATGTCAGACAAGGCGGTCACGCCGCATGCGACACCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTAAA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                                                                                                                                                                   <  2:193736/296‑1 (MQ=14)                                                                                                               agtaaaTTGACAACGCGGAAGTGAGATAAAGTGCCGAGAAGTGTGGCGCCCCGTATGTCAGACAAGGCTGTCACGCCGCATCTGACATCCAACGCCGGAGCCGGGTGCCTGAGGCGACGATGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGTTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTCGGATAAGGCGGTTAc                                                                                                                                                   <  1:744764/297‑1 (MQ=14)                                                                                                                  aaaTTGACAACGCGGTAGTTCGCTAAACTGCCGGGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                                                                                                               <  1:751386/301‑1 (MQ=255)                                                                                                                            cgcgGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTCGCCTGATGCGACCCTGGCGCCTCTTATCAGGCCTAAGCCGTCGGGCAGTGCTCACCCCCGTAGGACGGACAAGGCCGttcacgcacaaccgacaa                                                                                                                                     >  2:649905/1‑285 (MQ=14)                                                                                                                              cgGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACC‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGCCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGCTCGCACAAGGCGGTTACGCCGCATCCGCCac                                                                                                                                     >  1:461635/1‑298 (MQ=255)                                                                                                                                    ttCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGcc                                                                                                                             <  1:649929/301‑1 (MQ=255)                                                                                                                                    ttCCCTAAAATGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGGCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                                              <  1:566394/300‑1 (MQ=255)                                                                                                                                          acaccgctgTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCTCCACCCCCTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTACCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCCCTGTCGAGTGCCCCCACGCCTACGCCTGCTAACGCGGTTACCACGCCTCCgcccgccaccgccctagc                                                                                                                        >  2:744631/9‑283 (MQ=11)                                                                                                                                          aaaCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                                                                                                                                    >  1:865303/1‑240 (MQ=255)                                                                                                                                          aaaCCGCCGTGAAGTGACGCAACCCGTAGGTCAGACAAGGCGGTGACACCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGTTGGCGGGTCTTATCAG‑GCCTACA‑CCGC‑GGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                                                                                                                                    <  2:865270/240‑1 (MQ=14)                                                                                                                                                       gTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCCTAGGTCGGATAAGGGCGGTACGCCCCCTCCCACATCCCCCCCCCCACCCGGTTGccctcgg                                                                                                           >  2:195570/1‑295 (MQ=11)                                                                                                                                                            gCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGCCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCCCGTCTTATCAGGCCTACACCGCTTGGAAGTTCTCCACCCCGGAGGGAGGATAAGGCGGTTACGCCGCATCCGACCTCCACCGCCCCAGCCCGTTGCCTgcagccgcg                                                                                                      >  2:243376/1‑292 (MQ=11)                                                                                                                                                              aCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTGCGCCGCATCCGACATCCAACGCCCCAGCCGGTTCCCTGCTGCGCCg                                                                                                      >  2:826031/1‑298 (MQ=255)                                                                                                                                                                 ccGTAGGTCAGACAAGGCGGTCACGCTGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGGGCTCCACCCAGCACACCGGACAATGCGGACACGCCGCATCCGGCATCTCCCGCGCATCACGGCTGCCTTAGACGACGCTGgcg                                                                                                >  2:174673/1‑301 (MQ=14)                                                                                                                                                                 ccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTac                                                                                                                                                                                                >  2:610522/1‑205 (MQ=40)                                                                                                                                                                 ccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTac                                                                                                                                                                                                <  1:610540/205‑1 (MQ=40)                                                                                                                                                                 ccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCAACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTTCGCCGCATCCGACCCCTCAGCCCCAAGCCGTTTGGCCGCTCCGACGCttccg                                                                                                >  2:653787/1‑297 (MQ=14)                                                                                                                                                                     aGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCCACGCTGGCGCGt                                                                                             >  1:943802/1‑300 (MQ=255)                                                                                                                                                                      ggTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                             >  1:60000/1‑299 (MQ=255)                                                                                                                                                                      ggTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTCACGTGCGCCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGACGGTTGGCTGATGCAACGCTGGCGCGGCTTATAAGGCCCACACCGCTGTGATATGATCGACCCAGTCGGGACGACAAAGAGGCTCCGACGCATGCGACACCACACGCCCGAGCACGTTgaccgaatcagcgccggccacgc                                                                                            >  2:337250/1‑278 (MQ=2)                                                                                                                                                                                    cgtcccccgcATCTCGCATCCGACGCCCGCGGCGGTTGCCTGGTGGGTCGCGGGCGCGACTGA‑CAGCGCGTACA‑CCGG‑TGTGAAGGGCGGCACCCCGTAGGTAAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGGCTGATGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                             <  2:461628/295‑1 (MQ=14)                                                                                                                                                                                     ctgcgccgcATCCGACGTCCTCGTCCCGAGACGGTAGGCGGGTGCGACGTTGGAGCGTCGGATTGG‑GCCTACA‑CAGC‑TGGGAAGTCCGGCACCCCGCAGGTGTGGTAAGGCGGTCACGCCGCATCCGACATCAAACTCCCGAGCCGTTGGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                             <  2:745650/297‑1 (MQ=11)                                                                                                                                                                                     cagcgacccATCGGACAGCCAAGGGCGGAGCCGGGTGCGAGATGCGACGCTGGCGGGGCTGATCAG‑GCCTAAA‑CCGC‑TGTGACGGGAGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACACCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                             <  2:823583/292‑1 (MQ=14)                                                                                                                                                                                        cgccgcacgggtggcacaacgacctgggcGGTTGCCTGGTGGGAGCGTGGCGCGGCTTATCAG‑GCCTAAC‑CGGC‑CGTGAAGGGAGGCACCCCGTAGGTCCGGCAAGGCGGTCACGCCGCATCGGCCATCCACCGCCCGAGCCGGTTGCCTGATGGGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCAGGCGCGTCTTATCAGGCCTACACCGc                                                                          <  2:174337/272‑1 (MQ=11)                                                                                                                                                                                              ctCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                   <  1:174674/300‑1 (MQ=39)                                                                                                                                                                                              aTCCGACATCCAACGCCCGAGCCGGTTGCCGTAGGCGAGGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                   <  1:175649/301‑1 (MQ=255)                                                                                                                                                                                                        cAACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTGATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAccc                                                         <  1:94881/301‑1 (MQ=39)                                                                                                                                                                                                             cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAg                                                    <  1:457596/301‑1 (MQ=38)                                                                                                                                                                                                                      gTTGCTTGAGGCGAGGTGGGCCCGTCTTCTCAG‑GGCT‑CA‑CCCCGTGTGAAGTGCGCCACACGGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATa                                            <  2:60000/300‑1 (MQ=11)                                                                                                                                                                                                                      gTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATa                                            <  1:653811/300‑1 (MQ=38)                                                                                                                                                                                                                      gTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTAAA‑CGGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGACCGAGCGTGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGTTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCTACCCCGTAGGTTGGATa                                            <  1:27908/300‑1 (MQ=255)                                                                                                                                                                                                                                                 tATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTCACGCCCAAGCCGGTTACCTGATGCACCGCTGGCGCGTCCTATCTGGCCTAACCCGCTGTGACATGCCCAACCCCCTACGCAGGCTCACGCGTTTACGCCGCCTCCG‑ACATCGaa               >  2:591878/1‑301 (MQ=14)                                                                                                                                                                                                                                                 tATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTGCGCCGCATCCG‑ACATCCaa               >  1:17185/1‑301 (MQ=39)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGATAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCTCGTCTTGTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCTAGGTCGGATAAGACGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTTATGCAACACAGGCGCGTCTTATCCTGCCTTCACCGCGCCGCAGTTCCCCCCCCCTAAGTGTAGATAAGGCGGGTACGACCCACCCGCACACCCAGCGt            >  2:415545/1‑300 (MQ=11)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGGCGCCGGCGCCTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTACGTCGGATAAGGCGGTTACGACGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGCTGCGGCGTCGGCTCGTCTTATGCAGCCTACAACGCGGTGCAGTGCTCCCCCCAGCTGTGTGGCACAGGCGGATACGACGCCGCCG‑CCAACCAACgca           >  2:661816/1‑298 (MQ=11)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                              >  2:531147/1‑122 (MQ=40)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                              <  1:531158/122‑1 (MQ=40)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                             >  1:80468/1‑235 (MQ=39)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTc           >  1:432382/1‑300 (MQ=40)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTc           >  1:240361/1‑300 (MQ=40)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTACGTCGGATAAGGCGGTTCCGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGACCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTACGTTGGATAAGTCTGTTCCCCCGCATCCG‑ACCTCCacactc           >  2:943560/1‑295 (MQ=14)                                                                                                                                                                                                                                                      g‑gCCTACA‑CCGC‑TGAGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                             <  2:80468/235‑1 (MQ=255)                                                                                                                                                                                                                                                            aca‑cCGC‑TGTGAAGTGCTCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                                                                                                                                    <  1:765738/127‑1 (MQ=255)                                                                                                                                                                                                                                                            aca‑cCGC‑TGTGAAGTGCTCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                                                                                                                                    >  2:765709/1‑127 (MQ=255)                                                                                                                                                                                                                                                            aca‑cCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCCCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCTCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTCCGAGcc     <  1:337253/301‑1 (MQ=255)                                                                                                                                                                                                                                                            aca‑cCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATcc                                                                                                                                   >  2:414301/1‑175 (MQ=255)                                                                                                                                                                                                                                                            aca‑cCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATcc                                                                                                                                   <  1:414308/176‑2 (MQ=255)                                                                                                                                                                                                                                                                ccGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTCCTAGCCGGt  >  1:871851/1‑301 (MQ=255)                                                                                                                                                                                                                                                                ccGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTCCGAGCCGGt  >  1:478991/1‑301 (MQ=40)                                                                                                                                                                                                                                                                           gTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCATGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAg                                                                   <  1:148459/227‑1 (MQ=255)                                                                                                                                                                                                                                                                           gTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCAGGCCTACACCGCTGTGAAg                                                                   >  2:148458/1‑227 (MQ=255)                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                               CACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAG‑GCCTACA‑CCGC‑TGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCG‑ACATCCAACGTCCGAGCCGGT  >  NC_000913/2304255‑2304808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 32 ≤ 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.