Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,191	A→C	22.3%	intergenic (+397/+245)	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,191	0	A→C	22.3%	27.3 ‑2.5 ‑0.0	88	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (15/7):  ref base (29/36):  total (45/43)
Fisher's exact test strand distribution p-value = 8.35e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01

ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT  >  NC_000913/4295893‑4296476                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                  tttgttggtgtgtttTTAGGTGTGAAAGGGTGTTTAGGCCGCAGCCGACATCAACGCCTGTGGCGGCGAGTAACGGGTCTTAGTGGGACGACGCCAGACAGCGCAATAGCCTGATTGAGCGTGATTTTGTAGGTCGGATAAGGCGTTGACGCCGCACCCGACATCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                 <  1:1186012/289‑1 (MQ=11) aTTTAGAGTGAGTTGGTGGTCCGGGTAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTGAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                 <  2:851944/301‑1 (MQ=255)   ttagctggatttgttaggtccgtgagtgactttATGCCGCATCCGACATCACCGCCTTATGCTCTGCTTACCGCGTAGTAGCAGTCCTCCGCCTGACAGTGCAATAGCCTGATTTATCTTGATTTTTTAGGTCTGATAATCCGTTTACGCCGCTTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgcg                                                                                                                                                                                                                                                                                               <  2:1520814/271‑1 (MQ=11)          tgtttgtgaGGTAGGAGAAGGGGGTAATGCCGCATCGGACAAGAAGGCATGATGTGACGTGTAGCGCGGTTTATCAGGCCTAACACAGACAGCGCAATAGCCTGCTTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTGACGCCGCATCGGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                        <  1:1357302/293‑1 (MQ=11)          gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACGCCTCATGCGACGCTTCCCGCGTCTTATCAGGCCTACGCCAGACCGCGCCCTACCCCGCTTTAGCCTCATTTTCTATCTCGGATCAGCCGTTTATGCCCCCTCCGCCTTCAACGCCTCCGCCGACGCTTCACACGCCtcct                                                                                                                                                                                                                                                                                         >  1:1290758/1‑297 (MQ=11)          aaTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAGAAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                        <  2:1417070/300‑1 (MQ=25)           atgttgtaggtgggtgtgggcgtcgagtccGCATACGAGATCAACGCCTGATGCGCAGCTTGACGCGTCTTATCAGGCCTACGCCAGAACGCGCTATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                        <  2:947382/272‑1 (MQ=11)           aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGACTCGTACCGCTTCTTATc                                                                                                                                                                                                                                                                                        >  2:1426678/1‑300 (MQ=255)           aTTTTGTACGGCGGATATGTCGTTTATTCCCATTCAGTCATCACCGCCTGATGCGACGCTTAACGCGTCTTAGCAGGCCTACGTCAGAAAGCGAAATAGCCTGATGGAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                        <  2:523818/300‑1 (MQ=11)                     tcgattAAGGCGGTTCTGCCGCATCCCACATCAACGGCTGATGTGACCCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGcc                                                                                                                                                                                                                                                                             <  2:1384371/296‑1 (MQ=11)                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCCACGCCTGATGCGACCCTTAACGCGTCTTATCCGGCCTACGCCAGa                                                                                                                                                                                                                                                                          >  2:747448/1‑301 (MQ=255)                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCTGCCTCCGACATCAACTCCTGATGCGCCCCTTAACGCATCTTATCacgacgcacccag                                                                                                                                                                                                                                                                           >  1:1529517/1‑288 (MQ=255)                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAAGTCGGATAAGGCGTTTCCGCCGCATCCGACATCAACGCCTGATGCAACGCTTGCCGCGTCTTATCCGGCCTCGGCCGGACAGCGCCCTAGCATGACTTCGCCTAACTTTGGAGGTCGGAACAGGCATTTTAGCCGAATCCGGCATCAAGCACTGATCCAGCGCTTAGCCCGTTGTAACAGGGCTACCCAAGa                                                                                                                                                                                                                                                                          >  1:1419706/1‑301 (MQ=11)                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCCACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATACGGCGTTTACGCCGCACACGACCGCAAGGCCTGATGCGACGCTGGACGCGGCTTATCAGGCCAACGCCAAGAAGCGCAATAGCATGATTTAGCGTGAGTTTGTAGGGTGGAAAAAGACTTTAAGACGCAACAGGAAACGAAGCCTGCTGGCACACATAACGGGTTGTATCAGGCCTaagaacga                                                                                                                                                                                                                                                                          >  2:700030/1‑294 (MQ=2)                         aTAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                         <  2:206437/301‑1 (MQ=25)                         aTAAGGCGTTTATGCCGCACCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                         <  2:218166/301‑1 (MQ=25)                         aTAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGAGTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                         <  2:192568/301‑1 (MQ=21)                         aTAAGGCGTTGACGCCGCATGCGACATCAACGACTTATGCGACGCTTATCGCTGCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                         <  2:93549/301‑1 (MQ=255)                           aaGTCGTTTATGCCGCATCGGACATCAACGCCTGAGGCGAGGCTGAGCGCGTCGTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                         <  2:1290806/299‑1 (MQ=255)                               cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTATGCCGAATCCGACCTCAACGCCTGATGCGACGCTTCACGCGTCTTATCACGCCTACGCCAGACAgccc                                                                                                                                                                                                                                                                    >  2:1127132/1‑298 (MQ=255)                                 tgtagtccGCATCCTACATCACCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                 <  2:23446/295‑1 (MQ=21)                                   tatgcgtgagccgtaaTCAACGCCGGATGGCACGCTTAACGCTGTTTAGCAGCATGAGGAAAGACAGCGCATCAGCCTGATTTAGCGGGGTTTTGTAGGTCGGATAAGGCGCATACGCCGACGCCGGAATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                               <  2:827906/286‑1 (MQ=11)                                   tATTCCGCAACCGACATACACGCCGGATCGGCAGGTTAACTGGTCTTATCAGGTCTACGCCATACAGCGAAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                               <  2:1279049/301‑1 (MQ=11)                                   tATGCCGGAGCGGACAGAAACGCGTGAGGGGACGCGTAAGGCGGTGTAGTAGGCGTAAGTCAGACAGCGAAATGGCCTGATGTAGGGTGATTTGGTGGGGTGGTTACGGCGTTTACGCCGCATCGGATATCAACGCCTGATGCGACGCTTGCGGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                               <  1:1221107/301‑1 (MQ=11)                                   tATGCACCATTCGACACCAACTCCGTAGTCGATGCACAACGTGTTTTATCAGGCCTACGCCAGGCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTACGCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                               <  2:1210299/301‑1 (MQ=11)                                   tATCCCGCACCCTAACTCCACGCCTGGTGACACTCTTACCGCGTTTTATCAGGCCGACGCCAGACAGTGCCAGAGACGGCTGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCCACATCAACGGCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                               <  2:1529576/301‑1 (MQ=11)                                       gcgtgtctgatgctccggcctgggtctgcGCTTACCGGGGCTTTGCACGCATAGGTCAGAGTGCGCAATAGCCTGATTGAGCGTGCTTTGGTAGGTCGGCGAAGGCGGTTACGCCGCATCCGAGATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                                                                                                                                                                                            <  1:3539/272‑1 (MQ=11)                                       cccacctcgaacgcaacctcTGATGCCACTCTGAACGCGTCTTGTCAGGCGTAAGCAAGACAGCGCAATATCCTGATTTTGCGTTATTTTGGTGCTCGGATAAGGCGTTTACGCCGAATTCGCCATCACCGCCTGATGCGACGCTTGCCGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCCACGCCTGATTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                             <  2:1419760/280‑1 (MQ=11)                                               tgcatcgacggcgtaTGCGACGGTTAACGCGTGTGAGCAGGTGTAGGCCAGACAGCGCAATAGCGTGATTTAGGGTGATTTTGTAGGTCGGATACTGCGGTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGGCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                   <  1:1127089/287‑1 (MQ=11)                                               aaCATCAACGCCTGCGGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                   <  2:1086537/300‑1 (MQ=21)                                                agattaacgccggggtgctggCTGAACGAGTCTGATCAGGCTTAGGCACGACAGCGCAATAGCCTGATTAAGCGTGATTTTGTAGGTCGGATAAGGGGTGTACGCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                   <  1:342969/280‑1 (MQ=11)                                                     aaCGCCGGCTGCGACGCGTAACGCGGCTTAACAGGGCGCGGCCAGCAACCCAAATAGCGTAGTTGATCGTGATTTTGTAGGTCGGATAAGGCGTCTACGCCGCATCCTGCATCAACGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAACAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCTACATCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAATAGCCTGATTTAGCGTGAttt                                                                                                                                                                                                                                             <  2:42342/301‑1 (MQ=11)                                                       cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                                                                                                                                                                                           >  2:987768/1‑301 (MQ=34)                                                          tggggtggccgcttgcctgttCTTATCAGGCCGTGGCTGGGCGGTGCAATGGTCTGAATTGGCGGGATTTTGGAGGTCGGAGAGGGTGGTGACGCCGCATCCTGCATAAAGGCCTGATGCGACGCTTGCCGGGTCTTATCAGGCCGACGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                        <  1:1053893/281‑1 (MQ=11)                                                                    gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTCGCGTGATTTTGTAGTTCTGATAATGCGTTTATGCCGCCTCCGACATCAACGCCTGATTCTCCGCTTACCGCCTCTTATCAGGCCTACGCCAGACAGCGCAACACCCTgctctcccgtgtgttgtctgtgcggcgtag                                                                                                                                                                                                                              >  1:1094641/1‑272 (MQ=11)                                                                       aaactgttcgtttcagtcctagagtaGAGAGAGCAATGACCTGATGTAGTGGGATTTTGTAGGTCGGGTAAGGCGTTTGCGCCGCATAAGACATGAAAGCCGGATGCGACGCTTGGCGCGTCTTATCAAGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGACAGGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCg                                                                                                                                                                                                                           <  2:1042184/276‑1 (MQ=11)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAAGCTGCACGCGTCTTATCAGGCCTACGCCAAACCGCGCAAGAGCATGATTTAGCGTGCTTTTGTATGTCGCAGAAGACGTTTagccac                                                                                                                                                                                                                  >  1:503409/1‑295 (MQ=11)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCCAACCGCGCAATAGCCTTATTTAGCTTGATTTTTTTAGTCGGATAAGGCGTTTacgacg                                                                                                                                                                                                                  >  1:1030228/1‑295 (MQ=11)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGACTTATCGTGATTTTTTAGGTCGGATAAGGCGTTTACGccc                                                                                                                                                                                                                  >  2:1306970/1‑299 (MQ=12)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATCCGACGCTTACCGCGTCTTATc                                                                                                                                                                                                                                                                                        >  1:1301321/1‑230 (MQ=18)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                        <  2:1301370/230‑1 (MQ=33)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCGTTATTTTGTAGGTCGAATAAGTCGTTTACGCCg                                                                                                                                                                                                                  >  2:646367/1‑300 (MQ=18)                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGCCAGCGCAATGACCTCATTTAGCGTGATTTTTGAGGTCGACTAAGCCGTTTCGGCCGCATCCg                                                                                                                                                                                                            >  2:1119746/1‑301 (MQ=11)                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATAAGGCCTACGCCCGACAGCGCACTAGCCCGATTTCGCGTGCTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                                                                            >  2:1498192/1‑301 (MQ=255)                                                                                        cctctgccaaaccgctcataaGCCTGTTTGAGTGTTGTTTGTTAGGTCTGATAAGGATTTTATCCCGCATCCTACATCAACGCATGCGGCGCCGCGTACCGCGTCTTATCTGCCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                                                                                          <  2:1378661/281‑1 (MQ=11)                                                                                        ccTCCGCCAGACAGCGCAATTGCAGGATTTAGCGTGATCTTGTAGGTCGGATAAGGCGGTTACTCCGCACCCGACATCAACGCGTGATGCGACGCTTGCGGCGGCTTATCAGGCCTACGCCATACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                                                                                          <  2:503418/301‑1 (MQ=11)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCCCCAGACAGCGCCCTAGACTGATTTAGCGTGATTTTGTAGACCGGATAAGGCGTTCATGCCGCCTCCGACATCAACCCCTTATCCGACGCTTAACTCGTCTTCTACGGCCTACCCCAAACATCCCCATAACCTGATGTAGCTTGATTTTGTTGGTctgcaaagggttttcagcgccctcagacac                                                                                                                                                                                                        >  1:1402541/1‑271 (MQ=11)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGCATAAGGCGTTTACGCCGCATCCGAcgt                                                                                                                                                                                                        >  2:660530/1‑298 (MQ=255)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGCCGCATCCGACAta                                                                                                                                                                                                       >  2:1467096/1‑300 (MQ=25)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGCCGAATCCGACAt                                                                                                                                                                                                        >  2:1166888/1‑300 (MQ=25)                                                                                            tgcgggccAGCGCAATGGTCTGATTTAGCGTGATTAGGGAGGGCGGAAAATGTGTTACAGCCGGCTCCAACAACAAGGCGTGAGGGGACGGGTGCCGCGTCTCATCAGGCCGACGCGAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGGCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                                                                                                                                                       <  2:1071841/293‑1 (MQ=11)                                                                                                aaaCAGCGCAAGAGCTCGGTGTAGGGGGATTTTGTAGGTGTGATAAGGCGTGTACGCCGCATCCGAAAGCAACGCCTTATGCGACGCGTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                                                                  <  1:987751/299‑1 (MQ=11)                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGGCCAGACAGCGCAATAGCCTGAGTTCGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCCTCACACCTCAACGCCTGATGCGAGTCTTAACGCGTCTTATCCGGCTTACGCACAACACCGCAATATCCTCATTTAGACTCATTTTGTATGTGTGATCAGGCGtttcccgcccatccaacagccacg                                                                                                                                                                                                   >  1:722330/1‑279 (MQ=11)                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTAGCGTGTTTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGACTCTTACCGCGTCTTCTCAGGCCTCCCCCAGACCTACCCATAACATGATTTAGCCTATTTTTGTAACGCGGAACACGCTTTTCACCCCCATCCGACATCactg                                                                                                                                                                                                   >  1:1369553/1‑297 (MQ=11)                                                                                                     gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGCTTATCAGGCCTACGCCAGACAGGGCAATAGCCTGATTTAGCGTGAGTTTGTCGGTCGGATAAGGCGTTTCCGCCGCATCCCCCATAAATGCCTGag                                                                                                                                                                                             >  1:18821/1‑300 (MQ=255)                                                                                                      ctaaaTTGCTTGATTTAGGGTGATTTTGTCGGTCGGATAATGCTTTTACGCCGCACCCGACACCAACGCCTGATTCGACGCCTGCCGCGTCTTATCAGGCATACGCCAGCCCGCGCAATAGCCTGATTTATCGTGATTGTGTAGGTCGGATAAGGCGTGTAGGCCGCATCCGACAGCAACGCCTGAGGCGACGCTTAACGCGTCTCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                            <  2:357540/298‑1 (MQ=11)                                                                                                      cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACATCAACGCCTGATGCGACACTTAACGAGTCTTATCAGGCATCCTCCAGACAGCGCAATAGCCTACTTTAGCGTGATTTTTTACGTCGGATAAGGCGCTTCTGCCGCATCCGACATCAACTCCTg                                                                                                                                                                                               >  2:1276689/1‑298 (MQ=11)                                                                                                      cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTCATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTGATTTGGTAGGTCGGATAAGGCTTTTACCCCGCATCCCACATCacaccctcgtt                                                                                                                                                                                            >  1:398585/1‑291 (MQ=11)                                                                                                      cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTACGCCAGACACCGCAATATCCTGATTTCGCGTGATTTTGTCGGGCGGATAAGGCGTTGCTGCCGCATCCGACATCAACGCCTACAGAGACGCTTCACGCGTCTTATCATACCTACGCCCGAGCGCGAAAGACCGTGCTCTAGCGCGATGTTGTAGATCGGATGTGGCGTTGATAGAGCATACGACCTCGACACATGaaa                                                                                                                                                                                            >  1:1076959/1‑299 (MQ=2)                                                                                                      cgaaaTAGCGTGATTGAGCGTGATTGTTGAGGTCGGATAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                            <  2:373117/298‑1 (MQ=255)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGCGTATTATCACGCCTCCGCCAGACAGCGCACTATCCTGATCTCGCGTGATTTTGTCGGTCGGATCCGGCTTTTCCGCCGCCTCCGCCCTCACTGCCTGATGCgcct                                                                                                                                                                                       >  1:657087/1‑296 (MQ=11)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACATCGCACTAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTACGCCGCATCCGCCactcatgcctgcacaggcg                                                                                                                                                                                       >  1:275673/1‑281 (MQ=11)                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATAAACGCATGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTATGCAGCATCCGCCATCAACGCCTGATGAGACGCTTAACGAGTCTTATCAGCCCTACGCCAAACCACAACATAGACTGATTTGGCGTACTTTTGTAGATCGAGTAAGGCACTGAGGCCGCACCCGActagaatgccgtctgcgaaccttcc                                                                                                                                                                                  >  2:1570158/1‑276 (MQ=11)                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTGGTATCAGGGGTTTCTGCAGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAACGCAATAGCATGATTTAGCGTGCTTTTGTAGGGCGTGTAAGCCGTCTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTaac                                                                                                                                                                       >  2:1492254/1‑298 (MQ=11)                                                                                                                                       ggtactggcccgttcagccccACCCGACATACACGCTGAATGCACCGCTTGCCGCGTCTTATCAGGCATAGCCAACACAGCTCAATAGCCTGAGTTGGCGTGATTGTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                             <  1:37446/283‑4 (MQ=11)                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                                                                                                                                                                                            <  2:996715/236‑1 (MQ=33)                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                                                                                                                                                                                            >  1:996698/1‑236 (MQ=33)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAATTCCTGATGCGACTCTTGCCGCGTCTTGTCCGGCCTatcttac                                                                                                                                                          >  2:1453122/1‑294 (MQ=11)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTCAGCGTGATTTTGTAGGTTGGATCAGGAGTTTACGCCGCATCCGACCTCAATGCCTGATGCGGGGCTTGCCGCGTCTTAtacggtctctcataa                                                                                                                                                          >  2:1323812/1‑286 (MQ=11)                                                                                                                                              ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCtctct                                                                                                                                                             >  2:37448/1‑288 (MQ=9)                                                                                                                                                             tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                       >  2:253487/1‑267 (MQ=14)                                                                                                                                                                                                                                               tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAgg                                                 >  2:332469/1‑301 (MQ=18)                                                                                                                                                                                                                                                           ataactcggtgacgCCGCACCCGACATAAACCCCTGAGGCGACTCGTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCCGATTGAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGAGGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                     <  2:1369602/288‑1 (MQ=11)                                                                                                                                                                                                                                                             gacaggTTTACGCCCCATCGGACACCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                    <  2:381404/295‑1 (MQ=11)                                                                                                                                                                                                                                                             atgggactaagggacagtctgaactCAACGGCGGATGGGACGATGGCCGTGGCTGATCAGGCTGACGCCAGAGAGCGAAATAGTATGAGTTAGCGGGAGGTTGTAGGGCGGAAAAGGCGTTCAGGCCGCATCCGACATCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCGTAAcc                                     <  1:102157/275‑1 (MQ=11)                                                                                                                                                                                                                                                             aaGGCGTTTACGCCGCATCCGACATCAACGCTTGATTCGACGCGTGCCGCGTCTTATCAGGCCTACGCCATACAGCTGAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                     <  2:275677/299‑1 (MQ=11)                                                                                                                                                                                                                                                               ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTAGGATAAGGCGTTT‑‑ACGCCGCATCCACCATCAATGCCTGCTGCGACCGTTGCCGCGTCTTATCAGGCCCATCTTACCCGTTg                                 >  1:357821/1‑301 (MQ=11)                                                                                                                                                                                                                                                                 cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGTCCTATCTTAACCGTTGGt                               >  2:570851/1‑301 (MQ=255)                                                                                                                                                                                                                                                                            cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATTCCTGATGCGCCGCTTGCAGCGTCCTATCAGGCCTATCTTAACCGTTGGTTAATTTTTcca                    >  2:1311620/1‑299 (MQ=11)                                                                                                                                                                                                                                                                              tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTTTCCGGCCTATCTTACccctgtctcttctacccatc                   >  2:102160/1‑282 (MQ=11)                                                                                                                                                                                                                                                                              accgcaaTCAACGCCTGATGCGCCTCTTGCCGCGTCTTATCAGGCCGAAGCAAGATGGCGAAATAGCCTGATGTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTATGCCGCATCCAACATGAACGCCTGATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                  <  2:446051/295‑1 (MQ=11)                                                                                                                                                                                                                                                                                    aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCCACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGCCTAAGGCGTTTCCG‑‑CCGCCTCCCCCATCATTGCCTGCTGCGCCCCCTTCCTCCTCTTATCAGGCCTATCCTCCCCTTTCGTTCACTTTTACAActccca             >  1:548272/1‑297 (MQ=1)                                                                                                                                                                                                                                                                                     tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTGATCAGGCCTATCCTAACCGTTGGTTAATTTTTCAAACTCTTca            >  2:852433/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                     tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTca            >  2:423885/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        ccGCCGGGTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  2:18822/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                        acgcctgtgggtaCGCTTGCCGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  2:1030272/289‑1 (MQ=14)                                                                                                                                                                                                                                                                                         cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGGTGCGACGATTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATATAGCGTGATTTTGTCGGTCGGATCAGGCGTTT‑‑ACGCCGCATCCGCCATCCATGCCTGCACCGACTCTTTCCGCGTCTTAGACGGCCAAACTTAACAGTAGGTCTATTTTTCACACTCTTCCCAcc        >  1:485007/1‑300 (MQ=11)                                                                                                                                                                                                                                                                                              gATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGAAATCCATGCCTGATGCGACGCTTGCAGCGTGTTATCGTGCCTATCTTAACCGTTGGTTAATTTTTCAAAATCTTCACACCCGGTaa  >  1:545997/1‑300 (MQ=14)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                  ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT  >  NC_000913/4295893‑4296476

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.