Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,304	A→C	14.4%	intergenic (+510/+132)	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,304	0	A→C	14.4%	28.1 ‑0.6 ‑1.6	81	intergenic (+510/+132)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (4/8):  ref base (30/37):  total (35/46)
Fisher's exact test strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.41e-02

GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGG  >  NC_000913/4296011‑4296593                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         gCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCAGCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                                                                                    <  2:550245/300‑1 (MQ=12)     taTTTTGTTGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCTGGATGCGACGCTTGCCGAGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATc                                                                                                                                                                                                                                                                                              <  2:806165/300‑1 (MQ=255)      aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATc                                                                                                                                                                                                                                                                                              >  2:801289/1‑300 (MQ=12)      aTTTTGTAGGTAGGAAAAGGCTTTTACGCCGCATCCGTCATCAACGCCTGAGGCGACTGTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATc                                                                                                                                                                                                                                                                                              <  2:326542/300‑1 (MQ=11)          tggaaaaagggacaGGCGGTTCCCCCGCATCCGACATCAACGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATCAGGcc                                                                                                                                                                                                                                                                                         <  1:668888/288‑1 (MQ=11)          tGTAGGGCGGATAAGGCGTTTATGGCGCATCCGACATACACGCCGGATGCGACGCTTGCCCGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATCAGGcc                                                                                                                                                                                                                                                                                         <  2:345578/301‑1 (MQ=11)          aagagcagtgATAATGGGTTTACGCCGCATCCGACATTAACGGTGGATGCGCGGCTGGCCGGGCTTTAGCAGGCCAACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGTCGCTTGCCGCGTCTT‑ATCAGGcc                                                                                                                                                                                                                                                                                         <  1:32214/292‑1 (MQ=11)                  ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTTCCACATCAATGCCTGCGGCGACCCTCTCCGCTTCTT‑ATCAGCCCTAcattaa                                                                                                                                                                                                                                                                                 >  1:604458/1‑296 (MQ=11)                   gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGTCGCTTGCCGCGTCTT‑ATCAGGcc                                                                                                                                                                                                                                                                                         >  2:32214/1‑292 (MQ=17)                         gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTT‑ATCAGGcc                                                                                                                                                                                                                                                                                         >  2:668916/1‑286 (MQ=12)                                       tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTGTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTT‑ATc                                                                                                                                                                                                                                                                                              <  2:709863/267‑1 (MQ=9)                                                       gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTACGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCGGATTTACCGTCATTTTTTAGCTCGGATAAGGCGGACATGCCCCATCCAGCATCAACCCCTGCTTCGCACCGTCACGCCGATT‑ATCCTGCATACGCTCGACCGCGCAATCGccacctttcaccttatttactagat                                                                                                                                                                                                                                            >  1:690027/1‑278 (MQ=2)                                                                                                                                    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTAGGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                                                                                                                                             >  2:712284/1‑300 (MQ=255)                                                                                                                                      aaaGGAGTTTACGCCGCATCCGAAATCAACGCCGGAGGCGACGCTGGCAGCGTCTGATCAGGCCTACGCCAGACATCGCAATAGCCTGAGTGAGCGTGATTTTGTAGGTCGGATAAGTCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                                                                                          <  2:279309/300‑1 (MQ=11)                                                                                                                                       gcggcaggtgagagggttgtaaagggaaagCCTGATGCGACGCTTGCGGCGTCTTAGCAGGCGTACGCCAGACATCGAAATAGCCTGATGTAGCGTGATTTGGTAGGTCGGATAAGGCGGTGATGCCGCGTCCGACATAAACGCCTGATGCGACGCTTGACGCGTCTG‑AGCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                                                                                          <  1:691143/271‑1 (MQ=11)                                                                                                                                          gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑CTCCGACCTCCATTCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTACCCGTTgg                                                                                                                                                      >  2:686699/1‑301 (MQ=11)                                                                                                                                             tttACGCCGCATCCGCCAACGTCGCCGGATGCAACGCTTACCGCGGCTTAGCGGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                                                                                                    <  2:761358/300‑1 (MQ=11)                                                                                                                                                                        gATGGGGCGCTTGCCGCGTCTTTGCAGGCCTACGCCAGACAGCGCAATGGCCTGATTGAGCGGGATGTGTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                                         <  2:235238/300‑1 (MQ=14)                                                                                                                                                                         aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTTCAAACTCTTCACACCCGGTAt                                                                                                                        >  2:671934/1‑300 (MQ=255)                                                                                                                                                                              acgcttgcagggTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                   <  2:193060/289‑1 (MQ=255)                                                                                                                                                                              acgcatgctgcgTCTTATCAGGCATACTCCAGACATCGCAATAGCCTGATTGAGCGGGATTTTGTAGGTCGGAGAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                   <  2:597872/291‑1 (MQ=14)                                                                                                                                                                              aCGCTTGCCCCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                  <  2:649276/301‑1 (MQ=255)                                                                                                                                                                               cgctggccgcgTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                                                                                                 <  2:350862/294‑1 (MQ=255)                                                                                                                                                                                       cgttCTTATAAGGCTTACGCAAGACAGCGCAATAGCTTGAGTTAGGGTGATTTTGTAGGTCGGTTAAGGCGTTTATGCAGCATCGGACATCAACCCCTGATGCGACGGTGAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATGTTGTATGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  2:626308/298‑1 (MQ=17)                                                                                                                                                                                       cctgctcatcatccgtccgtcaacaaGCGCAATAGCCTGATTTATCGTGAGTGTGTAGGTCGGACAATGCACTTATGCCGCACCCGCCATCAACGCCTGATGCGACGCGTAACGCGTCGT‑ATCAGGCCGACGCCAGACAGCGCAATATTCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  2:656844/276‑1 (MQ=17)                                                                                                                                                                                        cttCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  2:264667/298‑1 (MQ=255)                                                                                                                                                                                        cggtggagagggCCAACGCCAGAAAGGGCAAAAGCGTGATTTAGGGGGATTTGGGAGGTCGGATAAGGCGTTTAGGCAGCAGCCGACATCAACGCCTGATGGGACGCTTAACGCGTCGT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  1:418756/290‑1 (MQ=17)                                                                                                                                                                                        cgTCTTATCAGGCCTACGTCAGAAATCGCAATAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGGGTTGAAGCCGCATCTGAAATAAACGCCTGATGCGACGGTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  2:943772/300‑1 (MQ=17)                                                                                                                                                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGGAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACCTCAATGCCTGCTGCGACGCTTGCCGCTTCTTATCAACCCTATCTTAACTGTTGGTTAATTTTTCACACTCTGCACACCCGGTATCAAACCCTTCCATAccga                                                                                                      >  1:137352/1‑298 (MQ=14)                                                                                                                                                                                           cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCACACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                                                      >  2:829500/1‑300 (MQ=255)                                                                                                                                                                                            ttAACAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGGAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                                                    <  2:40544/301‑1 (MQ=255)                                                                                                                                                                                               tCAGGCCTACGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                                                                                                 <  2:174424/301‑1 (MQ=255)                                                                                                                                                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACCGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAGCGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTAGATTTTGTGGGTCGGACAAGGCGTTTA‑‑CGCACC‑ATCCAGCATCACTGCCGTATGCGACCATTGCCGAGGCTTATCGGTCCTATCTGAACCGTTGGTTAATCTTTCAATCTGTTCCCTCACCGTATAAACCCCTACCATAAAGCTTCGGttt                                                                                              >  1:233267/1‑300 (MQ=11)                                                                                                                                                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTCTGTCGGTCGGATAAGGCGTTTACGC‑C‑GC‑CTCCGAACTCAAGGCCTGATGCCACGCTGCCCGCGTCTTATCAGGCCTATCTTAACCCTTGGTTACTCATTCAAAATCCTCACCCCTCGTTCCAATCCCTTCCATCCACCTGAGTtt                                                                                               >  1:242215/1‑300 (MQ=11)                                                                                                                                                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGTCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTcctgtt                                                                                               >  2:434631/1‑295 (MQ=255)                                                                                                                                                                                                  ggCCTACGCCAGACAGCGCAATAGCCCAATTTAGCGTGATTTTGTAGGTCGGAGAAGGCTGTGATGGCGCATCCGACATCAACGCCTGATGCGAAGCTTACCGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTGTGTAGGTCGGAGAAGGCGTTTACGC‑C‑GC‑AGCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGGCCATACAGCTCAGGtt                                                                                               <  2:233274/300‑1 (MQ=17)                                                                                                                                                                                                   gcttACGCCAGAATGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATCCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                                                                                               <  2:946427/296‑1 (MQ=255)                                                                                                                                                                                                      tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑Acc                                                                                         <  2:831141/301‑1 (MQ=255)                                                                                                                                                                                                      tACGCCAGAAAGCGCAATAGATTGAGTTGGGGGTATTGTGGAGGTTGGATAAGTCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GT‑AGCCAACATCAATGCCTGATGCGACGCTTGCCGAGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑Acc                                                                                         <  2:680873/301‑1 (MQ=17)                                                                                                                                                                                                       aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGATTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCCACCCGTTATCCACCCATTCCCTACCGCTCAGGTTCC‑ACCa                                                                                        >  2:826493/1‑301 (MQ=14)                                                                                                                                                                                                       aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGC‑C‑GC‑ATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATACAACCCTCTCATCCACCTCAGATTCA‑Acc                                                                                         >  1:532289/1‑300 (MQ=14)                                                                                                                                                                                                       aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑CTCCGACATAAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTACTTTTTCAAACCTTTCACCACCTGTATCCAACCCTTCCCTAACGCTCACGtctaccca                                                                                         >  1:587586/1‑294 (MQ=14)                                                                                                                                                                                                       aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAACCTCTTCACACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCA‑Acc                                                                                         >  2:4589/1‑300 (MQ=255)                                                                                                                                                                                                            agcaggcgaaagggctgtaTTTAGCGGGAGTTGGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCGT‑ATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGc                                                                                    <  1:400027/282‑1 (MQ=18)                                                                                                                                                                                                               cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCACACTCTTCACACCAGGTAACAGACCCTTACATACAGATCAAGtcacacacgttaccc                                                                                  >  1:576302/1‑286 (MQ=14)                                                                                                                                                                                                                 gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCACCCCGGTATCAAACCCTTCCCTACAGCCCAGGTTCA‑ACCAGTGCAgcgc                                                                               >  2:721286/1‑300 (MQ=255)                                                                                                                                                                                                                 gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGC‑CCAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACACTTCACTCCCGGTATCACACCCTTCCATTCACCTCACGTTCA‑Acgcgttcccccc                                                                               >  1:791879/1‑289 (MQ=14)                                                                                                                                                                                                                         gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGCCTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCAAGCTAAGGTTCA‑ACCAGTGCAACGCCTTTTGCtt                                                                      >  2:715589/1‑301 (MQ=255)                                                                                                                                                                                                                          gcagatgtatagttaTTTTGGAGGGCGGATAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATTCGACGCTTAACGCATCTT‑AGCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCt                                                                       <  2:372606/285‑1 (MQ=18)                                                                                                                                                                                                                          ccggATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTCCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCt                                                                       <  2:820115/296‑1 (MQ=255)                                                                                                                                                                                                                             gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGACGCCCCTTTTGCTTTttt                                                                  >  2:608732/1‑300 (MQ=255)                                                                                                                                                                                                                             gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                  >  2:677392/1‑254 (MQ=255)                                                                                                                                                                                                                                          ttttgaGGTCGGATAATGCTTTTATGCCGAATCCGACATCAACGCCTGATGCGACGCTTAACGCGTATT‑AGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                                                     <  1:721257/296‑1 (MQ=255)                                                                                                                                                                                                                                          tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                                                     <  2:604477/301‑1 (MQ=255)                                                                                                                                                                                                                                          ttggTAGGTGGGATAAGGCGTTTAGGCCGCACCCGACATCAACGCCTGATGCGACGCTTAACGCGTCGT‑ATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                                                     <  1:801260/298‑1 (MQ=255)                                                                                                                                                                                                                                                  tgggATAAGGCTTTTAGGCCGACACCGACACCGAAGCCGGATGCGACGCTTAACGGGGCTT‑ATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGGTAAGGCGTGTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTGCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                                              <  1:686671/298‑1 (MQ=21)                                                                                                                                                                                                                                                   cGGATAAGGCGTTTATGCCGCAGCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                            <  2:587603/301‑1 (MQ=255)                                                                                                                                                                                                                                                    ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                         <  2:22573/240‑1 (MQ=255)                                                                                                                                                                                                                                                    ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTATCTTAAAAGTTGGTTCATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACACCTCGTGTTCA‑CCCACTGAAGCTCCTTTTGATTGTTCTGCTAGCTATAACCTTCTTCAc                                            >  2:363366/1‑300 (MQ=17)                                                                                                                                                                                                                                                     gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGATTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGCAGC‑ATCCTGAATCGATGCCTGATGTGACGCTTGCCTAGACTTATCAGGCTTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACCCCAGGTATCACACCCTTCCATACCGCTCAGGTCCACACCA‑TCCAGCACCTTTTGCTCATGTTTCTCGCTACACCCTTCATCCcc                                           >  1:440916/1‑300 (MQ=17)                                                                                                                                                                                                                                                     gATAAGGCG‑TTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTCACCGTTGGTTCATTTTTCAAACTCTTAACACCCGGTATCAAACCCTTCCATAAAGCTCCGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTAAAACCTTCTTCAcctt                                         >  2:21114/1‑299 (MQ=18)                                                                                                                                                                                                                                                      gtaagcgggttAGGCCGCATCCGAAAGCAACGCCGGATGCGACGCTTAACGCGTCTT‑ATGAGGCCTCAGCCAGACAGCGCAATGGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCg                                          <  1:378142/291‑1 (MQ=21)                                                                                                                                                                                                                                                           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTCCAGCCCCTTTTGCTTGTTCTTCTCGATCAAACCTTCTTCCCCGTTTaca                                    >  1:826781/1‑299 (MQ=255)                                                                                                                                                                                                                                                            cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCGGACAGCGCACTAGCCTGATTTAGCGTGATGCTGTAGGTCGGATAAGTCGTTTACGC‑C‑GC‑ATCCGATATCAATGCCTGGTGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTCATTTTTCAAACTCTTGACCCCTGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCCTGATAAAACCTTCTTCACCGTTTCAGa                                   >  2:774174/1‑301 (MQ=21)                                                                                                                                                                                                                                                            cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCCATAGACTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACAGCACAGTTTCA‑ACCAGTGCAGCGCCTTTTTCTTTTTTTTCCTGATCACACCTTTTTCCCCTTTTacca                                   >  1:362888/1‑298 (MQ=18)                                                                                                                                                                                                                                                              tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAc                                 <  2:826809/301‑1 (MQ=255)                                                                                                                                                                                                                                                              ttgagcccGCACCCCAAATCACGGCCTGAGGCGACGCTTAACGCGGCTT‑ATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATAAGGCGGTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                  <  1:412349/294‑1 (MQ=25)                                                                                                                                                                                                                                                                      ccaccccaaatcacagCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCCACGCCAGACAGCGCAATAGCCTGTTTGAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTcc                         <  2:362899/286‑1 (MQ=25)                                                                                                                                                                                                                                                                      ccaaccGACCTTAATGCCTCATGCGACGGTTAGCG‑GTGTTGAACAGAGTTACGCGAGAAAGCGCAAGAGCCTGATTTTGCGGGATTTGGTCGGTCGGATAAGGCGTTTACGC‑C‑GA‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCATAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAGACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATCcc                         <  1:130879/297‑1 (MQ=25)                                                                                                                                                                                                                                                                       cATCCGACATCAACGCCTGATGCGACGCTAAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTcc                         >  2:241302/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCCTCACCGTTTAAGAACATCATTCCCg                       >  2:875843/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGACAGACAGCGCAATAGCCTGAATTAGCGTGATTTTGTAGATCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACTGCAAGGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCGATCTTAACCGTTGTTTAATTTTTCAAACTCTTCCCACCCAGTATCAAACCCTTCCTTACAGCTCAGCTCCA‑CCCCGTGCACCGCTTTTTACGTTTTTTTCTCGATAAACCCTTCTTCACCGTGTAAAAAAATACTTccac                       >  1:561791/1‑298 (MQ=17)                                                                                                                                                                                                                                                                             aCATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGTTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAACCGGTATCAAACCCTTCCATACAGCTAAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAACCCTTCTTCCCAGTTTAAGAACATCATTCCCGCCCat                  >  1:661527/1‑300 (MQ=255)                                                                                                                                                                                                                                                                             aCATAAACGCCGGCGGCCACGCTTAAGTCGTCTT‑AAGAGGCCTACGCCAGACAGCGCAATAGGCGGATTTAACGGTATTTGGGAGGGCGGAGAAGGCGTTTAAGC‑C‑GC‑AGCCGACATCAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAg                  <  1:826465/301‑1 (MQ=255)                                                                                                                                                                                                                                                                              cATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCC‑ACCAGTGAAGCGCCTTTTGCTTGTTCTTCTCGTTacaaa                                                     >  2:471559/1‑264 (MQ=255)                                                                                                                                                                                                                                                                                       cgggtgcgcaggttAACGCGTATT‑ATCAGGGCGAAGCCAGACAGCGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATACGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt        <  1:773403/289‑1 (MQ=37)                                                                                                                                                                                                                                                                                          aTGCGACGCTTAACGCGTCTT‑ATCAGGCATACGACAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGCAGC‑ATCCGACATCAATGCCTGAAGCGACGCTTGCCGCGTCTTATCAGGCGTATCTTAACCGGTGGTTAATTTTTCAAACTATTCACAGCCGGTATGAAAGCCTTCGCTACAGCTCAGGTTCA‑ACGAGTGCAGGTCCTTTTGCTTGGTCTTCTAGATAAAACCTGTTTCACCGTTTCAGAACATCATTTCAGACCAGCACTCGGAATaag     >  2:162029/1‑300 (MQ=18)                                                                                                                                                                                                                                                                                           tggcaCGCTTAACGCGTCTT‑ATCAGGCCTACGGCAGAAAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATAAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc    <  1:715560/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                           tGCGATGATTAACGCGTGTG‑ATCAGGCCTACGCCAGACAGAGCAATAGCCGGATTTAGCGTGATTGTGTGGGTCGGATAAGGCGTTTACGC‑C‑GC‑TACCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCGTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc    <  1:829472/301‑1 (MQ=37)                                                                                                                                                                                                                                                                                            gCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc    <  2:532305/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                              taCGCTTAATGATGCTT‑ATCAGGCCTACGCCAGACAGCGCAGTAGCCGGATTGAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTGCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCgg  <  1:608713/299‑1 (MQ=37)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         GCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTT‑ATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGG  >  NC_000913/4296011‑4296593

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