BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,880,794	A→G	3.3%	Y85H (TAC→CAC)	casD ←	CRISP RNA (crRNA) containing Cascade antiviral complex protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,880,794	0	A	G	3.3%	439.4 / 3.3	122	Y85H (TAC→CAC)	casD	CRISP RNA (crRNA) containing Cascade antiviral complex protein
Reads supporting (aligned to +/- strand): ref base A (65/53); new base G (2/2); total (67/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 122 total aligned reads displayed.

TCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAGGATTT  >  NC_000913/2880508‑2881081

  gccgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCACACCACGGTAATCTTCTCGCGCTCCAAGGCACGTACGGTAATc                                                                                                                                                                                                                                                                                               >  2:393875/1‑288 (MQ=255)

  gccgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATc                                                                                                                                                                                                                                                                                               >  2:426005/1‑288 (MQ=255)

  gccgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAAAACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATc                                                                                                                                                                                                                                                                                               >  2:896293/1‑288 (MQ=255)

  gccgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGCTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATCCGCGCCATATTTGAATCGTTTCATGCCTTTTCAAACTACGGTAAGCATCTCGCGCTCCCAGGCCTGTATGGTAATc                                                                                                                                                                                                                                                                                               >  2:759870/1‑288 (MQ=255)

    cgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCCGATAAATATTCGCGCCATGTTTGAATAGTTTCATGCCTTTTACAACCACGCTAACCGTCCCGCACTCCAAGGACTGCATGGTAATcca                                                                                                                                                                                                                                                                                             >  2:530135/1‑286 (MQ=255)

    cgccAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCC‑AAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGc                                                                                                                                                                                                                                                                                           >  1:729100/1‑289 (MQ=255)

    cgccAACGGGCTCATAATTTAACAGCGCCGTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTGTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTGGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCAc                                                                                                                                                                                                                                                                                             >  1:304840/1‑288 (MQ=255)

                              gcCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACacgacg                                                                                                                                                                                                                                                                    <  2:641050/287‑1 (MQ=255)

                               ccTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATc                                                                                                                                                                                                                                                                 >  2:114819/1‑289 (MQ=255)

                               ccTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATc                                                                                                                                                                                                                                                                 >  1:76152/1‑289 (MQ=255)

                               ccTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATc                                                                                                                                                                                                                                                                 >  1:341795/1‑289 (MQ=255)

                                   cTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCa                                                                                                                                                                                                                                                             <  1:89445/289‑1 (MQ=255)

                                              aTGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTc                                                                                                                                                                                                                                                  >  1:800089/1‑289 (MQ=255)

                                                           ccccAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTGTACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTACAGGAGGCATCACATAAATATTCGCGCCATGTTGGAATCGTTACATGACTTTTCAAACAACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCCCGCAACCCCGTTCAAGACACTCGACGATCGTCAAGAATGAGTTCatcgcatcgc                                                                                                                                                                                                                                        >  2:191519/1‑286 (MQ=255)

                                                                     aGCTGGTGTGTTAGTGGGCAACTTCTCCTGCCCAGGTCAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCCTCGTTGCATGGGGTGTTAACCAGGGAGCGACGGTAAAGGGGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACAACGGTAATCTTCTCGCGCTCCCAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACTATCGTCAAGAATGAGTTCATCGCAGCTCACTGCAAATTGca                                                                                                                                                                                                                           <  2:304840/289‑1 (MQ=255)

                                                                     aGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           <  1:318628/289‑1 (MQ=255)

                                                                     aGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           <  2:341795/289‑1 (MQ=255)

                                                                       gaccgggtgtAGTTGGCAACTTCTCCGCACCAAGTAAGGTGGATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGAATGGGGTGTTAACCAGAGCGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           <  2:599761/278‑1 (MQ=255)

                                                                        gggTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         <  1:747215/288‑1 (MQ=255)

                                                                        gggTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCGTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAAGATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         <  1:191519/288‑1 (MQ=255)

                                                                           tgtgtTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACActc                                                                                                                                                                                                                      <  1:247238/288‑1 (MQ=255)

                                                                            gtgtTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACActct                                                                                                                                                                                                                     <  1:113597/288‑1 (MQ=255)

                                                                                tAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGAt                                                                                                                                                                                                                <  2:76152/289‑1 (MQ=255)

                                                                                    gggCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAAAACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATaa                                                                                                                                                                                                              <  1:896293/287‑1 (MQ=255)

                                                                                      gCAACTTCTCCGGCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGc                                                                                                                                                                                                           <  1:596780/288‑1 (MQ=255)

                                                                                      gCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTGTCAAGTTCTTAGATAACCAACGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGc                                                                                                                                                                                                           <  1:572330/288‑1 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACa                                                                                                                                                                                                                                                                               <  2:596993/219‑1 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACa                                                                                                                                                                                                                                                                               >  1:596993/1‑219 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         >  2:338123/1‑273 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         <  1:338123/273‑1 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           >  2:203706/1‑271 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           <  1:203706/271‑1 (MQ=255)

                                                                                       cAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGAtcttca                                                                                                                                                                                                           >  1:599761/1‑282 (MQ=255)

                                                                                               ccGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATgaa                                                                                                                                                                                                 <  2:729100/289‑1 (MQ=255)

                                                                                                cGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATgaag                                                                                                                                                                                                <  2:649657/289‑1 (MQ=255)

                                                                                                cGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATgaag                                                                                                                                                                                                <  1:114819/289‑1 (MQ=255)

                                                                                                       ggTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           >  2:891441/1‑255 (MQ=255)

                                                                                                                 ataCCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTg                                                                                                                                                                                <  1:543613/288‑1 (MQ=255)

                                                                                                                  taCCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGa                                                                                                                                                                              >  1:352504/1‑289 (MQ=255)

                                                                                                                       aGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCa                                                                                                                                                                         >  2:823671/1‑289 (MQ=255)

                                                                                                                       aGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCa                                                                                                                                                                         >  1:603652/1‑289 (MQ=255)

                                                                                                                       aGGCTTTAATACTGCTTTTTCAAGTGCTGAGATAACCATCGTTGCACGGGGTGTTAACCACAGAGCGACGGTAAAGGAGGCATCACATCAATATTCGCGCCATGTGTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCACCCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCCCAGCGCACTGACAATTGCACACCCTCTTATACATCCGGTACCGAACAAATAACCACACGTTGGAACCCa                                                                                                                                                                         >  2:638153/1‑289 (MQ=255)

                                                                                                                       aGGCTCTAATACTGCTTTTTCAACTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGATAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTCGCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGCCTGTATGGTAATCACGCTACTCCGTTACATACACACCACGATCGTCAACAATGAGTTCATCTCAGCGCACTGCAAATCGCACTCTCTCTGATAACGCCTGTACTGAAGAAATATCATTCCCTTCGATCCCa                                                                                                                                                                         >  2:66428/1‑289 (MQ=255)

                                                                                                                        ggCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGTATCCCaa                                                                                                                                                                        >  1:676359/1‑289 (MQ=255)

                                                                                                                        ggCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCaa                                                                                                                                                                        >  1:618740/1‑289 (MQ=255)

                                                                                                                           tttAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGAc                                                                                                                                                                     <  1:393875/289‑1 (MQ=255)

                                                                                                                                tACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGcc                                                                                                                                                                <  1:131623/289‑1 (MQ=255)

                                                                                                                                   tGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAATGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGa                                                                                                                                                                                                                                                                             >  1:435419/1‑177 (MQ=255)

                                                                                                                                    gCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           >  2:706335/1‑226 (MQ=255)

                                                                                                                                    gCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGca                                                                                                                                                                                                                           <  1:706335/226‑1 (MQ=255)

                                                                                                                                    gCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGa                                                                                                                                                            <  2:603652/289‑1 (MQ=255)

                                                                                                                                     cTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCa                                                                                                                                                                                                                                                                                              >  2:499909/1‑158 (MQ=255)

                                                                                                                                     cTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCa                                                                                                                                                                                                                                                                                              <  1:499909/158‑1 (MQ=255)

                                                                                                                                     cTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcac                                                                                                                                                                                                                          >  1:430567/1‑226 (MQ=255)

                                                                                                                                     cTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcac                                                                                                                                                                                                                          <  2:430567/226‑1 (MQ=255)

                                                                                                                                     cTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAg                                                                                                                                                           >  1:531182/1‑289 (MQ=255)

                                                                                                                                      tttttCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAg                                                                                                                                                           <  2:618740/288‑1 (MQ=255)

                                                                                                                                                         aCCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGTCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCtt                                                                                                                                        >  1:534052/1‑288 (MQ=255)

                                                                                                                                                         aCCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTGATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGAc                                                                                                                                                                     <  2:683265/259‑1 (MQ=255)

                                                                                                                                                         aCCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTGATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGAc                                                                                                                                                                     >  1:683265/1‑259 (MQ=255)

                                                                                                                                                                     tGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGTAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATCACCCGCTTCGGGTTGGtaa                                                                                                                             >  1:575919/1‑287 (MQ=255)

                                                                                                                                                                      ggggTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTGATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATc                                                                                                                           <  2:597505/288‑1 (MQ=255)

                                                                                                                                                                        ggTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGAc                                                                                                                                                                     >  1:374770/1‑244 (MQ=255)

                                                                                                                                                                        ggTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGAc                                                                                                                                                                     <  2:374770/244‑1 (MQ=255)

                                                                                                                                                                                          aCGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCtt                                                                                                      >  2:238867/1‑289 (MQ=255)

                                                                                                                                                                                                      gCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACGCTCTCTGATAACGCCTGTAATGAAGAAGTatcat                                                                                                                                                                                       <  1:336136/196‑1 (MQ=255)

                                                                                                                                                                                                      gCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTatcat                                                                                                                                                                                       >  2:336136/1‑196 (MQ=255)

                                                                                                                                                                                                        aTCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTg                                                                                         >  1:268647/1‑288 (MQ=255)

                                                                                                                                                                                                        aTCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGCTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATAATCACTTTGCATCCCAAGCCAAGCCCCGAATTGCCCTAATAACCCGCCTCTGGTACGAAATCTTCCGGCAGGGCGTGCTCCTTCATAGGCCgggcg                                                                                         >  2:77613/1‑285 (MQ=255)

                                                                                                                                                                                                        aTCACATAAATATGCGCGCCATGTTTGAATCGTTTAATGACATTTCAAACCACGGTAATCTTCTAGCGCACCAAGGACTGCATGGTAATCACGCAACCCCGTTACAGACACACGACGACCGTAAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCCAGTCGCCCTAATAACCCGCTTCGGGTCGGAAATGTCACGGTAGGTCGCGTTCCTTCAAAGGTCGGCTg                                                                                         >  2:644225/1‑288 (MQ=255)

                                                                                                                                                                                                          cacaTAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcc                                                                                       >  2:205727/1‑288 (MQ=255)

                                                                                                                                                                                                           acaTACAGATTCGCGCCGTGTTTGAATCGTTTCATGGCTTTTCAAACCACGGTAATCTTCTCGCGCTCCAGTGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGGACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCGTCCGGGAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                     <  1:644225/289‑1 (MQ=255)

                                                                                                                                                                                                                 atatTCGCGCCATGGTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCAACGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGCAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTGCAAAGGTCGGCTGCCCCCAGGCt                                                                               <  1:638153/289‑1 (MQ=255)

                                                                                                                                                                                                                    ttCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGc                                                                                                                                                                                                                                   >  1:815661/1‑138 (MQ=255)

                                                                                                                                                                                                                     tCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCa                                                                           <  1:618381/289‑1 (MQ=255)

                                                                                                                                                                                                                      cgcgCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCAt                                                                          >  2:536557/1‑289 (MQ=255)

                                                                                                                                                                                                                             tGTTTGAATCGTTTCATGTCTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCACGAGTAGCCCTAATAACCCGCTTC‑GGTCGGGAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGGTGGCCGCCCCCAGGCTTACCTTGGCCCa                                                                   >  2:201022/1‑288 (MQ=255)

                                                                                                                                                                                                                                     tCGTTTCATGACTTGTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCAGCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCg                                                           <  1:269550/289‑1 (MQ=255)

                                                                                                                                                                                                                                      cGTTGCATGACGTTTCAAACCACGGTAAGCTTCTCGCGCTCCAAGGACAGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGc                                                          <  2:49465/289‑1 (MQ=255)

                                                                                                                                                                                                                                       gTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACATACACACGACGATCGTTAAGAATGAGTTCATCGCAGCGCCCTGCAACTTGCACACTCTCTGCTCACGCCTGTAATGACCAAGTATCATCCCGTTCGATCCCAAGACAAGCCCCGAGTAGCCCTCATAACCCCCTTCGGGTGCGACATCTTCCGGTAGGTCGCGTCCATTCAACGGTCGGCTGCCCCCCGGCTTCCCTTCGACCAGCAATCCGCa                                                         >  2:269550/1‑289 (MQ=255)

                                                                                                                                                                                                                                          ttaTGACTTTTCAAAACACGGTAATCTTCTCGCGCTCCAAGGACTGTTTGGGAATCACGCAACCCCGTTACAGACAAACGACGATAGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGCTAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGa                                                      <  2:352504/287‑1 (MQ=255)

                                                                                                                                                                                                                                           cATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGAt                                                     <  1:238867/289‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                     <  2:30948/249‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                     >  1:30948/1‑249 (MQ=255)

                                                                                                                                                                                                                                                   ttCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCAAGCAAGCCGCAAGACCAAAAAAg                                             >  2:905251/1‑289 (MQ=255)

                                                                                                                                                                                                                                                           aCGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAACTCTTCCGGCAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGCTCAAATAAGAttcactt                                     >  2:781482/1‑283 (MQ=255)

                                                                                                                                                                                                                                                             ggTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCt                                                                                          >  1:9898/1‑234 (MQ=255)

                                                                                                                                                                                                                                                             ggTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGtt                                  >  2:768479/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                      ctcGCGCTCCAAGGACTATATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGACCAAATAA‑ACCTCCACGTTCACGCCTc                          >  2:120371/1‑288 (MQ=255)

                                                                                                                                                                                                                                                                                aaGGACTGTATGGTAATCACGCAACCCCTTTACAGCCACACTACGATCGTCAATAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTCCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTAc               <  1:530135/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 aGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACg              <  1:651472/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                                      tGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAAccc         <  1:768479/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                                       gTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         <  1:429451/81‑1 (MQ=255)

                                                                                                                                                                                                                                                                                       gTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGcaca                                                                                                                                                                                                                         >  2:429451/1‑81 (MQ=255)

                                                                                                                                                                                                                                                                                       gTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCa        >  1:44779/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                                       gTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCa        >  2:45636/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                                         aTGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGTCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAgg      >  1:216567/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                                         aTGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAg       <  1:823671/288‑1 (MQ=255)

                                                                                                                                                                                                                                                                                             tAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAGGAttt  >  1:868026/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                                             tAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAAACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAGGAttt  >  2:132851/1‑289 (MQ=255)

TCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCT‑CATGTTCACGCCTCGCCATTATTACGAACCCAGGATTT  >  NC_000913/2880508‑2881081

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