BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,380:1	+CG	100%	intergenic (+586/+56)	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.3%	300.1 / 0.6	106	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 (51/49); total (54/52)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*	NC_000913	4,296,380	2	.	G	96.3%	307.9 / 0.5	106	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 (53/48); total (54/52)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 106 total aligned reads displayed.

CGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGA  >  NC_000913/4296210‑4296609

||

          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATTAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAAGCCTACGCCAGACGGCGCCTAAGCCTTATTTAGCGTGATTCTGTAGGTCGGGTAAGGTGTTG‑‑ACGC‑‑CGCATCCGACATCAATTCTTGATGCCGCGCTTCCCGCGCTTTTTCAGCCCTATATTACCAGTGGGTTAATTTTTc                                                                                                                                                             >  2:867306/1‑239 (MQ=11)

          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGTGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                             >  1:830347/1‑239 (MQ=255)

          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                             >  1:186330/1‑239 (MQ=255)

          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATCCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                             >  1:785993/1‑239 (MQ=255)

          cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCACATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                             >  1:27065/1‑239 (MQ=255)

             aaTGGCCTGATTTAGCGGGATGTGGTAGGTCGGATAAGGCGTTTTTGCCGCCTCCGACATCAACGCCTGATGCGACCCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                                                                                          <  1:662633/239‑1 (MQ=14)

             aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                                                                                          <  1:1140772/239‑1 (MQ=255)

             aaTAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                                                                                          <  1:1151623/239‑1 (MQ=255)

                    tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGCTTTAGCGTGCTTTTTGTGGTCGACTAAGTCGTTTCCG‑‑C‑‑CCCACCTGACATCAATCCCGGAAGGACCGCTTGCCGCGTCTTACCCGGCCTATCTTAACCCTTGGTTACATTTTCCAAcctttc                                                                                                                                                    >  2:206474/1‑233 (MQ=11)

                    tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                                   >  1:15714/1‑239 (MQ=255)

                    tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                    >  1:213700/1‑238 (MQ=255)

                    tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTAATTCACACCTTCCCCCGTCTTTTCAACCCAATCTTAACCGTTGGTTAATTTTTTAAACCCttaa                                                                                                                                                   >  2:381100/1‑237 (MQ=11)

                    tGATTAAGCGTGTTTTTGTAGAGCTGATAGGGCGTTTATGCCGCATCCGACCTGAACGCCTGATGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                                   <  1:867306/239‑1 (MQ=14)

                      aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGTCCGGCTACGGCGTTT‑‑ACGC‑‑CGCATCCTACATCAATGACCGATCCGACGCGTGCCTCGTCTTATCAACCCTCTCTTAACCCTTTGTTAATTTTTCAAActcccaaaa                                                                                                                                                 >  2:855230/1‑233 (MQ=11)

                        ttAGCGGGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                               <  1:1116432/239‑1 (MQ=255)

                         tAGGGTGGTTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                                                                              <  1:584777/239‑1 (MQ=14)

                         tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                                                                              <  1:341475/239‑1 (MQ=255)

                            cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCTTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGCCTTTTCAGGCCTATCTTAACCGTTGGTTAATTTTTTAAACCCTTCCAACCCgga                                                                                                                                           >  2:732649/1‑238 (MQ=14)

                            cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACTTCAATGCCTGATCCGACGCTTGCCGCGCCTTTTCAGGCCCTTCTTAACCGTTGGTTAATTTTTAAAAACCTTCACACCCGGt                                                                                                                                           >  2:6048/1‑239 (MQ=14)

                            cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                                                                                                                           >  1:1293443/1‑239 (MQ=255)

                                  tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                                     <  1:1087067/239‑1 (MQ=255)

                                  tttGGAGGTTGGAGAAGGCGTTTATTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                                     <  1:801218/239‑1 (MQ=255)

                                          tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGGTGCGACGCTTGCCCCGGCTTACCAGGCCTATCTTAACCTGTGGTTAATTTTTTAAACTCTTCAAACCCGGGACCAAA‑ACCTTCCa                                                                                                                            >  2:573101/1‑239 (MQ=14)

                                            cgATAAGGCGTTTAACGCGCATCCGGCATCAACGGCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATa                                                                                                                          <  1:1291072/238‑1 (MQ=17)

                                                aaGGGGTTTAGTCCGCATCCGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGc                                                                                                                      <  1:423576/239‑1 (MQ=255)

                                                 aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:381100/239‑1 (MQ=255)

                                                 aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:6048/239‑1 (MQ=255)

                                                 aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:302990/239‑1 (MQ=255)

                                                 aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAACAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:605589/239‑1 (MQ=255)

                                                       ttctGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACTGGTTTTATCAGGCCTACGCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGt                                                                                                                <  1:896769/235‑1 (MQ=18)

                                                        tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCACTGCCTGATGCGACGCTTGCCGCGCCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTCAAAACTCTTCACACCCGGAATCCAA‑CCCTTCCATGAACGTCAGGTTc                                                                                                              >  1:225117/1‑239 (MQ=14)

                                                        tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAAACCCGGGATCAAA‑CCCTTCCATACACCCCAGGTTc                                                                                                              >  2:711794/1‑239 (MQ=255)

                                                           gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAAc                                                                                                           <  1:516013/239‑1 (MQ=255)

                                                           gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAAc                                                                                                           <  1:206474/239‑1 (MQ=255)

                                                                 tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAATTTTTCAAACTCTTCACACCCCGTACCAAA‑CCCCTCCATACACCCCCGGTTCAACCaccc                                                                                                      >  2:516013/1‑235 (MQ=17)

                                                                     aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAgt                                                                                                  >  1:271410/1‑237 (MQ=255)

                                                                        tccaCGCCTGATGCGATGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                              <  1:1030514/236‑1 (MQ=255)

                                                                        tCAACGTCTGATGCGATGCTTAACGCGTCTTAGCAGGCCTACGGCAGGCAGCGCAATAGCCTGATTGAGCGTGAGTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAATCTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                              <  1:284294/239‑1 (MQ=18)

                                                                            cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtttt                                                                                           >  1:289607/1‑238 (MQ=255)

                                                                                  aTGCGAGGCTTAACGCGTCTTATCAGGCCTACGCCAGGCAGCGCAATACCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTT‑‑ACGC‑‑CGCATCCGCCATCAATGACTGCTGCGCCGATTACCGCGTCTTATCGGACCTATCTTAACCGTTGCTTCATTTTTCAACCTCTTCACACCCGCCATCAAA‑CCCTTCCATACAGCTCACGTTAACCCATTGCACCGCCTTTTAC‑Ttttt                                                                                   >  2:291106/1‑239 (MQ=17)

                                                                                      gaccgttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGttctt                                                                                <  1:82217/233‑1 (MQ=255)

                                                                                      gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGttctt                                                                                <  1:825055/238‑1 (MQ=255)

                                                                                      caCGCTTTAGGTGTCTTATTAGGCCGACGGGAGGCACCGCAATGGCCTGATTTCGCGTGAGTTTGTGGGTCGGAGAAGGGGTTT‑‑ACGC‑‑CGCAACCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGttctt                                                                                <  2:107693/237‑1 (MQ=21)

                                                                                       ccGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGCCATCAATGCCTGATGCGCCGCTGGCAGCGTCTTATCAGGCCTATCTTCACCGTTGGTTAATTTTTCTAACCCTTCCCACCCGTGCTCCAC‑CCCTTCCATACAGCTCAGGTTCACCCAGTGCATGCCGTTTTGT‑TTGttctttt                                                                              >  2:836007/2‑237 (MQ=17)

                                                                                       aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGttcttc                                                                               <  1:1240982/238‑1 (MQ=255)

                                                                                        cgcgttaagcgTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTtctc                                                                             >  1:475143/9‑239 (MQ=25)

                                                                                              aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATaaa                                                                       <  1:882845/239‑1 (MQ=255)

                                                                                              aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATaaa                                                                       <  1:732649/239‑1 (MQ=255)

                                                                                                  gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACct                                                                   >  1:595442/1‑239 (MQ=255)

                                                                                                   tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATCCAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGCTAAAACct                                                                   >  1:581844/1‑238 (MQ=255)

                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTGTTCACACCCGTGATCAAC‑CCCTTCCATACACCTCAGGTTCAACCAGTTCAGCCCCTTTTTCTTTGTTC‑TCTCGAAAAAACct                                                                   >  2:1286397/1‑237 (MQ=18)

                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:912472/188‑1 (MQ=255)

                                                                                                    cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     >  2:912472/1‑188 (MQ=255)

                                                                                                        ttaGGCCCTCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTATCCGTTGGTTAATTTTTCAATCTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCa                                                             <  1:505803/237‑1 (MQ=25)

                                                                                                        tCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGTATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCttcttcg                                                             <  1:1180425/239‑2 (MQ=255)

                                                                                                         cAGGCCTACGCCAGACAGCGTAATAGCCGGATTTTGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCAc                                                            <  1:455019/239‑1 (MQ=255)

                                                                                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                  >  2:795129/1‑136 (MQ=255)

                                                                                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                  <  1:795129/136‑1 (MQ=255)

                                                                                                            gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGAAGCTTGCCGCGTCTTATCAGTCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATAAAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTTTTCTTGATAAAAATTTCTTAACCGt                                                         >  2:231924/1‑239 (MQ=25)

                                                                                                               tACGCCAGACAGCGCAATAGCCTTATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa                                                      <  1:231924/239‑1 (MQ=255)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCGGCTACGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCTTCAAA‑CCTCTCCATACCTCTCCGGTTCACCCAGTCCAGCGCCTTTTGT‑TGGTTCTTCCTGAAAAAACC‑TCATCACCAGTTAAGa                                                   >  2:110741/1‑239 (MQ=18)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCATTGCAGCCGCTTTTGC‑TTTTTTTTCACCACAAAACCTTCTTCACCGTTTAAg                                                    >  2:470572/1‑239 (MQ=21)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCCTACAGCCCCGGTTCAACCAATGCAGCGCCTTTTTC‑TTGGTCTTCCCGATAAAACCATCTTCACCGTATaaa                                                    >  2:409092/1‑238 (MQ=21)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCt                                                                                                                     <  1:993971/175‑1 (MQ=255)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                                    >  1:264091/1‑239 (MQ=255)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCCCGATAAAACCTTCTTCACCGTTTAAg                                                    >  2:538334/1‑239 (MQ=255)

                                                                                                                 cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACAATTGCACCGCCTTTTTC‑TTTTTCTTTTCGATAAAACCTTTTTCACCGTTAAAg                                                    >  2:796038/1‑239 (MQ=25)

                                                                                                                      gccAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGCTT‑‑ACGC‑‑CGCTTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTATCTTAACCGGTGGTTAATTTTTCAAACTCTCCACACCCGCTATCAAA‑CCCTTCCCTACCGCCCAGGTTCAACCAGTGCAACGCCTTCTGC‑GTTTTCCTCCCAATCAAACCTTCTTCACCGTTCAAGAAaat                                               >  2:1231457/3‑239 (MQ=21)

                                                                                                                      gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                               <  1:57014/145‑1 (MQ=255)

                                                                                                                      gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                               >  2:57014/1‑145 (MQ=255)

                                                                                                                      gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑AAGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGCTAAAACCTCCTTCACCGTTTAAGAAca                                                >  1:745445/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCATCAAA‑CCCTTCCATACAGCTCAGGTTCCACCAGTGCACCGCCTTTTGC‑CTGTTCTTCTCCATAAAACCTTCTTCACCATTTAAAAACATCCtt                                           >  2:1086307/1‑238 (MQ=25)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCTGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:1232043/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:448139/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:1275122/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:85144/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:103384/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:1017289/1‑238 (MQ=255)

                                                                                                                           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTCCAACCAGGGCAGCGCCTTTTGC‑TTGTTCTTCTCGTTAAAACCTTCTTCACCGTTTAAGAACATCAtt                                           >  1:558876/1‑238 (MQ=255)

                                                                                                                            gcacttgcttGATTGTGCGTGATTGTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑‑CCCATCCTAACTCAATGCCGGAGGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTc                                          <  2:27065/229‑1 (MQ=37)

                                                                                                                            gcAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCCACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat                                               <  1:111806/233‑1 (MQ=255)

                                                                                                                            gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTc                                          <  1:466871/238‑1 (MQ=255)

                                                                                                                            gaaaTAGGCTGATTTAGGGTGATTTTGGAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTc                                          <  1:1084456/236‑1 (MQ=255)

                                                                                                                             cAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATTATTccc                                        <  1:1272598/239‑1 (MQ=255)

                                                                                                                                tagtgtGATTTAGCGTGATTTTTGAGGGCGTCGAAGGCTTTT‑‑ACGC‑‑CGCATCCGACATCAATGCTGGATGCGACGCTTGGCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                     <  2:785993/234‑1 (MQ=255)

                                                                                                                                tAGCCTGGTTTAGCGTTGTGTGTTGGGTCGGAAAAGGCGTTT‑‑ACGC‑‑CGCCTCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                     <  1:780777/239‑1 (MQ=38)

                                                                                                                                tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCTGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATa                                                                                                                          >  1:147492/1‑155 (MQ=255)

                                                                                                                                  gCCTGTTTTAGCGTGTTTTTGTAGGGCGGATAAGGCGTTTCCG‑‑C‑‑CGCATCCGACATCAATGCCTGATTCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCa                                   <  2:1017289/239‑1 (MQ=255)

                                                                                                                                     tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTa                                <  1:538334/239‑1 (MQ=255)

                                                                                                                                      taTTTAGGGTGTTTTTGTAGGTCGGATTAGGCGTTT‑‑ACCC‑‑CGCATCCGACATCAATGCCTGATGCGACTATTGCCGCGTCTTATCAGGCCTATCTTAACCTTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                               <  2:871965/238‑1 (MQ=255)

                                                                                                                                       aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACt                              <  1:1108800/239‑1 (MQ=255)

                                                                                                                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACt                              <  1:782537/238‑1 (MQ=255)

                                                                                                                                        tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTCTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACt                              <  1:592995/238‑1 (MQ=255)

                                                                                                                                            tgtggtttttGTAGGTCGGATAAGGCTTTTTCG‑‑C‑‑CGCACCCTAAATCAATGCCTGAGGCGACGCTTGCCTGGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCCTTCCATACAGCTAAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCGCCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                          <  2:1020593/233‑1 (MQ=255)

                                                                                                                                                gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAAATCATTCCCGCCCAGTACTCGGATTac                      >  1:209801/1‑237 (MQ=255)

                                                                                                                                                gATTTTGGAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc                     <  1:20087/239‑1 (MQ=255)

                                                                                                                                                 aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATaaaa                                                                      >  1:662720/1‑189 (MQ=255)

                                                                                                                                                       tAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg              <  1:938815/239‑1 (MQ=255)

                                                                                                                                                            cGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCCGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGCTTCGGGaaa         >  1:833691/1‑239 (MQ=255)

                                                                                                                                                                   ggCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGCTAAAACCTCCTTCACCGTTTAAGACCTTCATTCCCGCCCAGTAACCGGAATAACCGGTTCGGGAAATTGCaaa  >  2:421176/1‑237 (MQ=39)

||

CGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA‑CCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGC‑TTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGA  >  NC_000913/4296210‑4296609

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