BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	374,345	T→G	7.9%	intergenic (+240/‑338)	dmpG → / → Y75_RS01815	4‑hydroxy‑2‑oxovalerate aldolase/3‑(3‑hydroxy‑phenyl)propionate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	374,345	0	T	G	7.9%	287.5 / 23.9	118	intergenic (+240/‑338)	dmpG/Y75_RS01815	4‑hydroxy‑2‑oxovalerate aldolase/3‑(3‑hydroxy‑phenyl)propionate transporter
Reads supporting (aligned to +/- strand): ref base T (56/51); new base G (4/6); total (61/57)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 118 total aligned reads displayed.

TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGT  >  NC_007779/374109‑374580

tgaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGGAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCCCGGGGATTCGTAGGCCTGATAAGACGCG‑TTAGCAGCGCATCAGGGATCTgc                                                                                                                                                                                                                                             >  1:2108715/1‑239 (MQ=255)

 gaCTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTgc                                                                                                                                                                                                                                             >  1:2087735/1‑238 (MQ=255)

    tGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCa                                                                                                                                                                                                                                          >  2:1963527/1‑238 (MQ=255)

                gCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGACGCG‑GTAGCGGCGCATCAGGCATCTGCGCACGACTGCCgggt                                                                                                                                                                                                                              >  2:347974/1‑236 (MQ=255)

                  ttgttgtcggCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCGTTAGCCGGGCTAAGGATGGCGTGAGAATTTGTAGGTCTGATAAGAGGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgc                                                                                                                                                                                                                            <  2:62423/229‑1 (MQ=14)

                       gTATGCCCGTCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGGTAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTg                                                                                                                                                                                                                      <  1:329324/239‑1 (MQ=255)

                         aTGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGa                                                                                                                                                                                                                     <  1:498970/238‑1 (MQ=255)

                             ccATCCCAGCACGACGGCAGGATGGGGCGGGAACGCCGTATCCGGCCGACGGATGGCGTGAGAATTTGTAGGTCTGAGAAGACGCGTTAGCGTCGCATCGGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGc                                                                                                                                                                                                                 <  2:1975368/238‑1 (MQ=14)

                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAt                                                                                                                                                                                                                              <  2:1374490/222‑1 (MQ=40)

                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt                                                                                                                                                                                                               >  1:693897/1‑237 (MQ=40)

                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                             >  1:532462/1‑239 (MQ=40)

                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                             >  1:1726945/1‑239 (MQ=40)

                                tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCCGGCCTCTGCGCACGACTGCCGGATGCGGCGGTAAGGCCTTa                                                                                                                                                                                                             >  1:2188759/1‑239 (MQ=255)

                                  gcaggccggcTGGCGGATGGGTCGGGAACGCCTTATCGGGCCTGCGGATGGCGGGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                                                                           <  2:375501/230‑1 (MQ=11)

                                       ACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                                                                      >  2:303911‑M1/1‑239 (MQ=255)

                                                ggATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGGTAAGACGCG‑TTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGGTGCGGTGTGAACCCCTTCTCCCGGCTACgggggg                                                                                                                                                                                             >  2:680242/1‑239 (MQ=14)

                                                        cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGTATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGGAAGCCTGATAAGACGCG‑TTAGCGTCGCCTCAGGGATTCGCGCACGACTGCCGGATGCGGCGTTCACGCCCTCTGCCGCCCACggggcgcgggggga                                                                                                                                                                                     >  2:675795/1‑228 (MQ=11)

                                                        cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑CTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCCGCGTGAACCCCCTTACCCGCCTACGGGTGCCGCGGGaa                                                                                                                                                                                     >  2:1379549/1‑239 (MQ=11)

                                                        cgTGAACACCTGTTCCGGCCTACGGCTGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCgaga                                                                                                                                                                                      <  2:1338770/238‑1 (MQ=255)

                                                             aCGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAGCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTg                                                                                                                                                                                 <  2:619956/238‑1 (MQ=255)

                                                              cGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtag                                                                                                                                                                               <  1:428750‑M1/239‑26 (MQ=255)

                                                                     tCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGa                                                                                                                                                                        >  1:2270113/1‑239 (MQ=255)

                                                                          ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGCATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGGGGCGCGCGGGGATTTTTAGGCCTGGTAAGa                                                                                                                                                                   >  2:1418727/1‑239 (MQ=11)

                                                                          ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGGTAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCa                                                                                                                                                                                                                                          >  1:610791/1‑168 (MQ=255)

                                                                          ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAg                                                                                                                                                                    <  2:532462/238‑1 (MQ=255)

                                                                           gtACGGAGGGCGTGATAGCTTGTAGGTCTGATAAGACGCGGTGGCGTCGCATCAGGGATCTGGGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGACTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                   <  2:1028076/237‑1 (MQ=14)

                                                                           cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                   <  1:1913680/238‑1 (MQ=255)

                                                                           cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                   <  1:726354/238‑1 (MQ=255)

                                                                           cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGGGAATTTTGAGGCCCGATAAGAc                                                                                                                                                                  >  1:1438024/1‑239 (MQ=255)

                                                                              cGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcctacacttcgcccgtaagccgtaggcctgataagacgcgc                                                                                                                                                               >  1:1274921‑M1/1‑198 (MQ=255)

                                                                               ggATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑t                                                                                                                                                             <  1:347974/239‑1 (MQ=255)

                                                                                gATGGCGGGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑t                                                                                                                                                             <  2:2097931/238‑1 (MQ=255)

                                                                                 aTGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑tt                                                                                                                                                            <  2:333997/238‑1 (MQ=255)

                                                                                   ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGc                                                                                                                                                         >  1:189725/1‑239 (MQ=255)

                                                                                    gCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCg                                                                                                                                                        >  1:755642/1‑239 (MQ=255)

                                                                                      gcgggaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                            >  2:10924/5‑186 (MQ=14)

                                                                                          gaATTTGTAGGTCTGAGAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCAGGACTGCCGGATGGGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑a                                                                                                                                                  <  2:1949466/238‑1 (MQ=255)

                                                                                          gaATTTGGAGGGCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑a                                                                                                                                                  <  2:1993685/238‑1 (MQ=255)

                                                                                           aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACggat                                                                                                                                                                                               >  1:628129/1‑192 (MQ=25)

                                                                                           aaTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACggat                                                                                                                                                                                               <  2:628129/194‑3 (MQ=25)

                                                                                                     tCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAAGCG‑TTAGCGTCGCATCCGGCATCTGCGCACCACTGCCGGACGCGGCGTGAACCCCTTATCCGGCCCACGGGTGGCGCCAGAATTTGTAGGCCCGATCAGACGCG‑TTACCCGCGC‑ATCAGGCagccg                                                                                                                                       >  2:1837418/1‑234 (MQ=11)

                                                                                                      cTGGTAAGAGGCGTTGGCGTCGCGTCAGACATCTGGGCGCGACTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGGATGGCGCGGGAATGTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTgcg                                                                                                                                     <  2:1698582/239‑1 (MQ=11)

                                                                                                       tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGCCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGGGCGGCGCGGGAATTTGGAGGCCTGATAAAACGCG‑TTAGCG‑CGCGATCGGGGATCTgcgc                                                                                                                                    >  2:1602035/1‑239 (MQ=11)

                                                                                                       tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTgcgc                                                                                                                                    >  1:1400364/1‑239 (MQ=38)

                                                                                                                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCCTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACGGGTGgcg                                                                                                                                                                                           >  2:839766/1‑174 (MQ=255)

                                                                                                                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGCACGCCTTATCCGCCCTACGGGGGGCGCGGGAATTTGTAGGCCTGTTAAAACGCG‑TTAGCCTCCC‑CTCCGGCACCCGCGCACGACCGGCGGa                                                                                                                        >  2:1557103/1‑239 (MQ=11)

                                                                                                                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcg                                                                                                                                                                                           <  1:839766/174‑1 (MQ=38)

                                                                                                                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGa                                                                                                                        >  1:1861660/1‑239 (MQ=38)

                                                                                                                       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATgc                                                                                                                     <  2:189725/238‑1 (MQ=38)

                                                                                                                         tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATg                                                                                                                                                                                                                             >  2:1943491/1‑134 (MQ=34)

                                                                                                                         tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATaa                                                                                                                                                                     >  1:885895/1‑190 (MQ=39)

                                                                                                                         tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATaa                                                                                                                                                                     <  2:885895/190‑1 (MQ=39)

                                                                                                                             aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                               <  2:693897/238‑1 (MQ=40)

                                                                                                                             aTCAGGCACCTGCGCGAGACGGCCGGATGCGGCGGGAACGCGTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTa                                                                                                               <  2:1313503/238‑1 (MQ=255)

                                                                                                                                         cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGCGAACGCCTTATCCGGCCTCCGGGTGGCGCGAGCATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCCC‑ATCCAGCCACTGCGTCCGCCCGGCGGATGCGGGGCAAACGCCttttgcc                                                                                                  >  2:463211/1‑234 (MQ=14)

                                                                                                                                           cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGGAACGCCTTTTCCGGCCTAAGGGTGGCGCGAGAATTTGTAGGCCTGCTAAGACGCG‑GTATCGTCGC‑ATCAGGCATCTGCGCACGACTTCCGGATGCGGCGGAAACGCCTTATCCgg                                                                                                 >  2:2085496/1‑238 (MQ=11)

                                                                                                                                              gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGCAAAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataa                                                                                                                                                                      >  2:346087‑M1/1‑134 (MQ=255)

                                                                                                                                               aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtcttcacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     <  1:1394597‑M1/167‑36 (MQ=255)

                                                                                                                                               aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     >  2:1394597‑M1/1‑132 (MQ=255)

                                                                                                                                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGCAAAGCGTCGC‑ATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGTCTAcacttcg                                                                                      <  2:1815432/238‑7 (MQ=11)

                                                                                                                                                        tGCGGCGTGAGCGCGTTATCCGGCTTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCGCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGTGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                    <  1:2236080/238‑1 (MQ=255)

                                                                                                                                                        tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     >  1:101378‑M1/1‑123 (MQ=255)

                                                                                                                                                        tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     <  2:101378‑M1/158‑36 (MQ=255)

                                                                                                                                                        tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcgc                                                                                   >  1:1087714/1‑239 (MQ=255)

                                                                                                                                                                gAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGcctacacttcgcccgtaagccgtaggcctgataagacgcgcaaagcgtcgcatcaggcatggcgcacgactgccggatgcggcgtgaacgccttatccggtctacacttcgcccataaaccg                                                                              <  1:119301‑M1/238‑123 (MQ=255)

                                                                                                                                                                     ccTTGTCCGGCTTACGGAGCGCGCGGGAATTTGTAGGCGTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                       <  2:1316283/238‑1 (MQ=255)

                                                                                                                                                                      cgcatcaGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataa                                                                                                                                                                      >  2:45790‑M1/4‑110 (MQ=255)

                                                                                                                                                                      cgcatcaGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataa                                                                                                                                                                      <  1:45790‑M1/141‑35 (MQ=255)

                                                                                                                                                                              gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCCTGTCCGCGCTACGGATGGCGCGGGAATTTGTAGGCCTGACAAg                                                             >  2:1922140/1‑239 (MQ=255)

                                                                                                                                                                                ctacggggtgcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCGGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                            <  2:2270113/229‑1 (MQ=255)

                                                                                                                                                                                cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGa                                                            <  1:197703/238‑1 (MQ=255)

                                                                                                                                                                                 tACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGtg                                                                                                               <  1:10924/186‑2 (MQ=38)

                                                                                                                                                                                   cGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGACGGCGCGGGAATTTGTAGGCCTGATAAAAcccg                                                        >  2:827243/1‑239 (MQ=255)

                                                                                                                                                                                   cGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                        >  1:1288665/1‑239 (MQ=255)

                                                                                                                                                                                    ggATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑CTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGATGGCCCGGGAATTTGGAGGCCCGGTAAAAcggcg                                                       >  2:1878070/1‑236 (MQ=14)

                                                                                                                                                                                    ggATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TCAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGCACTTGTAGCCACGCTAAGAcgcct                                                       >  2:2268981/1‑237 (MQ=14)

                                                                                                                                                                                         gcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     >  1:36792‑M1/1‑91 (MQ=255)

                                                                                                                                                                                         gcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataag                                                                                                                                                                     <  2:36792‑M1/126‑36 (MQ=255)

                                                                                                                                                                                         gcgcGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACCCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACCGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAAGGCCCGGGCATTTGGACGCCTGGTAAGACGGGtttacc                                                  >  2:2102134/1‑235 (MQ=14)

                                                                                                                                                                                             gggAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataa                                                                                                                                                                      >  2:929087‑M1/1‑86 (MQ=255)

                                                                                                                                                                                             gggAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCAT‑GGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGtctacacttcgcccataaaccgtaggcctgataa                                                                                                                                                                      <  1:929087‑M1/120‑35 (MQ=255)

                                                                                                                                                                                               gAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcctacggatggcgcgggaatttgtaggcctgataagacgcgcaaagcgtcgcatcaggcatggcgccgactgccggatgcggcgtgaacgccttatccggtctacacttcgcccataaaccgtaggcctgataagacgcgcatagcgtcgcat                                               <  1:250567‑M1/238‑154 (MQ=255)

                                                                                                                                                                                                       tggctgtagaaGGCGCG‑TTAGCGTCGGCACGGGCATCTGCGCACTACGGCCGGATGCGGCGTGAACGCCTTTTCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                     <  2:1516532/229‑1 (MQ=11)

                                                                                                                                                                                                       aGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                     <  2:1861660/238‑1 (MQ=38)

                                                                                                                                                                                                         gCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                    <  1:159913/173‑1 (MQ=37)

                                                                                                                                                                                                         gCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGATGGCGCGGGGATTTGTAGGCCTGATAAAACGCGTTCGCGTCCCACCTGGCATcc                                   >  2:2044225/1‑237 (MQ=255)

                                                                                                                                                                                                         gCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCCTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                    >  2:159913/1‑173 (MQ=37)

                                                                                                                                                                                                            tataagagccag‑TTAGCGTCGCATCAGGCATCTGCGCACTACTGCCGTATGCGGCGTGAACGCATTCACCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGca                                <  1:1970622/227‑1 (MQ=11)

                                                                                                                                                                                                                aaGACGCGCATAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACg                            <  1:1512914/239‑1 (MQ=255)

                                                                                                                                                                                                                aaGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGa                           <  1:2104379/239‑1 (MQ=38)

                                                                                                                                                                                                                    cgcg‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTACCCGGCCTACGGACGGCGCGGGCATTTGTAGGCCTGATAAGACGCCTCAGCGTCGCATCTGGCATCTGCAcaccacct                        >  2:1185606/1‑236 (MQ=255)

                                                                                                                                                                                                                    cgcg‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCCTCGCCTCAGGCATCTGCACACGACTg                        >  1:1755843/1‑238 (MQ=255)

                                                                                                                                                                                                                       g‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCgg                    <  1:58651/239‑1 (MQ=38)

                                                                                                                                                                                                                         ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGTCTACGTACGGCGCGGGAATTTTTAGGCCCGATAAAGCGCGGTTGCGTCGCATCAGGCAGCGGCCCACGACTGccctg                   >  2:1234165/1‑236 (MQ=11)

                                                                                                                                                                                                                          tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGGCGCGGGAATTTGCAGGCCTGATAAGACGCGGTAGCCTCGCCACCACCATCCGCACACGACCGCCGGac                  >  2:951471/1‑238 (MQ=11)

                                                                                                                                                                                                                           aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATg                 <  1:1421407/239‑1 (MQ=25)

                                                                                                                                                                                                                             cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCCCATCAGGCATCTGCACACGACTGCCGGATGCg               >  1:858616/1‑239 (MQ=255)

                                                                                                                                                                                                                                 gCATCAGGCATCTGCGCACGGCTGCCGGATGCGGCGTGAGCGCCTTAGCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa           <  1:951471/239‑1 (MQ=255)

                                                                                                                                                                                                                                  cATCAGGAATCTGCGCACGACTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACGGGGGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa           <  2:2188759/238‑1 (MQ=255)

                                                                                                                                                                                                                                          aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt  <  1:2102134/239‑1 (MQ=255)

TGACTGCCGAGAACGTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCG‑TTAGCGTCGC‑ATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGT  >  NC_007779/374109‑374580

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: `‑`