BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	96.5%	344.0 / 4.5	115	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (49/61); total (50/65)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	4,296,380	2	.	G	96.5%	336.8 / 4.2	115	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand): ref base . (1/3); new base G (48/61); total (50/65)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAA  >  NC_000913/4296169‑4296621

||

tgaagtcgtacggcggcgGTTAGCGTGTCTTGTCAGGCCAACGCCACACAGCGCACTAGCATGATCT‑GGCGTGGTTTAGTCGGTGGGGTAAGGCGTTTATGCCGAATCCGCCGTCTACGCCTGCTGCGACCATTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAACCTGACGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                                                                                              <  2:605956/239‑1 (MQ=11)

                             ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                                 <  1:967481/251‑1 (MQ=17)

                                  aGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                            >  1:831678/1‑251 (MQ=21)

                                       tACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGATCATTTTTc                                                                                                                                                                       >  1:401824/1‑251 (MQ=11)

                                        gccccAGACTGCGCAATAGCCTGGTTT‑GGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                      <  2:919041/248‑1 (MQ=255)

                                        aCGCCCGACAGTGCAATAGCCTGAGTT‑AGCGTCATTTTGTAGGTCGAATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                      <  1:896284/251‑1 (MQ=11)

                                        aCGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                      <  1:219067/251‑1 (MQ=255)

                                        aCGCCAGACAGCGCAAAAGCCTGTTTT‑CGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATTTTAACCGTTGGTTAATTTTTCa                                                                                                                                                                      <  2:347102/251‑1 (MQ=255)

                                         gaccAGACAGCGAAGTGGCCTGGTTG‑AGCGTGATTTTCGAGGTAGGATAAGGCTTTTCTGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                     <  2:357666/249‑1 (MQ=11)

                                         cGCCAGTAAGCGCATGAGGCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATCCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                     <  2:753520/251‑1 (MQ=11)

                                             gtaCAGCCCAATGGTCTGCTTT‑GGCGAGATTTTCTAGGACGGACAAGGCCTTACGGGCGCATCCTAAATCAAACGCTGGTGCGACGCTTACCTCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCAGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                                                                                                                                 <  2:491186/249‑1 (MQ=14)

                                                  gcgcAAGTGCCAGATTT‑AGGGTGGTGTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAGCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTGcac                                                                                                                                                            <  2:1057521/251‑1 (MQ=14)

                                                   ccgccggagcctgctttagatgtgttttttAGCTAGGTTAAAGTGTTTATACCTCTCCCAACATCAACGCCTGGGGCGACGTTTTACGCGTCTTATCCCGCCTGCGCCATGCAGCGCAATTGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCTTCCGACATCAATGCCTGATCCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTcac                                                                                                                                                            <  2:302910/227‑1 (MQ=9)

                                                      aaTAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacactc                                                                                                                                                        >  1:94998/1‑249 (MQ=255)

                                                      aaTAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCCACCCAGCGCTATAGCCTGCTTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTT‑‑ACGC‑C‑GCGTCCAACATCACTGCCAGAGTCATGGTCTGACGCGTCTTAGCAGGCCTCTCTTAAACGTTGGTTAATATTTCCAGCTCTTCACAccc                                                                                                                                                        >  2:253848/1‑251 (MQ=11)

                                                       aTAGCCTGATTT‑GGCGTGATTTTGTAGATCGGATAAGGCGTTTATGCCGCATCCGACATCGACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                                       <  1:193371/251‑1 (MQ=255)

                                                       aTAGCC‑GATTTGAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACAGCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                                       <  1:746596/251‑1 (MQ=14)

                                                           ccTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGACTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCTAGACAGCGGAATAGCCTGCATTAGCGTGATTTTGTAGGGCGGAGAAGGCGTTT‑‑ACGC‑C‑GCGTAAGAAATGAACGAGTGATACGACGCTTTCCGCGTCTTATTAGGACTATCCTAACCGTTGGATAATATTTCATAATCCTCacactaggcca                                                                                                                                                   >  2:42617/1‑244 (MQ=11)

                                                            cTGATTT‑TGCGTGATTTGTTGTGTCGCATAAGTCGTTTATGCCGCGTCCGAGCTCAACGCCTGTTGCGGCGCTTAGCTCGTGTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                                  <  2:95044/251‑1 (MQ=14)

                                                            cTGATTT‑AGCGTAATTTTGTATGTCGGATAAGGCGTTTATGTCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                                  <  1:141318/251‑1 (MQ=255)

                                                                ttt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                            >  1:694200/1‑237 (MQ=255)

                                                                ttt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                              >  1:787492/1‑251 (MQ=255)

                                                                ttt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                              >  1:728520/1‑251 (MQ=255)

                                                                ttt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGACGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAAc                                                                                                                                              >  1:268554/1‑251 (MQ=255)

                                                                ttt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGAAATCAATGCCTGATGTGACGCTTTCCGCGTCTTATCAGGCCTATTTTAAGCGTTGGTTAATTTTTCAAACTCTTCACACCAGGTATCAAAc                                                                                                                                              >  2:70216/1‑251 (MQ=14)

                                                                 tt‑AGCGTGTTTTTGTAGGTCGGATAAGGCGTTTGGACCGCATCCGCCATCACCGCCTGAGTCGACTCTTAACGCGGTTTATCAGGCCGAGGCCAGACAGCGGAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCGATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  2:1028810/251‑1 (MQ=17)

                                                                 tt‑AGCGTGATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATACGACATCAACCCCTGATGCGACGCGTAACACGTCTTATCAGGCCTACGGCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  1:560473/251‑1 (MQ=17)

                                                                 tt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  1:147161/251‑1 (MQ=255)

                                                                 tt‑AGCGTGATTTTGTAAGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  1:316653/251‑1 (MQ=255)

                                                                 tt‑AGCGTGATTTTGGAGGTCGGGTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  2:86898/251‑1 (MQ=255)

                                                                 tt‑AGAGTGATTTTGTAGGTCGGTTAAGGCGTTTATTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  2:185688/251‑1 (MQ=255)

                                                                 tt‑ACCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  1:587021/251‑1 (MQ=255)

                                                                  ttcgCGTGATTTTTGAGGTGGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             <  2:340513/248‑1 (MQ=255)

                                                                           ttttGTAGGTCGGATAAGGCGTCTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAAGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATCCAGCGTGACTTTGTAGGTCGGATAAGGCGTTT‑‑ACTC‑C‑GAATCCGACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTATCCGTTGGTTAATCTTTCAAAATCTTCACACCCGGTATCCAACCCTTACATAc                                                                                                                                    >  2:929817/1‑251 (MQ=14)

                                                                            ttttaaGGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑C‑C‑GCATCCAACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                                                                                                   <  1:122406/246‑1 (MQ=17)

                                                                            tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑A‑‑CGCTGCATTCGACCTCAATGACTGATGCGACGCTTGCCGTCTCTTATCAGGCCTATCTTAACCTTTTGTTCATTTTTCCAACTCTTTACACACTGTAACATATCCTTCCATaaa                                                                                                                                   >  2:820380/1‑249 (MQ=14)

                                                                            tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTGGCCGCGTCGTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCGGGTATCAAACTCTTCCATACa                                                                                                                                   >  2:254842/1‑251 (MQ=17)

                                                                            tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                                                                                                   <  1:146154/251‑1 (MQ=255)

                                                                                  ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAg                                                                                                                             >  1:185632/1‑251 (MQ=255)

                                                                                   gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GAATCCGAAATTAATGCCTGATGGGTAGCTTTACGCGTCTTATCAGGCCTAGCTTAACCGGTGGTTAATTTTTCAAACCCTTCACACCCGGTAAACAACCCTTCAATACAGCTCAgg                                                                                                                            >  2:358177/1‑251 (MQ=14)

                                                                                       gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                                                                                                                        >  1:903900/1‑251 (MQ=255)

                                                                                       gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACAGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCATTCCATACAGCTCAGGTTCa                                                                                                                        >  2:106716/1‑251 (MQ=255)

                                                                                       gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCCTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                                                                                                                        >  2:319597/1‑251 (MQ=255)

                                                                                             gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt                                                                                                                  >  1:1017920/1‑251 (MQ=255)

                                                                                              cGTTTATGCCGCATCCGACATCAACGCCTGATGCGAAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg                                                                                                                 <  1:482568/251‑1 (MQ=255)

                                                                                                    tGCCGCATCCGACATCAACGCCTGCTGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                                      >  1:643541/1‑176 (MQ=11)

                                                                                                      ccgacaccGACATCACCGCGCGATGCGACGCTTAACGCGTCTTATCAGGCCTACGAGAGACGACGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTT‑‑ACTC‑C‑TCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                                         <  2:904066/245‑1 (MQ=21)

                                                                                                       cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                                                                                        <  1:305038/251‑1 (MQ=255)

                                                                                                               aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                                >  1:895194/1‑251 (MQ=255)

                                                                                                               aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                                <  1:106670/251‑1 (MQ=255)

                                                                                                                 tataagagacagggCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                                                                              <  1:560715/238‑1 (MQ=255)

                                                                                                                 aTCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                     <  2:150429/180‑1 (MQ=255)

                                                                                                                      cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCC‑ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                        <  1:645291/251‑1 (MQ=255)

                                                                                                                         tgtaTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGAt                                                                                      <  1:880811/248‑1 (MQ=255)

                                                                                                                          ttaTGCGCCGCTTGACGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATCAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa                                                                                     <  2:498304/249‑1 (MQ=25)

                                                                                                                                 acgcgtatgggttctgagcAGGCCTCAGCCAGATAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTCCG‑‑C‑C‑GCCTCCGACATCAATGCTTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACctt                                                                              <  2:848764/233‑1 (MQ=25)

                                                                                                                                  cggttAACACTTGTTGTCAGGCCTACGACAGACAGCGCAAGAGCCTGTGTTAGCGTGATTTGGTAGGTCGGCTAAGGCGTTT‑‑ACGC‑C‑GGATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGGCTATCTTGACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCATGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                             <  2:502634/248‑1 (MQ=25)

                                                                                                                                  cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                             <  1:669150/251‑1 (MQ=255)

                                                                                                                                  cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                             >  1:696724/1‑251 (MQ=255)

                                                                                                                                    cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttctt                                                                           <  1:70172/251‑1 (MQ=255)

                                                                                                                                     ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                                                                          <  1:171890/251‑1 (MQ=255)

                                                                                                                                       aaCTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAGCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                        <  1:208150/251‑1 (MQ=255)

                                                                                                                                       aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                        <  2:895360/251‑1 (MQ=255)

                                                                                                                                       aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTACTGTTCAACCAGTGACGCGCCTTTTGCTTGTTCTTCTCGCTAAAACTTTCTTCAc                                                                        >  2:124899/1‑251 (MQ=21)

                                                                                                                                        aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAcc                                                                       >  2:646256/1‑251 (MQ=255)

                                                                                                                                        aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCAAACTCTTGCCACCCGGTATCAAACCCTTCCATCCAGCCCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAAACCTTCTTCAcc                                                                       >  2:378930/1‑251 (MQ=21)

                                                                                                                                         cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCATCGCAGCGCCTTTTGCTTGTTCGTGTCGATAAAACCTTCATCACCg                                                                      >  2:171943/1‑251 (MQ=21)

                                                                                                                                         cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCg                                                                      >  1:366122/1‑251 (MQ=255)

                                                                                                                                         cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                             >  1:150377/1‑180 (MQ=255)

                                                                                                                                              cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                                 >  1:796671/1‑251 (MQ=255)

                                                                                                                                                 accAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                              <  1:42575/249‑1 (MQ=255)

                                                                                                                                                 aaagagacagcGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCTCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                              <  2:503495/241‑1 (MQ=255)

                                                                                                                                                 aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                              <  2:1018107/251‑1 (MQ=255)

                                                                                                                                                     ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTTTTCTCGATAAAAACTTCTTCACCg                                                                      >  2:560826/1‑239 (MQ=255)

                                                                                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATTAGGCCTATTTTAATCGTTGGTTAATTCTTCAAACTCCTCACACCCGGTATCAAACCCGTCCATACAGCCCAGCCTCACCCAGTTCGGTTCCTTTTGCTTGTTTTTCCCGCAAAATCCTCCTTCACCGTTTCAGAACTTCAtttc                                                     >  2:325299/1‑249 (MQ=18)

                                                                                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                                                                                                                                  <  1:789329/174‑1 (MQ=255)

                                                                                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                              <  1:604928/242‑1 (MQ=255)

                                                                                                                                                            gccagccaACGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACCTCCATTCAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCg                                                   <  1:253775/251‑1 (MQ=255)

                                                                                                                                                              cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                                 >  1:1026680/1‑251 (MQ=255)

                                                                                                                                                                      gccaTAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAATGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACATCCTGTATCAAACCCCTCCATACCGCCCAGGTTCAACCACTGCAGCCCCTCTTGCTTGTTCTTCTCGATAAAATCCTCTTCACCGTTTAAGCATATCATTCCCGCCACGTACTc                                         >  2:77465/4‑251 (MQ=21)

                                                                                                                                                                       cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTCAAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCg                                        >  2:757248/1‑251 (MQ=255)

                                                                                                                                                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGGACTCgg                                       >  2:2104/1‑251 (MQ=255)

                                                                                                                                                                        aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCCATAAAACCTTCTCCACCGTTTAATAACATTATTCACGCCCAGTaactgg                                       >  1:498790/1‑246 (MQ=37)

                                                                                                                                                                            gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                        <  1:421610/230‑1 (MQ=255)

                                                                                                                                                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCTTATCTTAATCGTTGGTTAATTTTTCAAATTCTTCACACCCGGTATCAAATCCTTCCATACAGCTCAGGTTCATCCACTGCAGCGTCTTCTGCTTTTTCTTCTCGATAAAACCTTCTTCACCGTTTAATAACATCATTACCGCCCAGTACTCGGAATaa                                  >  1:384424/1‑251 (MQ=37)

                                                                                                                                                                                  tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGACCAGGTTCAACCAGTGCAGCGCATTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTACGAACACCATTCCCGCCCAGTACTCGGAATAAACGGt                             >  1:246445/1‑251 (MQ=37)

                                                                                                                                                                                  tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACCCGGa                                      >  2:605040/1‑242 (MQ=255)

                                                                                                                                                                                  tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAACAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                    >  2:909995/1‑148 (MQ=255)

                                                                                                                                                                                   ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCTCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACCCGGa                                      >  1:503387/1‑241 (MQ=255)

                                                                                                                                                                                     aGCGTGATTTTCTAGTTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCCTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTAGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg                          <  1:205625/251‑1 (MQ=255)

                                                                                                                                                                                      gCGTGAGTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGGCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAACCGTTTGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCgg                         <  2:818457/251‑1 (MQ=255)

                                                                                                                                                                                           aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                                                                <  1:765766/175‑1 (MQ=255)

                                                                                                                                                                                            ttttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAAtt                   >  1:602927/1‑251 (MQ=255)

                                                                                                                                                                                              ttGTAGTTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGc                 <  2:603039/251‑1 (MQ=255)

                                                                                                                                                                                              ttGTAGATCGGATAAGTCTTTT‑‑ACGC‑C‑GCAGCCGCCATCAATGCCTGAGGCGACGCTTGCCCCCTCTTATCAGGCCTATCTTAACCTTTGGTTAATTTTTCAAACTCTTCACTCCCTGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCTGGAAATTGc                 <  2:617006/251‑1 (MQ=255)

                                                                                                                                                                                                   ggTCGGATAAAGCGTTT‑‑ATGC‑C‑CCAGCCGACAGCAACGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGc            <  2:1049756/251‑1 (MQ=255)

                                                                                                                                                                                                    gTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACCCGGAATAAccctctcctttaaacaactctga           >  2:421707/1‑230 (MQ=38)

                                                                                                                                                                                                     tCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTg          <  1:319518/251‑1 (MQ=255)

                                                                                                                                                                                                        gATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGt       <  1:757109/251‑1 (MQ=255)

                                                                                                                                                                                                             ggCGTTT‑‑ACGC‑C‑GCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGaa  <  1:358096/251‑1 (MQ=255)

||

TCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAA  >  NC_000913/4296169‑4296621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`