BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	9,907	A→C	5.0%	intergenic (+14/+21)	mog → / ← satP	molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	9,907	0	A	C	5.0%	450.6 / 14.9	142	intergenic (+14/+21)	mog/satP	molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Reads supporting (aligned to +/- strand): ref base A (73/59); new base C (4/3); total (77/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-02

Only 100 of 144 total aligned reads displayed.

TGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCAC  >  NC_000913/9635‑10180

tGTACCAACTGCGATCCTCTCGCGTCAGGTGGGCGTGATTCTTAAACAGGCGCTGATCCTTCACTTACCCGGTCAGCCGAAGTCTATTAAAGAGAAGCTGGAAGGGGTGAAGGACGATGAGGGTAACGNTGTGGTACACGGTATTTTTGCCAGCGTACCTTATTTCATTCTGTTGCTGGAAGGGTAATCCGTTGAAACNGCACCGGACGTGGTTGCAGACTTGCTAACGACGAGTGCAAGACACNANGTTAGCGAAAAAAAAAACACACCCGAGc                                                                                                                                                                                                                                                                                   >  3:53593/1‑275 (MQ=255)

   accgatttcgtgccTTTCGCGTGCGGGGGNCGGTTTTCCCAAACAGGCGCTGATCNTTAACTTACCCGGTCATCCGAATTATGTTAAAGAGACGCTTGAAGGGGTGAAGGACGCTGAGGGTAACGTTGTGGTACACNGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCGGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCggg                                                                                                                                                                                                                                                                                <  3:104246/263‑1 (MQ=255)

       aCTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATCCAGACCGAAGAGTGCAAGACGCGACGTTAGACAATAAAAAAATCCCCCCGACCGGGGGGa                                                                                                                                                                                                                                                                            >  4:430682/1‑275 (MQ=255)

             aTCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaa                                                                                                                                                                                                                                                                      <  4:300804/275‑1 (MQ=255)

                                                  cgcTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGCCGCTGGAAGGTGTGAAGGACGATGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTAACGTACTGCATTCAGTTCGTGGAAAGGACATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGGCCGAAGAGTGCAAGAGGACGCGTTAGCGAATAAACAAATGCACCCGAACGGGGGGATC‑CAAAACACATTCGTGGGATT‑CACCAATCGGCAGCACGGtac                                                                                                                                                                                                                               >  4:248984/1‑273 (MQ=255)

                                                                 tACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGGTCCCAAACCA‑ATTAGTGGGATT‑CACCAATCCGCAGAAAGATGCGACAACACTGCTCGt                                                                                                                                                                                                                >  2:455952/1‑275 (MQ=255)

                                                                   cccGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGCCGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACC‑ATTAGTGGGATT‑CACCAATCGGCAGACCGGGGCGACCAAACTGCCCGtt                                                                                                                                                                                                               >  4:284447/1‑274 (MQ=255)

                                                                   cccGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGCCCAAACTTCTCGtt                                                                                                                                                                                                               >  2:507814/1‑274 (MQ=255)

                                                                   cccGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                                                                                                                                                                                               >  2:189723/1‑274 (MQ=255)

                                                                   cccGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCACGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATTACACC‑AGCGGCAGAACGGTGCGACCAAAATGCCCGtt                                                                                                                                                                                                               >  4:469985/1‑274 (MQ=255)

                                                                         ccgCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGAAGCTGAGTGTAACGTGTTGGTAAACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑AACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACt                                                                                                                                                                                                        <  4:403733/273‑1 (MQ=255)

                                                                               aagttgaTTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTAAACGGTATTTTTGCCAGCGTACCGTACTGAATTCAGTTGCTGGAAGGGCCATACGTTGAGACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                                                                                                                                                                                                   <  1:53593/268‑1 (MQ=255)

                                                                                atgttATTAAAGAGGAGCTGGAAGGTGGGAAGGAAGCGGAGGGTAACGTTGTAGTACACGGTATGGTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCAATCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                                                                                                                                                                                                 <  4:201872/271‑1 (MQ=255)

                                                                                      ttAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTTCTCGTTCAGTACTTACCCCATCGCCa                                                                                                                                                                                           >  4:42832/1‑275 (MQ=255)

                                                                                              aCGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt                                                                                                                                                                                   >  2:196251/1‑275 (MQ=255)

                                                                                                                gccGCTGAGGGTAACGTTGTGGGAAACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGAt                                                                                                                                                                  <  4:507813/272‑1 (MQ=255)

                                                                                                                   gCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCCCCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAg                                                                                                                                                               >  4:249436/1‑274 (MQ=255)

                                                                                                                                  agggTAAACGGTATTTTTGCCAGCGTACCGAACTGCTTTCAGTTGCTGGAAGGGCCAAACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCgg                                                                                                                                                 <  4:3303/271‑1 (MQ=255)

                                                                                                                                   tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATCGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCaa                                                                                                                                              >  1:145381/1‑275 (MQ=255)

                                                                                                                                   tgGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCACTCAAGCCGGCaa                                                                                                                                              >  1:98729/1‑275 (MQ=255)

                                                                                                                                      tACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGa                                                                                                                                                                                                                                                                            <  7:207938/148‑1 (MQ=255)

                                                                                                                                            gTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa‑a                                                                                                                                                                                                                                                                <  4:64940/153‑1 (MQ=255)

                                                                                                                                               tttttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGcatgctg                                                                                                                                   >  2:3303/1‑271 (MQ=255)

                                                                                                                                                ttttGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTgc                                                                                                                                 >  2:402584/1‑275 (MQ=255)

                                                                                                                                                       aGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGNGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTNACCCATCGCCAGATAGATTGCGATGGCCCCGCAGATCCGCCCACCNCNCCCGGCAAAGTGGATGATTgcgggggt                                                                                                                           >  3:51791/1‑270 (MQ=255)

                                                                                                                                                         cgtaccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGTCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCACCGAGCGGGGGGATCTCAAAACA‑ATTAGTTGGATT‑CACCAGTCGCCAGAACGGTGCGACCAACATACTCGTCCAGTACTTCNCACAGCGCAGGATAGATGGAGCTGGCCCCGCAGACCAGCCACATCCAGCCGGCCCAGTGGACGAATTCGACGTCAca                                                                                                                        >  3:122738/1‑274 (MQ=255)

                                                                                                                                                          gtaccgtacTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTAcc                                                                                                                        <  4:189722/274‑1 (MQ=255)

                                                                                                                                                                 acTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGt                                                                                                                <  2:249437/275‑1 (MQ=255)

                                                                                                                                                                   tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGAGACGTTAGCGCATAAAAAAATCCCCCCGAGCGGGGGGATCTCAACACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTAGTTCAGTACTTCACCCATCGCCAGATAGATTGAGCTGGCACCGCCACTAAGCCCACTCCAGCCGGACAAGAGGATGCTTGCGGCGTTACCGGaaacatta                                                                                                              >  4:14337/1‑270 (MQ=255)

                                                                                                                                                                     cATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCACAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAACGTGGATGCGTGCGGGGTTACCGGCAATGTTAcc                                                                                                            >  4:149505/1‑275 (MQ=255)

                                                                                                                                                                     cATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTgg                                                                                                                                                                                                                                                         >  2:180959/1‑135 (MQ=255)

                                                                                                                                                                       ttCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCCCCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGa                                                                                                          >  2:376195/1‑275 (MQ=255)

                                                                                                                                                                          aGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                                                                                                                                                                                                       >  4:169658/1‑179 (MQ=255)

                                                                                                                                                                          aGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCg                                                                                                       >  2:308072/1‑275 (MQ=255)

                                                                                                                                                                              gCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAcc                                                                                                                                                                                                   <  2:169659/179‑1 (MQ=255)

                                                                                                                                                                               cTGGAAGGGCCATAGGTTGAAACGGCACCGGAAGTGGTTGCAGCATTAAGACCGACGAGTGAAAAACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGTGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCAAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCcagc                                                                                                  <  1:96871/275‑1 (MQ=255)

                                                                                                                                                                               cTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAACCTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCcagc                                                                                                  >  2:476453/1‑275 (MQ=255)

                                                                                                                                                                                  gAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCTACGTTAGCGAATAAAAAAATCCCCCCGAGCgggg                                                                                                                                                                                                                                                                               <  6:1059148/101‑1 (MQ=255)

                                                                                                                                                                                           tACGTTTATAGGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAAAGGTCAAAACGCGACGTTAGCGCAGAAAAAAATCCCCCCGAGCGGGTGAATCTTCAAACA‑CGTAGTAGGATG‑CACCAGTCGGAAGAACGGTGCGGCCAAACTGCTCATTCAGTACTTCACCCATCGCCAGATAGATTTCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGAGTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGc                                                                                      <  4:473379/275‑1 (MQ=255)

                                                                                                                                                                                                tGAAACGGCACCGGAAGTGGTTGCAGCACTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACTCATCGCCAGCTAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGt                                                                                 >  1:37114/1‑275 (MQ=255)

                                                                                                                                                                                                 gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCaaa                                                                                                                                                                                                                        >  4:240304/1‑139 (MQ=255)

                                                                                                                                                                                                 gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACCGCACGGTc                                                                                >  4:505633/1‑275 (MQ=255)

                                                                                                                                                                                                 gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGt                                                                                 <  2:42833/274‑1 (MQ=255)

                                                                                                                                                                                                 gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGt                                                                                 <  2:320005/274‑1 (MQ=255)

                                                                                                                                                                                                 gAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGAACGGTc                                                                                >  2:528570/1‑275 (MQ=255)

                                                                                                                                                                                                  aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                                                                                                                                                                                                          <  7:61419/90‑1 (MQ=255)

                                                                                                                                                                                                  aaaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                                                                                                                                                                                                          >  8:61419/1‑90 (MQ=255)

                                                                                                                                                                                                   aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                                                                                                                                 >  6:200008/1‑82 (MQ=255)

                                                                                                                                                                                                   aaCGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                                                                                                                                 <  5:200008/82‑1 (MQ=255)

                                                                                                                                                                                                     cGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGc                                                                                                                                                                                                                    <  2:240305/139‑1 (MQ=255)

                                                                                                                                                                                                                    ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCGCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                                                                                                                                                                                            >  6:917319/1‑148 (MQ=255)

                                                                                                                                                                                                                    ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                                                                                                                                 <  5:423319/65‑1 (MQ=255)

                                                                                                                                                                                                                    ttGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCgg                                                                                                                                                                                                                                                                                 >  6:423319/1‑65 (MQ=255)

                                                                                                                                                                                                                         gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑Ca                                                                                                                                                                                                                                                  >  5:107007/1‑89 (MQ=255)

                                                                                                                                                                                                                         gcaTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CAc                                                                                                                                                                                                                                                 <  5:741785/90‑1 (MQ=255)

                                                                                                                                                                                                                          caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCg                                                                                                                                                                                                                              >  8:1115412/1‑108 (MQ=255)

                                                                                                                                                                                                                          caTTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCg                                                                                                                                                                                                                              <  7:1115412/108‑1 (MQ=255)

                                                                                                                                                                                                                            ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa‑a                                                                                                                                                                                                                                                                >  5:51093/1‑73 (MQ=255)

                                                                                                                                                                                                                            ttCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa‑a                                                                                                                                                                                                                                                                <  6:51093/73‑1 (MQ=255)

                                                                                                                                                                                                                                 aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTg                                                                                                                                                                                                                                                          <  6:501964/74‑1 (MQ=255)

                                                                                                                                                                                                                                 aCCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTg                                                                                                                                                                                                                                                          >  5:501964/1‑74 (MQ=255)

                                                                                                                                                                                                                                    gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                                                                                                                                                           >  8:923084/1‑117 (MQ=255)

                                                                                                                                                                                                                                    gAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                                                                                                                                                                                           <  7:923084/117‑1 (MQ=255)

                                                                                                                                                                                                                                      agagTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCAGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAACACGAATTGCAGAAACCGTGCGCCtt                                           >  4:165066/1‑275 (MQ=255)

                                                                                                                                                                                                                                                aCGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTg                                   <  2:149506/273‑1 (MQ=255)

                                                                                                                                                                                                                                                 cgcgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATATCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAg                                                                                                                                                               <  5:917319/148‑1 (MQ=255)

                                                                                                                                                                                                                                                  ccaaagagacggAAAAAAAAAACCCCCCCNGGGGGGGGGATCACAAAACA‑ATTCANGGGTTT‑CACCAATCCGCAGACGGGTGCGACCAAACTTCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCANCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCga                               <  3:98729/264‑1 (MQ=255)

                                                                                                                                                                                                                                                   cgACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGAGTGCGCCTTTCAGCGTGCCgaa                              >  2:445398/1‑275 (MQ=255)

                                                                                                                                                                                                                                                   caaCGTTAGCGAATAACAAAACCCCCCCGAGCGGGGGGATCTCAAATCA‑CTTAGTCGGATT‑CACCAATCGTCAGAACGGTGCTACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCga                               <  2:537569/272‑1 (MQ=255)

                                                                                                                                                                                                                                                    gACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCCAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCACACAGCa                                                                                     >  2:358551/1‑219 (MQ=255)

                                                                                                                                                                                                                                                           gCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATa                       <  2:430683/274‑1 (MQ=255)

                                                                                                                                                                                                                                                            cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGGTT‑CACCAATCGGCAGAACGGTGc                                                                                                                                                                                                                               <  5:766294/73‑1 (MQ=255)

                                                                                                                                                                                                                                                            cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGGTT‑CACCAATCGGCAGAACGGTGc                                                                                                                                                                                                                               >  6:766294/1‑73 (MQ=255)

                                                                                                                                                                                                                                                            cGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATaaa                     <  2:469986/275‑1 (MQ=255)

                                                                                                                                                                                                                                                               aTAAAAAAATCCCCCCCAGCGGGGGGGTATCAAAACA‑ATTAGGGGGATT‑CACCAAACGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTGACCCGTCGCCAGATAGATTGCGCTGGCACCCCGGATCAGCCCAATCCAGCCGGCAAAGTGGAGGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATaaa                     <  4:469845/272‑1 (MQ=255)

                                                                                                                                                                                                                                                                    aaaaTCCCCCCGGTGGGGGGGGTCTAAAAACA‑ATTATTGGGATT‑CACCCAGCGGCAGAACGGTGCAAGCAAACTGCCCGTCAAGTACTTCGCCCGTCGCCAGATGGATGTGGCTGGCACCGCAGATCGGGCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTaa             <  4:402583/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                       aTCCCCCCGCTCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                                                                                                                                                                                            >  8:428207/1‑81 (MQ=255)

                                                                                                                                                                                                                                                                       aTCCCCCCGCTCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                                                                                                                                                                                            <  7:428207/81‑1 (MQ=255)

                                                                                                                                                                                                                                                                        tCCCCCCGCTCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcg                                                                                                                                                                                <  6:967576/108‑1 (MQ=255)

                                                                                                                                                                                                                                                                        tCCCCCCGCTCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTgcg                                                                                                                                                                                >  5:967576/1‑108 (MQ=255)

                                                                                                                                                                                                                                                                         ccccccGAGCGGGGGGATCTCAAAACA‑ATTAGTTGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCGAt                                                                                                                                                                                                >  7:585062/1‑91 (MQ=255)

                                                                                                                                                                                                                                                                         ccccccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                                                                                                                                                                                                <  8:585062/91‑1 (MQ=255)

                                                                                                                                                                                                                                                                          cccccGCTCGGGGGGATTTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCaa                                                                                                                                                          >  8:168915/1‑128 (MQ=255)

                                                                                                                                                                                                                                                                           ccccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg                                                                                                                                                                                  >  5:288885/1‑103 (MQ=255)

                                                                                                                                                                                                                                                                           ccccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTg                                                                                                                                                                                  <  6:288885/103‑1 (MQ=255)

                                                                                                                                                                                                                                                                           ccccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACg        <  2:317598/273‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTgg                                                                                                                                         <  6:675332/143‑1 (MQ=255)

                                                                                                                                                                                                                                                                             ccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                                                                                                                                                                                            >  7:128991/1‑91 (MQ=255)

                                                                                                                                                                                                                                                                             ccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                                                                                                                                                                                            <  8:128991/91‑1 (MQ=255)

                                                                                                                                                                                                                                                                             ccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatagat                                                                                                                                                                                    <  7:424174/99‑1 (MQ=255)

                                                                                                                                                                                                                                                                             ccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCagatagat                                                                                                                                                                                    >  8:424174/1‑99 (MQ=255)

                                                                                                                                                                                                                                                                             ccGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAg                                                                                                                                            <  5:868431/139‑1 (MQ=255)

                                                                                                                                                                                                                                                                              cGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                                                                                                                                                                                               <  6:776447/71‑1 (MQ=255)

                                                                                                                                                                                                                                                                              cGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                                                                                                                                                                                                      <  5:1004806/80‑1 (MQ=255)

                                                                                                                                                                                                                                                                              cGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                                                                                                                                                                                                      >  6:1004806/1‑80 (MQ=255)

                                                                                                                                                                                                                                                                               gAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAg                                                                                                                                                               <  5:1107395/118‑1 (MQ=255)

                                                                                                                                                                                                                                                                               gAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAg                                                                                                                                                               >  6:1107395/1‑118 (MQ=255)

                                                                                                                                                                                                                                                                               gAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGga                                                                                                                                        >  7:877837/1‑141 (MQ=255)

                                                                                                                                                                                                                                                                               gAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGga                                                                                                                                        <  8:877837/141‑1 (MQ=255)

                                                                                                                                                                                                                                                                                aGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCcac  >  4:37243/1‑274 (MQ=255)

TGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACA‑ATTAGTGGGATT‑CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCAC  >  NC_000913/9635‑10180

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