Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,191	A→C	14.4%	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	14.4%	18.8 ‑1.5 ‑0.0	86	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (8/5):  ref base (29/42):  total (38/48)
Fisher's exact test strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.97e-01

ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC  >  NC_000913/4295893‑4296481                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                     ttttTGCGTGATTTTGCTGGTCGGATATGTCGTTTATGCGGTAGCGGTCAGCACCTCCGGGTGGTCAGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GAGACGCTTACGGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                    <  1:709092/300‑1 (MQ=11) ttttAGTGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                    <  2:683940/300‑1 (MQ=25)          tattttttcgagtggAAAAGTCGGTCATGCCGCAGCCGCCACCAAGGACGGAGTGGACGCTTGCCGCGTCGTATAAGGCCTACGCCAGAAAGCGCTATAGCCTGATTTAGCGTTATTTTGTAGGTCGGAGAAGGCGTCTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                           <  1:661502/288‑1 (MQ=11)           agtttgaggttCGGATAAGGCATTTAGCAGGAAGCCGGCATAAACGACTGAGGAGACCGTAGACGCGTGTTATCAGGCGTACGCTAGACAGGGCAATAGTCTGATTTAGGGTGATTTTGTAGGGCGGATACGGCGTTTACGCCGCACCCGGCATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                           <  1:1153018/290‑1 (MQ=11)                   ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTAc                                                                                                                                                                                                                                                                                   >  2:624690/1‑300 (MQ=35)                      cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTCGGTC‑GGATAAGGCGTTTATGCCGCATCCGCCATCAACTC‑CTGATGCGACGCTTCACGCGTCTTATCAGGCCTAcaca                                                                                                                                                                                                                                                                                >  1:494678/1‑297 (MQ=255)                       ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCCTCCGAAATCAACAC‑CTGATGCGACGTGTAACGAGTCTTATCAGGCCTACGcccg                                                                                                                                                                                                                                                                              >  2:1112590/1‑299 (MQ=255)                        gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATATCCTTATTTAGCGTGATTTTGTAGGTC‑GGATCAAGCGTTTCTTCCGCATTCGACATCCACCC‑CCGTATAGACGCTTATCGCGTCTTATCAGGCCTACTCCAGa                                                                                                                                                                                                                                                                             >  1:588463/1‑301 (MQ=11)                         aTAAGGTGTTGATGTCCCCTCAGACATTAACGCCTGAGGCGACGCTTAACGCGTCTTATCAGGCCTACGCCTGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:962454/301‑1 (MQ=255)                         aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:89890/301‑1 (MQ=35)                         aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:21337/301‑1 (MQ=35)                         aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:421174/301‑1 (MQ=35)                         aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:759459/301‑1 (MQ=35)                         aTAAGGCGGATAAGCCGACACCTACATCAAAGCTGGCTGCGACGCTTGACGCGGTGTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  1:87237/301‑1 (MQ=11)                         aTAAGCCGTTCACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:811959/301‑1 (MQ=14)                              gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTATCCTGCTTTAGCTTGATTTTGTATGTC‑GGATAAGGCCTTTATGCCGCATCCGAACTCCACGC‑CTGATCTGACGCTTAACGCGTCCTATCAAGCCTCCGCCACAccggt                                                                                                                                                                                                                                                                        >  1:937732/1‑296 (MQ=11)                               gggttggtgcgcaccgagacacAACGTCTGGTGCCTCGCGTAACGCGTCCTATCAGGCCTACGACAGACAGCGCAAGTGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                                       <  1:352829/279‑1 (MQ=11)                                   tATGCCGCCCCCCGCATCAACACCTGATGCCACGTTTACCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAGCGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                                  <  2:191352/301‑1 (MQ=255)                                   tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                                  <  2:1174037/301‑1 (MQ=25)                                     ccggcgactgcagatgtgtataagcaccagcGCTTAGCGCTTCTAACCAGGCCACCGCCAACCGTCGCTAGAGGCTTATTTAGCGTGAGTTTGTAGGACGGGTATGGCGTTTACGCCGCATCCGACATCAACGCCTGTT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                                <  1:1074016/271‑1 (MQ=11)                                       ccgaagcggacacaaatgacggATGCGCCGCTTAACCGGTCGTATCAGTCTTATGCCAGACAGCGCAATAGCCTGAGTTGTCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGAGGCTTGACGCGTCTGATCAGGCCTACGCCAAAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                                  <  1:530371/276‑1 (MQ=11)                                            tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAACCAGCGCAATAGCCTGCTCTCGCCTGATTTTGTCGGTC‑GGATCCGGCGTTCAGGCCGCATCCGACATCAACGT‑CTGATGCGACGCTCACCGCGACTTATCAGGCCAACCcctacaactgcattatcttgatttt                                                                                                                                                                                                                                                         >  1:346145/1‑278 (MQ=11)                                                aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCgg                                                                                                                                                                                                                                                     >  2:86178/1‑300 (MQ=21)                                                  aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTTATGCGACGATTAACGCCTCTTATCAGGCCTACGCCAGACGTCGCAATAGCCTGATTTAACGTg                                                                                                                                                                                                                                                    >  1:1107473/1‑300 (MQ=255)                                                            taagattgttacaaagtagagtaccAGCCCTACGCCAGCAAGCGCAATATCCTGATTTAGCGTGATTTGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCAACGGCTGAT‑GCGACGCTTGCCGCTTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTC‑GGATAAGGCGTTTATGCCGCATCCGGCATCAACGC‑CTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg                                                                                                                                                                                                                                          <  1:954109/276‑1 (MQ=11)                                                            gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATATAGCGTCATTTTGTAGATC‑GGATAAGGCGTTTATGCCGCATCCGACCTC‑ACGCACTGATGCGACGTTTTACGCCTCTTTTCAGGCTTCCTCCCGACTACGCAACAACCttccttagctgtcttatgtatgt                                                                                                                                                                                                                                         >  1:866704/1‑279 (MQ=11)                                                             aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:961161/265‑1 (MQ=35)                                                                   cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CCGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                                >  2:1074054/1‑271 (MQ=25)                                                                    gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:20753/258‑1 (MQ=35)                                                                            cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                                                                                                                                                                                          >  2:905549/1‑300 (MQ=17)                                                                                 tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg                                                                                                                                                                                                                                          >  2:954139/1‑279 (MQ=33)                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTGACGCGTCTTATCAGGGCTTCGCCAGACAGCGCACTAGCCCGATTTTGCCTTGTTTTGTATGTCGGAGAAGGCGCTTACGCCGCATcg                                                                                                                                                                                                                >  1:551941/1‑299 (MQ=11)                                                                                      ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                                  <  2:599593/250‑1 (MQ=33)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTACGCCGCCTCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                           <  2:4331/268‑1 (MQ=18)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                                                                                                                                                                                          <  2:1004725/237‑1 (MQ=33)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGACTGATTTAGCGTGATTTTGTAGTACGGATAAGGCGTTTACGCCGCATCCGACAt                                                                                                                                                                                                           >  1:400126/1‑300 (MQ=255)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTACCGCAGCTTATCAGGCCTACGCAAGACAg                                                                                                                                                                                                                                                                          >  1:1004694/1‑237 (MQ=16)                                                                                           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                                                                                                                                                          <  2:140478/301‑1 (MQ=25)                                                                                            cccaagactggtctagtgtcagttttTGCCTGATTTTGTAGCGGGTAAATGGCGTTCATGCCGCATCCGTCATC‑ACTCCCGTTGGTGACATTTTCTGCGTCTTACAAGACCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGTTTTAGTAGGTC‑AGAAAAGGCGTTTATGCCGCCTCCGACAGCAACGC‑CTGATGCGACGCTTAACGCGTCGTGTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                                                                                                                                                          <  1:1112552/277‑1 (MQ=2)                                                                                              caagcgaGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCTTTGAGCCGGAATCCGACATCAACGGCTGTG‑GCGACGCTTGGCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                                                                                                                                                                                                        <  1:1124602/294‑1 (MQ=11)                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                                                                                 <  2:747166/225‑1 (MQ=33)                                                                                                aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                                                                                 >  1:747149/1‑225 (MQ=33)                                                                                                      ataagggaacgtaattAGAGTTATTTGTTAGGTAGGATAAGGCATTTACGCCGCAACCGAGATCAACACCTGAT‑GAGACGCTTGCCGCGTCTTATCGGGCCTACGCCAGGCAGCGAAAGAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                               <  1:179989/287‑1 (MQ=11)                                                                                                       tccaaggttgataagagacacaTTTTGGAGGTCGGATAATGCGTTTACGCCGCATCCGGCATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGTCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                                                                              <  1:683026/280‑1 (MQ=11)                                                                                                        cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATCGCCTGTATCACCTTGATTTTGTATGCCGGATAAGCCGTTTACCCCGCATTCGACATccagtcctgctt                                                                                                                                                                                               >  1:889534/1‑289 (MQ=11)                                                                                                         aaTGGCCTGATTTGTCGTGATTTTGTAGGGCTGATAAGGCCTTTACGCCGCACCCGACATGAACGCCTGAT‑GCGACGCTTACCGGGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                                                                            <  1:623775/301‑1 (MQ=11)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTCTGTAGGTC‑GGATCAGGCGTTTCTGCCGCTTCCGACATCAACGC‑CTGATTAGCCGCTTCACTCGTCGTATAAGCACTACCCCAGCCATCGCAATATCCTACCTGATCGTGCTTTTGTCGGTCCGATCAGGCGATTACGCCGCCTCCGCCCTCaatcacttgtgcgacg                                                                                                                                                                                          >  1:615825/1‑285 (MQ=11)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAATTCCTGATGCGACg                                                                                                                                                                                          >  1:807103/1‑299 (MQ=18)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:183488/217‑1 (MQ=34)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            >  1:183486/1‑217 (MQ=34)                                                                                                             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGAGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                                                                          >  2:662156/1‑299 (MQ=21)                                                                                                                 gagttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  1:379514/210‑1 (MQ=21)                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACTCCAGACCGCGCAATCGCCTGCTTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCttgc                                                                                                                                                                                     >  1:475142/1‑298 (MQ=255)                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            >  2:379521/1‑213 (MQ=34)                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                               >  2:179991/1‑290 (MQ=21)                                                                                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                                                                                 >  1:23001/1‑202 (MQ=34)                                                                                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                                                                                 <  2:23001/202‑1 (MQ=34)                                                                                                                          tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCATCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                           >  2:727352/1‑301 (MQ=17)                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                                                                              >  2:683038/1‑280 (MQ=21)                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                          <  2:807125/300‑1 (MQ=17)                                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                          <  2:551949/300‑1 (MQ=17)                                                                                                                                 ttagatcggtgtaggactttACGCTCCAACCGACCTCATAGCCGGCT‑GCGAATCTTGCCGCGTCTTATCAGTGCTAAGCCGAAAAGCGCAATAGCCTGATTTGACGTGATTTTGTTGGTC‑GGGTAAGGCGTTTATGCCGCATCCGACATCACCAC‑CCGCTGAGACGCTTAGCGAGTGTTAGAAGGCCAACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                     <  1:630287/283‑1 (MQ=2)                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                                                                                                                                                          >  2:202374/1‑257 (MQ=37)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                    >  2:523850/1‑197 (MQ=34)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                    <  1:523843/197‑1 (MQ=34)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttaa                                                                                                                                                             >  2:354805/1‑294 (MQ=2)                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttaa                                                                                                                                                             >  2:365528/1‑294 (MQ=2)                                                                                                                                          taaaGGCGTTCACGCCGGCACGGAAGTCAACGTCCGTC‑GTGATGCTTGGCGTGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCAGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttaac                                                                                                                                                            <  1:6546/298‑8 (MQ=11)                                                                                                                                          aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATACGGCGTTTTCGCCTCACACGACATCACTGCCTGATTCGACGcgtgccggtgattatccggcctatcttaac                                                                                                                                                            >  1:939291/1‑271 (MQ=11)                                                                                                                                              ggCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                               <  2:923857/261‑1 (MQ=21)                                                                                                                                              ggCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                                                                               >  1:923829/1‑261 (MQ=21)                                                                                                                                                cGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAATGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            >  1:985874/1‑182 (MQ=21)                                                                                                                                                cGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                            <  2:985905/182‑1 (MQ=34)                                                                                                                                                      cgccgcATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGTCCTACGCCAGACAgc                                                                                                                                                                                                                                                                         >  1:2721/1‑179 (MQ=18)                                                                                                                                                      cgccgcATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgc                                                                                                                                                                                                                                                                         <  2:2721/179‑1 (MQ=33)                                                                                                                                                            acccGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                         <  2:1016393/299‑1 (MQ=21)                                                                                                                                                                                                               cctggACAGTGGAATGGGCTGTGTTAGCGTGATTTTGTATGTC‑GGATAAGGAGTGTATGCCGCATCCGCGATCAACGG‑CTGGTGAGACGCTTAACTCGTCTTATAAGGCCTAAGTCATATAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                                                                                       <  1:505698/297‑1 (MQ=11)                                                                                                                                                                                                                                                             atacagcagtgccgacggtccagACAGCAAGGC‑ATGATGCGACGCTGTCCGGGTCTTATCAGGCTTACGCCAGGCAGCGAAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCAGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTATCTTAAc                                           <  1:86177/278‑1 (MQ=1)                                                                                                                                                                                                                                                             aTAAGGCGTTTACGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                           <  2:939319/300‑1 (MQ=2)                                                                                                                                                                                                                                                              aaaGGCGTTTACGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                         <  2:494685/300‑1 (MQ=2)                                                                                                                                                                                                                                                              aaaGGCGTTTACGCCGAATCCGACATCAACGC‑CGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCAACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                         <  2:838848/300‑1 (MQ=2)                                                                                                                                                                                                                                                                              cATCCGACATCAACGC‑CTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAAACAGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCTTTCGCCATCAATTCCTTATGCGCCGCTTGCCTGTTATTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                         >  1:7590/1‑301 (MQ=11)                                                                                                                                                                                                                                                                                tCCGACATCAATGC‑CTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCcaa                        >  2:8877/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                      aTCAACGC‑CTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGCCATCCACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACCCCAGACAGCGCAATAGCCTGCTTTCGCGTGCTTTTGTCGGTCGGACACGGCGTTAACGCCGCATCCGACATCGATGCGTGATGCGACGCTTGCCGCGTCTTATCGGGCCTATCTTAACCGTTGATTAATTTTTCAAACTCTTca                 >  1:827362/1‑301 (MQ=11)                                                                                                                                                                                                                                                                                                       aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc  <  2:7590/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc  <  2:915699/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                     ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTC‑GGATAAGGCGTTTATGCCGCATCCGACATCAACGC‑CTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC  >  NC_000913/4295893‑4296481

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