Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,380	+CG	100%	intergenic (+586/+56)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.→C	100.0%	11.9	73	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (29/34):  ref base (5/1):  total (35/38)
Fisher's exact test strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.66e-01
*	NC_000913	4,296,380	2	.→G	100.0%	11.8	73	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (30/35):  ref base (5/1):  total (35/38)
Fisher's exact test strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.53e-01

TAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTT  >  NC_000913/4296144‑4296666                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                                                                  tAAGGCGTTTACGCCGCATGCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGCTAAGACGTTGATGCCGCATCCGACATCAACACCTTATGCGGCGGTTAACGCGTCTTATCAGGCCTACGCCAGACGCCCAAATAGCCGGATGTAGCGTGAATTTGTAGGTCGGATCAGGCGATA‑‑C‑‑A‑‑C‑CAGCATCCGACATCAATGCCTAGTGGGACGATTGCCGCGTGTTCTACGGCGTATGTTAAcc                                                                                                                                                                                                                                    >  1:739110/1‑300 (MQ=11)  aggtgtttgtggagcgtaacGACATGATCGCCGGATGCGACGTTGGACGCGTCTTAGCAGGCCTACTGCAGCAAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGGCGGAGAAGGCGTTGATGCGGCATCCGCCAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                                                                                    <  2:530762/280‑1 (MQ=11)    ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATACGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTCATGAGGCCTATGTTAACCGTTg                                                                                                                                                                                                                                >  2:257762/1‑301 (MQ=11)                       ccATCCAGGCCCGATGCGACGCTTGACGCAGCTTATGTGTCCGACGCAAGCCCACGCAAGGGCCTGATTTAGCTTAATTGTGGAAGTCGGATAAGGCGTGTATGCCGCATACTACATAAGCGCTTGAAGCGACGCTTAACGCTGCTGATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                                                                                                                 <  1:582756/296‑1 (MQ=11)                            ccgtcctggggcgcagctggctgtgtgtcaTCAGGACTAGGCCAGACAGCGCAAAAGGTTGATGTAGACGGAGTGTGTCGGCAGGATAAGGCGTGTATGCCGCAGCCGACAGCAACGCCTCATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                                                                                                                                                                        <  2:429627/272‑1 (MQ=14)                              cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTCGCGTGATTTTGTCGGTCGGATCAGGCGTTTATGCCGCATCCGCCCTCAACGCCTGATGCTACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACCGCGCCGTCGCCTGATTTATCGTGCTTTTATCGGTCGGATACTGCG‑TT‑‑T‑‑ACGC‑C‑GCCTCCGCCGTCACTGCCTGATGCGCCCCTTCCCTCATCTTCTCCGGCCTATCTTCACCTTTGGTTAATTTTTCAAACTCTTccctcc                                                                                                                                                                                                       >  1:450962/1‑295 (MQ=11)                                  cgaggcgcagctggcggctttttATCAGGCCTACGCCAGAACGCGCAAAGGCCTGGTTTAGCGTGAATGTGTAGGGCGGATAAGGAGTTTATGCCGCATCCGACATCAACGGCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                                                                                                                  <  2:38625/280‑1 (MQ=14)                                     tGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCTTCTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGtcctt                                                                                                                                                                                               >  1:38627/1‑297 (MQ=14)                                         acgctggccgcgTCTTACAAGGCCGACTCAAGACAGCTAAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGTGTTCATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                                                                            <  1:648335/292‑1 (MQ=14)                                          ccctgtcgtcgTTTTCTCAGCTCTCCGCCATAACGCGTAATCGCCTTATTTAGCGTGATTTTGTTGGTCGGATAAGGCGTTTATGCCGCACCCGTCATCAACGCTTGATTCGAAGCTTACCGTGTCTTATCAGCCCTCCGCCCGACAGCTCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATTAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAAcc                                                                                                                                                                                           <  1:356066/291‑1 (MQ=255)                                          cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCACCGCCTGATGCGCAGCTTCACGCGTCTTCTCCTGCCTACGCCCGACAGCCCCATCGCCTTCTTTAGCGTGATTTTGTATTTCGGATAAGGCG‑GT‑‑T‑‑A‑‑CGCATCATACGACATCTATGCCTTCTGCGCCGCTTACCCCGTCTTATCAGTCCTCTCTTGACCGTGGGGTACTTTTTCATCCTCTTCACACCCGTTCTCAACcc                                                                                                                                                                                           >  2:684066/1‑300 (MQ=2)                                          cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                                                                           >  1:614966/1‑300 (MQ=255)                                          cGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCCTCCGACCTCAATGCATGATTCTACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTCGGTTCCTTTTTCCAACTCTTCCTACCCGCTacccaccc                                                                                                                                                                                           >  2:101851/1‑293 (MQ=11)                                                   ttggttATCAGGCCGTCGCAGGACACCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATCAGGCGTTTATGCCCCATCCGACATCAACGACGGACGCGACGCTTAACGCTTCTGATCAGTCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                                  <  1:263273/296‑1 (MQ=17)                                                   gtgctgaacAGTCCTACGCGATACATCGCAGTGACATGGTGGAGCGTGGTTGTGGAGGTCGCATAATGCGTGTCTGCCGCATCCGCCATCTATGCCGGGTGCGTAGCTTAACGCGTCTTATCAGGTCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCG‑TT‑‑T‑‑ACGC‑C‑TCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACACTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                                  <  1:446629/292‑1 (MQ=17)                                                   gggtttaaccggcgcaggCCAGACAGCGAAAAAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCACCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                                  <  1:236663/283‑1 (MQ=17)                                                   gagcTTAGCAGGCCTAAGCCATACAGCAGAATAGCCTGATTTAGCGTGAGTTTGTAGGCAGGATAAGGCGTGTATGCCGCAACCGACATGAACGCCTGATGCGACGCGTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                                  <  1:460341/297‑1 (MQ=17)                                                   cgtctgatgaggctgtcggcagcacgcgcAATAGGCTGATTTAGCGTTATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACTCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                                  <  2:300214/275‑1 (MQ=17)                                                        tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCATGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTAATTTTGTAGGTCGGATAAGGCG‑CT‑‑TGCG‑‑C‑C‑GCATCCGACCTCAATGCCTGATGCGACATTTGCCGCGTCTTATCAGCCCTATCTTACCAGTTGGTTTATCTTTCCAGCCCCTCAACCCCGGTATCAGACCATTCCATTCAGCTc                                                                                                                                                                             >  2:108430/1‑300 (MQ=11)                                                                  aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTCGCGTGATTTTGTAGGTCGGATACGGCG‑TT‑‑TCCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACCCCCGATATCAAACCCTTCCCTACCGCTCCGTTTCA‑Aca                                                                                                                                                                  >  1:565694/1‑299 (MQ=14)                                                                  aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATACGGCG‑TT‑‑T‑‑ACGC‑C‑CCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTCTCAGGCCTATCTTACCCGTTGGTTACTTTTTCACACTCTTCACACCCGGTACCCAACCCTTCCATCCACTCCCGGTTCA‑ACCa                                                                                                                                                                 >  2:102714/1‑301 (MQ=14)                                                                      acggcAGCGAAATCGCCGTGTTTTGGGTGATTTGGGAGGGCGGATAAGGCGTTTAGGCCGAATCCGAAATCAACGCCGGAGGGGACGCTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGc                                                                                                                                                             <  2:335581/297‑1 (MQ=18)                                                                            gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTATGTCGGGTACGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGCTGCGCGGCTTGCCGCGTCTTGTCAGGCATATCTTAACCGTTGGTTAATTTTTCCACCTCTTCACACCCGGTATCAAACACTTCCACACAGCACCGGTCCA‑ACCCATGCAccga                                                                                                                                                        >  2:694041/1‑299 (MQ=14)                                                                             cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGACAAACAGCGCAATAGCCTGAATTAGCGTGATTTTTTAGGTCGAATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCAGCCGACCTCACTGCCTGCTGAGACGATTGCCGCGTCTTATCAGGCCTCTATTCACCGTTGTGTACTTTTTCCCACTCTTCACACCCGATCTACAACCCTGCCAtcacgaccaccttaccccccttccccgcc                                                                                                                                                        >  2:492174/1‑272 (MQ=11)                                                                                  tggcctggtttgacggtaTTTTGTAGGTCGGATAAGTGGGTGAGGCCGCATCCGACACCAACGCCTTATGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGc                                                                                                                                                 <  1:257768/284‑1 (MQ=18)                                                                                  tacctcgatgtatcttGAGTTTGTAGTTCGGAGAAGGCTTTTCTTCCCCAGCCGACTGCAACGCCGTATGCTACGCTTAACCCGTCTTATCAGGCTGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACAGCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGc                                                                                                                                                 <  1:213936/286‑1 (MQ=18)                                                                                    gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTAGCGTGATGTTGTAGGTCGGGTAAAGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCGGGCCTATCTTAACCGTTGGTTAATTTGTCAAACTCTTCACACCCGGTGTcca                                                                                                                                                                                              >  1:190918/1‑253 (MQ=14)                                                                                     ccTGATTGAGTGTGAGTTTGGAGGTGTGGTAAGGCGTTTAGTCCGCATCCGACATCAAAGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCt                                                                                                                                                <  2:126663/299‑1 (MQ=18)                                                                                        taTTTAGCGTGATTTTGGAGGTCGGAGAATGCGTTTATTCCTCATCCGACATCAACTCCTGATGCTGCGCCTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGTCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGttc                                                                                                                                           <  1:419299/300‑1 (MQ=18)                                                                                                     tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCA‑ACCAGTGCAGCCCCTTTTGCTTGTTCTTCTCGATaaaaa                                                                                                                              >  2:306614/1‑300 (MQ=255)                                                                                                      tggttgggctgATAAGGCGTTTAAGCCGCATCCGACATCACAGCCTGATGCGCCGCTTAACTGGTTGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                                                                                                                             <  2:342325/291‑1 (MQ=21)                                                                                                                gATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGTCCTACGCCAGACCGCCCAATAGCCTGATTTAGCGTGATTTTGTAGTGCGGATACGGCG‑TTTCT‑‑G‑‑C‑C‑GCATCCGACATCAATGCCTGATTCGCCGCTTGCCGCGTCTTATCACGCCTATCTTAACCGTTGGTTACTGTTTCACACTCTTCACACCCGGTCTCAACTCCTCCCATACACCAC‑TGTTCACACCTGTGCAGCGCCTTTTCCTTGTTCTTCTCTATAAAACCTTCTTCTCCg                                                                                                                   >  2:549795/1‑301 (MQ=17)                                                                                                                      gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACCGCTCAGGTTCA‑ACCATTTCACCGCCTTTTGCTTTTTCTTCTCGATAAAACCTTCTCCACCTTTTaaa                                                                                                             >  1:587919/1‑300 (MQ=18)                                                                                                                       cgttcaagcccgtgcgtcaatacaCGCCTGATGCGACGCTTAACGCGTTTTATGTGGCCTAAGACAGACAGCGCAATAGCCTGTTTTAGCGTGCTTGTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                                                                                             <  1:733439/277‑1 (MQ=25)                                                                                                                         tgtATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑G‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAc                                                                                                          <  1:574436/299‑1 (MQ=255)                                                                                                                         gttacgcccactcagACATCAACGCCTGAAGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCAGCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTGCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                                                                           <  2:565671/286‑1 (MQ=25)                                                                                                                            aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                                                                                                                               <  2:506508/230‑1 (MQ=255)                                                                                                                            aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                                                                                                                               >  1:506531/1‑230 (MQ=255)                                                                                                                               ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATGGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCG‑TT‑‑TTCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                                                                     <  2:184061/300‑1 (MQ=255)                                                                                                                                  cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCACTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTcc                                                                                                  >  1:446622/1‑300 (MQ=255)                                                                                                                                       gtttgaatagacagATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                                                                             <  2:472321/287‑1 (MQ=255)                                                                                                                                                cctggaggtgACGCTGAACGCGTCTTAGTAGGCCGAAGCGAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCT‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                                                                                   <  2:61118/292‑1 (MQ=37)                                                                                                                                                cccgtaTGCGACGCTTAACGCGTCTTATCAGGCCGCCGACAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                                                                                   <  1:102718/296‑1 (MQ=255)                                                                                                                                                 ggggaggaggggctgaaggcggattATGAGGGCGACGCAAGAAAGGGGAATAGCCGGATTTAGCGGGAGTTGGGAGGGCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATGCGACATGAATGCCGGATGCGACGCTTGCCGCGGCTTATCAGGCCGATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                                                                                   <  2:566462/277‑1 (MQ=255)                                                                                                                                                 ccTGATGCGACGGTTAACGCGTCTTATCAGGCTGACGACAGACAGCGCAATAGCCTGAGTTAGCGTGATTTGGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                                                                                   <  1:587205/300‑1 (MQ=255)                                                                                                                                                     aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGCACATCATTCCCGcccctgttctttataca                                                                               >  1:472335/1‑286 (MQ=255)                                                                                                                                                      tggaaggattAACTGGTCTGATCAAGCATAAGACAGACAGCGAAATAGCCTGATTGGTCGTGATTGTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑AAGC‑A‑GAAACCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCAATACAGCGAAGGTTCA‑ACCAGTGCAGCGCCTTTTGCGTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc                                                                              <  1:333430/292‑1 (MQ=255)                                                                                                                                                                cattgtttttactagtgcttcaccAGACAGCGCAATCGCCTGATTTAGCGTGATTTGTGAGGTCGGAGAAGGCC‑TT‑‑T‑‑ACTC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGaa                                                                   <  2:446608/279‑1 (MQ=37)                                                                                                                                                                    cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTCAAACCTTCTTCACCGTTTAAGCACATCCTTCCCTCCCCGTCCTAGGAATAACCGGTTCGGGAAAttt                                                               >  2:231547/1‑300 (MQ=25)                                                                                                                                                                    cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTGGGTTAATTTTTCACACTGTTCAAACACGGTATCAAACCCTTCCATAAAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACGGTTTAAGAAAATCATTCACGCCACGTACTAGGAATCAACGGTTCGGGAAAtga                                                               >  2:147163/1‑299 (MQ=25)                                                                                                                                                                      ttttATCAGGCCCACGCAAGAAAGCGAAATGGCCTGATTTAGCGTGATTTTTGAGTTAGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCATGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGc                                                              <  2:375508/298‑1 (MQ=37)                                                                                                                                                                          tgaaGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGc                                                         <  1:306620/298‑1 (MQ=255)                                                                                                                                                                              ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTACGAAAATCATTCCTGCCCAGTATCCGGCATAAACGGTTAGGGAAATTggcggcctcag                                                     >  1:214888/1‑291 (MQ=25)                                                                                                                                                                                   aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAGTACCTGATGCGGCCCTTGCCGCGTCTTATCAGGCCAATATTAACCGTTGGTTAATTGTTCAAACTCTTCACACCCCGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTCCAGCGCCTTTTTCTTGTTCTTCTCGATACAACCTTCTTCAACGTTTAACAACCTCATTCCCGCGCACTCCTCGCACTCACCGGTTCgtgcacttgaaatgctatttttga                                                >  2:437901/1‑278 (MQ=21)                                                                                                                                                                                   aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCCAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATACAACCTTCTTCACCGTTTCAGAACATCATTCCCGCCCATTACTCGGAATCACAGGTTCGGGAAATTTCAGAGCTTCGtttca                                                >  2:618563/1‑299 (MQ=37)                                                                                                                                                                                   aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCCCCTTTTGCTTGTTCTTCGCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCATCCCAGTACTCGGAAATACCGGTTGGGGAAATTGACGACGTGCGTTTGa                                                >  1:347233/1‑301 (MQ=37)                                                                                                                                                                                   aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGATATCAAACCCTTCCATACAGCTCAGGTGCA‑ACCAGTGAAGCGCCTTTTGCTTGTTCTTCTCGATAAAAACTTATTCACCGTTTAAGAACATCATTCCGCCCCAGTAGTCGGAATCACCGGTATCGGAAATTGCAGAACTTCGTTTGa                                                >  1:599916/1‑301 (MQ=25)                                                                                                                                                                                   aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATATTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCA‑ACCCGTGCAGCGCCTTTTGCTTGTTCTTAGCGATCAAACCTTCTTCACCGTTTAAGAACATCATTACCGCCCAGTACTCAGAATACACTGGTCGGAAAAtgtgagacacgcgtttaa                                                >  1:617287/1‑284 (MQ=25)                                                                                                                                                                                         aaCAGCGCAATAGTCTGATTAAGCGTGAGTTTGAAGGACGGAGAAGGCG‑TT‑‑T‑‑ACTC‑C‑GCATCCGACATCAAGGCCTGATGCGACGCGTGCCGCGTCGTATCAGGCCTATCGTAACCGTTGGTTAATTTGTCAAACTCCTCACACCCGGTAGCAAACCCGACCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTGTTCTCGATAAAACCGTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATAcc                                          <  1:108434/300‑1 (MQ=255)                                                                                                                                                                                                  aTAGCCTGACTGATCGTGATTGTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGATATCAAGGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTTGTTAATTCTTCCAACTCTTCACACCCGGTATCAAACCCCTCCAGACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAACCCTTCTTCACCGTTTAAGAACACCATGCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCt                                 <  1:549818/301‑1 (MQ=38)                                                                                                                                                                                                       tttaTTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATAAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCa                            <  2:614942/298‑1 (MQ=255)                                                                                                                                                                                                         gATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTAACCGTTTGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCAGGTTCA‑ACCAGTGCAGCCCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCTGTTCGGGAAATTGCAGCGCTGCGTTTGACTTCCCAGCTTGCCTTGTCAt                           >  1:201949/1‑300 (MQ=38)                                                                                                                                                                                                         gATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGATCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGTTTGCCTTGTctt                           >  1:252110/1‑298 (MQ=255)                                                                                                                                                                                                                 gggATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCAGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCATTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGcc                   <  2:587894/298‑1 (MQ=255)                                                                                                                                                                                                                  tGATTTTGTAGGTCGTATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCGTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCtt                 <  2:701042/301‑1 (MQ=255)                                                                                                                                                                                                                    aTTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ATGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGCACATCATTCCCGCCCAGTACTCGGAATAACCGTTTCGGGCACTTGACGACCTGCGTTTGAAATACCAGGTTCCCTTTTCATCCTCGGCCTTaa               >  2:744487/1‑301 (MQ=38)                                                                                                                                                                                                                    aTTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑CCCAGTGCAGCGCCTTTTGCTTGTTATTCTCGATAAAAACTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTTCTCGGAATAACCGGTTCGGGAAATTGTAGCGATGAGTTTGAAATACCAGGTTGCATTGTCATCGTCGGCCTTaa               >  1:511248/1‑301 (MQ=38)                                                                                                                                                                                                                      tttGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGTGTCTTATCAGGCCTATCTTAACCGTTGGTTCATTTTTCAAACTCTTCACACCCGGTATCAAAACCTTACATACAGCTCAGGATCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTGTATATAAAACCTTCTTCACCGTTTAAAAACATCATTCCCTCCCAGTCATCGGAATAACCGGTTCGGGAAATTGCAAAACTGCCTTTGCAATACCAGATTCCCTGTACATCTTCGCCCTTCAt              >  1:697617/1‑300 (MQ=37)                                                                                                                                                                                                                      tttGTAGGTCGGATAAGGCG‑TT‑‑T‑‑AAGC‑C‑GCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATg             <  2:347226/301‑1 (MQ=255)                                                                                                                                                                                                                      ttgggaGGTCGGAGAAGGCG‑TT‑‑T‑‑ATGC‑C‑GCAACCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTGATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGGATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGGGCAGCGCCTGTTGCTTGTTCTTCGCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCTTTAATg             <  1:744514/296‑1 (MQ=255)                                                                                                                                                                                                                                ggATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCAAGTACTCGGAATAACCGGTTCGGGAAATTGCAGATCTGCGTTTGAATTACCAGGTTGCCTTGTCATCGTCGTCCTTACTGccacctccat   >  2:210300/1‑293 (MQ=39)                                                                                                                                                                                                                                 gATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCCCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAAACCCGGTATCAAACCCTTCCCTACAGCTCAGGTCCA‑ACCAGTGCCGCGCCTTTTGCTTGTTCTCCTCGCTAAAACCTTCTTACCCGTTTCCGAACATCATTCCCTCCCACTATTCGGACTACCACGTTTGGGAACTTGCAGAGCTGCGTTTGACATACCAGGTTCCCtttactccctcgcgcttatttcccacccagtt  >  1:334762/1‑271 (MQ=255)                                                                                                                                                                                                                                 gATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTACGAACATCATTCCCGCCCAGTACTCGGCATAACCGGTTCGGGCAATTGCAGAGCTGCTTTTGAAATACAAGCTTGCCTTTTCTTCGTCGGCCTTAATGCCCACGCCGtt  >  2:523157/1‑301 (MQ=255)                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                                                                  TAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTT  >  NC_000913/4296144‑4296666

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