Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,190	A→G	16.3%	intergenic (+396/+246)	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,190	0	A→G	16.3%	18.8 ‑1.1 ‑0.1	69	intergenic (+396/+246)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (5/6):  ref base (34/21):  total (41/28)
Fisher's exact test strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.04e-01

GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT  >  NC_000913/4295892‑4296476                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                 gATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAGGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGCCCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTCACGACACGGCGTTTCTGCCGCGTCGGACTTCAACGACTGCTGCGACGTTTaa                                                                                                                                                                                                                                                                                                 >  2:619256/1‑300 (MQ=11)            aTTTTGTAGTGCGTATAAGGGCTTGCTGCCGCAGCCGCCATCAACGCCTGATGCGACGCGTAACGGTTATTATAAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                      <  2:659173/300‑1 (MQ=11)            aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACCCCTGATGCTACGCTTAACGCCTCTTAt                                                                                                                                                                                                                                                                                       >  2:508323/1‑299 (MQ=255)            aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGCCGCTTAAcactacttctc                                                                                                                                                                                                                                                                                      >  2:437338/1‑290 (MQ=255)            aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCACCGCCTGATGCGACGCTTACCGAGTCTTaac                                                                                                                                                                                                                                                                                      >  2:455412/1‑298 (MQ=255)                tttAGGTCGGATAAGGCTTTTATGCCGCAACCTACACCAACTCTTGATGCTACGCTTAACGCGTGTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCAGCACCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                 <  2:968503/299‑1 (MQ=11)                ttatgtgaggATAAGGCGTTTATGCCGCATCCGAAATCACCTGCTGGTGGGGTGTTTACGGCGTCTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGTTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                 <  1:250521/293‑1 (MQ=11)                tGTAGGTCTGTTAAGGCGTTTATGCCCCATCCGACATCAACGCCTGATGCGACGCGTAACGCGGCTGATCAGGCCTACGCCAGACAGCGAAATAGCCTGATGTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTGGACGCGTCTGAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                 <  2:532980/301‑1 (MQ=11)                 ggtggTCGGAAAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGTTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                 <  2:751744/297‑1 (MQ=255)                         gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATAAGTCCTTTATGCCTACTCCGACATCAACGCCTGATCCCAGGCTTCACGCGCATTATCAGGCCTACCCCaa                                                                                                                                                                                                                                                                         >  1:1093825/1‑299 (MQ=11)                         gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTGGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTATAGGTCGGGAAAGGGCTTTTTGCGCGATCCTGAACCAACGCCTGTGGCTCCCCTTAAGCGGTCTTATCGGGCCCACCCCAGa                                                                                                                                                                                                                                                                        >  1:440153/1‑301 (MQ=11)                          aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                       <  2:444374/301‑1 (MQ=35)                           gaaGGCGTTTAAGCCGCATCCGGCATCAACGCCTGATGCGACGCTTAACGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                       <  2:1107302/299‑1 (MQ=25)                                    ttgccacgagcagaaattcaggcgggATGCGTCGTGTAGCGGGTCTTAGCAGCGCAAAGGCGGACAGGGCGAGTGCGTGAGTGAGCGTGATTTTGGAGGTGGGATAAGGCGTTTACGCCGCACCCTACAGCAACGCCTGATGCGACGGTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                             <  1:825387/276‑1 (MQ=11)                                      tGCCGCATCCGACATCAAATCCTGAGGCGACGCTTAACGCGTCGGATCAGGCAGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCAACGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                           <  2:79027/301‑1 (MQ=11)                                            cgccaaCATCAACGCCTGATGCCACTCTTACTGCGTCTTATCAGGCCTAGGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGGATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                                                                                                                                     <  2:243997/296‑1 (MQ=255)                                                      cccgcgttaTGCGTCGCTTAAGCTTTCTTATCAGGTGTACGGCAGACAGCGCAATAGCCTTATTGGTCGTGATTTTGGAGGTCGGATAAGGCGTTTACGCCGCATCCGTCATGAACGCCGGAGGGGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAttt                                                                                                                                                                                                                                           <  2:600907/293‑1 (MQ=11)                                                        cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACCGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGAATCAGGCGTTTCTGCCGCATCCTCCATCAACGCCTGATGCGCCGCCTCACGCGTCTTAGCAGCCCTACACCCACCCACGCATTCCCCTGCTTCAACGTGATTTTg                                                                                                                                                                                                                                         >  1:1230591/1‑301 (MQ=11)                                                        cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACGTCGATAAAGCCTGATAAAGGGTGATATTg                                                                                                                                                                                                                                         >  2:171210/1‑301 (MQ=11)                                                          ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCCTTTATGCCGCATTCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGCCCTCCGCCAGCCACCGCacacccatactttaccgtaattttgtt                                                                                                                                                                                                                                       >  1:68243/1‑275 (MQ=11)                                                          ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTa                                                                                                                                                                                                                                       >  2:102124/1‑301 (MQ=21)                                                           ctgtggCGACGCTTAACGCGTCTTGTCAGGCCTAGGCCAAACAGCGCATTAGCCTGATTTAGCGTGATTTGTTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                      <  2:440163/296‑1 (MQ=255)                                                            tGATGCGACGGTTAACGCGTCTTAGCAGGCCGACGCCAGACAGCGCAATGGCCTGATTTAGGGTGATTTTGTAGGTCGGTTGAGGCCTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                      <  1:544405/300‑1 (MQ=255)                                                              aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                       <  2:1244016/265‑1 (MQ=25)                                                                     gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCCTGATTTTGTAGGTCGGATAAGCCGTTTATGCCCAATCCCACCTCCACCCCTGATGCGACGCTTAACGACTCGTCTCAGGCCTCACCCAGAACGCGCATCTGCCCGATGTCGGGTAACTTTGTGCGTCGGACaaa                                                                                                                                                                                                                            >  1:1076978/1‑300 (MQ=11)                                                                                cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGc                                                                                                                                                                                                                  >  2:295277/1‑300 (MQ=17)                                                                                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCCAGACCGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGGGCTTTACGCCg                                                                                                                                                                                                                >  2:1216211/1‑300 (MQ=255)                                                                                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGTCCTCCGCCAGCCCGCGCACTAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                >  2:1216266/1‑300 (MQ=255)                                                                                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGCCATCGCACTAGCCTGATTTCGCGTGCTTTTGTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                >  2:590908/1‑300 (MQ=255)                                                                                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                >  2:239976/1‑300 (MQ=18)                                                                                  tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                >  2:547361/1‑300 (MQ=18)                                                                                            aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACCCGTCTTATCAGGCCTACGCCCGACCGCGCAATAGCCTGATTTCGCGTGATTTTGTATGTGGGATAAGGCGTTTACGCATCATcctccat                                                                                                                                                                                                      >  1:85214/1‑295 (MQ=11)                                                                                            aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCCGATGCCCCGCTTAACCGCCCTTATCAGGCCTCCTCCACCCACCGCCAGATCATTATTTAGCCTAATCTTTAACGTCGGATCAGGCGtttccgctgcagcagccat                                                                                                                                                                                                      >  1:562008/1‑284 (MQ=11)                                                                                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGTCTCTTAACGCGTCTTATCAGGCCTCCCCCAGCCAGAGCAACACCATGCTTTACCTGGACTGTCAAGGTGGGATCAGTCGtttcgcccggctacaccccccatg                                                                                                                                                                                                 >  1:891387/1‑279 (MQ=11)                                                                                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGCCAgcgc                                                                                                                                                                                                                                                                  >  2:334279/1‑235 (MQ=255)                                                                                                      gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGCTGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                      >  1:622789/1‑210 (MQ=21)                                                                                                      gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCCATAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATACACGCCTggt                                                                                                                                                                                           >  2:640640/1‑299 (MQ=255)                                                                                                      gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                      <  2:622802/210‑1 (MQ=34)                                                                                                      gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCCACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                                                                                                                                     <  2:709384/243‑1 (MQ=18)                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                                                             >  2:508491/1‑298 (MQ=21)                                                                                                       cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCAGACATCAATGCCTGATg                                                                                                                                                                                          >  2:100024/1‑301 (MQ=21)                                                                                                        gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTTATTTTGTAGGTCGGATAAGGCGTTTCCTCCGCATCCGACCTCAATTCCTGATg                                                                                                                                                                                          >  1:113151/1‑300 (MQ=255)                                                                                                         cAATAGCCTGATGTAGCGTGATTTTGTGGGTCGGATAAGGCGTTGACGCCGGATCCGAAATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                                                                         <  2:731334/300‑1 (MQ=255)                                                                                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                                                                     >  2:455059/1‑299 (MQ=21)                                                                                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGTCGTTTATGCCGCATACGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTCCGCCAGACAGCGAACTAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATACGGCATTTCCGCCTCATCCGCCATCaattccgtatgcgccacttgc                                                                                                                                                                                >  2:971836/1‑281 (MQ=11)                                                                                                                  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCACCGCTTAACGCGTCTTATCAGGCCTTCGCCAGACAGCGCACTCGCTTGATTTACCGGGATTTTTTACGTCGGATAAGGCTTTTACCCCGCATCCTACCTCAATGCCTGATGCGACGTttgc                                                                                                                                                                                >  2:643281/1‑298 (MQ=11)                                                                                                                           tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACCTCAATGCCTGATGCGACCCTTGCCCCGTCTTAt                                                                                                                                                                      >  2:386006/1‑301 (MQ=12)                                                                                                                            gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                      >  2:1207574/1‑300 (MQ=17)                                                                                                                             aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                     >  2:335900/1‑300 (MQ=17)                                                                                                                               tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAgcc                                                                                                                                                                                                                                                                       >  2:442440/1‑198 (MQ=21)                                                                                                                               tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATATGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                                                                  >  2:704551/1‑301 (MQ=12)                                                                                                                                ttgggtgcgggtgaaGGGGTTGCCGCCGCATCAGCCGTGAACGCCTGATGAGAGGCTGGCGGGGTCTGATCAGGCTTAAGCGATACAGCGCACTAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTGAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                 <  1:110433/288‑1 (MQ=11)                                                                                                                                ttGTAGGTCGGATAAGGCGGTTACGCCGCAGCCGACATCAACGCATGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                 <  2:795059/301‑1 (MQ=12)                                                                                                                                 tggaGGTCGGTGAGGGCGTTGACGCGGCAGCGGACACAAACGCCTGATGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCGGAAAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCAACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                <  1:335892/298‑1 (MQ=11)                                                                                                                                 tGTAGGTCGGATAAGGCGTTTCGCCCGCATCCGCAATCAACGCCTGATTGGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                <  2:1230632/301‑1 (MQ=11)                                                                                                                                        cGGATCAGGCGTTTACGCCACATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATACCCTGCTTTAGCGTGATTTTGTAGGTCGTCTCACGCGTTTATGCCGCATACGACATCCTCGCCCGATGCGACGCTTAACGCGTCTTCCCCGGCATACGCCTGCCAGCGCACACCCCTGCTTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTCTGCCGCCTCCGCCCGCAATCCCTGCTGCGCCGCATGCCGCCTCTTCTCTGGcttctctta                                                                                                                                                         >  1:915675/1‑293 (MQ=2)                                                                                                                                         ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAAGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtctcttaa                                                                                                                                                        >  2:83036/1‑294 (MQ=255)                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAACCGACATCAATGCCTGGTGCGACGCtttcggctgcttatgagtgctttattaa                                                                                                                                                        >  2:667486/1‑274 (MQ=11)                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCGTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCAAGACAGCGCAATACCCTGATGTAGCGTGATTTTGTAGCGCGGATAAGGCTTTTACGCCGAATGCGAAATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGTCCTatattaa                                                                                                                                                        >  2:670425/1‑294 (MQ=11)                                                                                                                                                                                                                                                             ctaggagaggtgtataagagacataTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                     <  2:32990/276‑1 (MQ=11)                                                                                                                                                                                                                                                             aaaggtttttACGCCGCTACCAACATGAACGTCGTATGCTGCGCTTGCCGCGTTTGATCAGGCCTACGCCAGACAGCGCAAAAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                    <  2:882015/294‑1 (MQ=1)                                                                                                                                                                                                                                                                        cgCCGCATCCGACATCAACGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                         <  2:85216/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                     aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGCTTTTGTGGGTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCTTCTCTTCAGGCCTACCTTAAccctgtatcttatcaacatctgccgaa            >  1:32989/1‑276 (MQ=11)                                                                                                                                                                                                                                                                                     aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc             >  2:275975/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                     aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATAACTTCCTGGTGCGACGCTTGCCGCGTCTTATGCGGACTATCTTCACCGTTGGTTAATTTTTCAACATCTTca            >  2:441966/1‑301 (MQ=11)                                                                                                                                                                                                                                                                                         acgcggtatggtaCGCTTTCCGCGTCTTATAAGGCCTACGGCAGACAGCGCAATAGCCTGATTTATCGTGATGTTGTATGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  2:771912/289‑1 (MQ=14)                                                                                                                                                                                                                                                                                         aCGTCTGATGCGACGTTTGCCGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGTATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  2:259549/301‑1 (MQ=14)                                                                                                                                                                                                                                                                                         aCGCCGGATGCGACGATTGCCGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  2:975075/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                          agcttgagggtaCGCTTGCCGCGGCTTATAAGGCCTACGCCAGAGAGCGCAATAGCATGAGTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTATGCCGGATCTGACATGAACGCCTGATGCGAAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc        <  1:522929/289‑1 (MQ=14)                                                                                                                                                                                                                                                                                                aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTCATTTTTCAAACTCTTCACACCCGTTAt  >  2:458961/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                 GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAT  >  NC_000913/4295892‑4296476

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