BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	254,228	A→T	5.2%	intergenic (+26/+31)	prfH → / ← pepD	pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	254,228	0	A	T	5.2%	408.9 / 9.1	115	intergenic (+26/+31)	prfH/pepD	pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand): ref base A (46/63); new base T (3/3); total (49/66)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

Only 100 of 127 total aligned reads displayed.

CGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGA  >  NC_000913/253968‑254502

cGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCg                                                                                                                                                                                                                                                                      >  2:68920/1‑275 (MQ=255)

  aCTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                                                                                                                                                                                                    >  2:69156/1‑275 (MQ=255)

   cTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgca                                                                                                                                                                                                                                                                   <  2:334307/275‑1 (MQ=255)

           tACGCGCCACTCATTTGGCATCCGGTATTAGCTGGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatca                                                                                                                                                                                                                                                           <  1:20222/275‑1 (MQ=255)

                                            tGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCAGAAGGCGATTAAACGCACTCCGGACGCAAATCAAAGAATTACTTCGCCGGAATTTCTTTcagcag                                                                                                                                                                                                                           >  4:141615/1‑274 (MQ=255)

                                              aaGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGACTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAATGGCCGGAAGGCGATTAAACGCCATCAGGCAGCAAATCAAAGAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                                                                                                                                                                                                        >  4:413455/1‑275 (MQ=255)

                                               aGGTTCAGTCAGGGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                                                                                                                                                                                                                          <  1:75238/272‑1 (MQ=255)

                                                    cagtcagAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGGGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag                                                                                                                                                                                                                  <  2:4461/275‑1 (MQ=255)

                                                      tgcagGGCGTAGTTAGGATGCGACACAGCGGCTGGCTCGATTGCTGTTTGCCTGGAAGCTGGAGCCACAGCGACAGGAAAATAGCGCGGCGCTGAAAGCGCAGCGGCGAATGGTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGAGTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                                                                                                                                                                                                                <  4:278185/273‑1 (MQ=255)

                                                           agCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTCTCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCCGCAAATCAAATACTTACTTCGACGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                                                                                                                                                                                           >  2:10263/1‑275 (MQ=255)

                                                           agCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                                                                                                                                                                                           >  2:206008/1‑275 (MQ=255)

                                                            gCGTAGTCAGCAAGCAAACAAAGGGGTGGCACGAATGCGGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAGCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                                                                                                                                                                                           <  4:333426/274‑1 (MQ=255)

                                                            gCGTAGTCAGAATGCTAACAAGGGGCTGGAACGATTGCTGATTGCCTGGAAGCTGGGGCAACAGCAACAGGAAAATAGAGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                                                                                                                                                                                           <  4:163909/274‑1 (MQ=255)

                                                                 gTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTaa                                                                                                                                                                                                     >  4:6315/1‑275 (MQ=255)

                                                                 gTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGATTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTAAGTCAGCAGTTGCCAGTa                                                                                                                                                                                                      >  4:93657/1‑274 (MQ=255)

                                                                          cTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCTTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa                                                                                                                                                                                              <  2:64384/273‑1 (MQ=255)

                                                                                     cTGGCACGATTGCTGATTGCCTGGAAGCTGGTGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATAACCAGATTGACCGTGGCAACCCGCGACGGACATTCAAAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCAAGTAATGACCTACGCTTTCGAtgta                                                                                                                                                                                 >  4:291259/1‑274 (MQ=255)

                                                                                       ggCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGa                                                                                                                                                                                >  4:72207/1‑274 (MQ=255)

                                                                                           cGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTg                                                                                                                                                                           <  2:309136/275‑1 (MQ=255)

                                                                                             gttGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCt                                                                                                                                                                         <  2:244035/274‑1 (MQ=255)

                                                                                             aTTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGa                                                                                                                                                                                                  >  4:139427/1‑250 (MQ=255)

                                                                                             aTTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCt                                                                                                                                                                         <  2:141616/275‑1 (MQ=255)

                                                                                                 cTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa                                                                                                                                                                                              <  2:139428/250‑1 (MQ=255)

                                                                                                      tggcgggAAGCTTGAGCAACAGCAACGNNAAGATCGCGCGGCGCTGAAATCGCAGCGGCNAATGTTCCATCCGCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATNAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCgg                                                                                                                                                                  <  3:67560/268‑1 (MQ=255)

                                                                                                      tGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGga                                                                                                                                                                 <  2:6315/274‑1 (MQ=255)

                                                                                                        ccTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGagag                                                                                                                                                              >  4:345490/1‑275 (MQ=255)

                                                                                                              aGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGTAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGa                                                                                                                                                        <  2:93658/275‑1 (MQ=255)

                                                                                                                cTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGTCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGTGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGAc                                                                                                                                                       <  1:128359/274‑1 (MQ=255)

                                                                                                                   gAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATCAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGt                                                                                                                                                   >  2:254615/1‑275 (MQ=255)

                                                                                                                   gAGCAACAACAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAACGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGt                                                                                                                                                   <  2:444388/275‑1 (MQ=255)

                                                                                                                     gcaacagcaacaGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa                                                                                                                                                                                                                                                                                <  8:642316/148‑1 (MQ=255)

                                                                                                                      caacagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                                                                                                                                                                                               <  7:342790/148‑1 (MQ=255)

                                                                                                                       aacagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc                                                                                                                                                                                                                                                                              <  5:923149/148‑1 (MQ=255)

                                                                                                                         cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                                                                                                                                                                                                    >  7:577874/1‑140 (MQ=255)

                                                                                                                         cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                                                                                                                                                                                                    <  8:577874/140‑1 (MQ=255)

                                                                                                                         cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa                                                                                                                                                                                                                                                                                >  8:108475/1‑144 (MQ=255)

                                                                                                                         cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGt                                                                                                                                             >  2:343362/1‑275 (MQ=255)

                                                                                                                              acagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                                                                                                                                                                                                        >  8:331422/1‑147 (MQ=255)

                                                                                                                              acagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                                                                                                                                                                                                                                                                        <  7:331422/147‑1 (MQ=255)

                                                                                                                              acagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGcc                                                                                                                                        <  4:68919/275‑1 (MQ=255)

                                                                                                                               cagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGcc                                                                                                                                        <  2:345491/274‑1 (MQ=255)

                                                                                                                                  gAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                                                                                                                                                                                                    <  8:780692/147‑1 (MQ=255)

                                                                                                                                  gAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc                                                                                                                                                                                                                                                                    >  7:780692/1‑147 (MQ=255)

                                                                                                                                          cgcgGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCa                                                                                                                            >  1:123305/1‑275 (MQ=255)

                                                                                                                                            cggcgCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                                                                                                           >  2:297765/1‑274 (MQ=255)

                                                                                                                                            cggcgCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                                                                                                           >  2:504471/1‑274 (MQ=255)

                                                                                                                                               cgcTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGCTAGAAACCATGTcc                                                                                                                       >  1:48311/1‑275 (MQ=255)

                                                                                                                                                 cTGGAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTTGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCa                                                                                                                      <  2:17909/274‑1 (MQ=255)

                                                                                                                                                    aaaTCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                                                                                                                                                                                                    <  6:89604/113‑1 (MQ=255)

                                                                                                                                                    aaaTCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa                                                                                                                                                                                                                                                                                    >  5:89604/1‑113 (MQ=255)

                                                                                                                                                        cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                                                                                                                                                                                                          >  5:388128/1‑135 (MQ=255)

                                                                                                                                                        cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa                                                                                                                                                                                                                                                          <  6:388128/135‑1 (MQ=255)

                                                                                                                                                                       tCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCAGCAGCAAATCAAATAATTACGTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCGGTAATGACCTACGCTGTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAAc                                                                                               <  2:413456/275‑1 (MQ=255)

                                                                                                                                                                         cATCAGCGGGTTGAACGTGGCAACCCACGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACa                                                                                              <  4:297764/274‑1 (MQ=255)

                                                                                                                                                                          aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTaat                                                                                                                                                                                                                                                                               >  8:136624/1‑95 (MQ=255)

                                                                                                                                                                          aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTaat                                                                                                                                                                                                                                                                               <  7:136624/96‑2 (MQ=255)

                                                                                                                                                                              ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa                                                                                                                                                                                                                                                         <  5:1022495/114‑1 (MQ=255)

                                                                                                                                                                              ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa                                                                                                                                                                                                                                                         >  6:1022495/1‑114 (MQ=255)

                                                                                                                                                                              ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca                                                                                         <  2:194879/274‑1 (MQ=255)

                                                                                                                                                                                   ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa                                                                                                                                                                                                                                                                               <  7:1034647/87‑1 (MQ=255)

                                                                                                                                                                                   ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTcc                                                                                   >  2:136583/1‑275 (MQ=255)

                                                                                                                                                                                      aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                                                                                                                                                                                                           <  4:552341/152‑1 (MQ=255)

                                                                                                                                                                                      aaCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAgg                                                                                <  4:186751/275‑1 (MQ=255)

                                                                                                                                                                                       aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                                                                                                                                                                                               <  6:959639/99‑1 (MQ=255)

                                                                                                                                                                                       aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa                                                                                                                                                                                                                                                               >  5:959639/1‑99 (MQ=255)

                                                                                                                                                                                       aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAg                                                                                 >  2:393911/1‑273 (MQ=255)

                                                                                                                                                                                          tGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGGAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCcg                                                                            <  2:456303/275‑1 (MQ=255)

                                                                                                                                                                                                accGAGACGGACAGTGAAAGGGATGGCTTTTATAGAAGGATAATGAAGGATGCAACGGCCGGAAAGCGATTAAACGCNATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTTTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGa                                                                      <  3:123305/274‑1 (MQ=255)

                                                                                                                                                                                                     gACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCNGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCACGATAGAAANCATGTCCATTTCCGGATACGGTTTTTTGAANNGACCACATTCCAGGTCCCTGTGGAGaaac                                                                 >  3:68114/1‑273 (MQ=255)

                                                                                                                                                                                                            aTTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCgg                                                                                                                                                                                                                                                                     >  4:533764/1‑72 (MQ=255)

                                                                                                                                                                                                                aCAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagc                                                                                                                                                                                                                                                                 <  2:533765/72‑1 (MQ=255)

                                                                                                                                                                                                                 cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                                                                                                                                                                                                   <  8:183110/85‑1 (MQ=255)

                                                                                                                                                                                                                 cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                                                                                                                                                                                                   >  7:183110/1‑85 (MQ=255)

                                                                                                                                                                                                                    ggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAGATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCCGGATGTTCGGAGTc                                                  >  2:89487/1‑275 (MQ=255)

                                                                                                                                                                                                                      aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGt                                                   >  1:17300/1‑272 (MQ=255)

                                                                                                                                                                                                                            ttaataaagggATAATGAGGGAGGAAACNGNGGGAAGGAGATTAACCGCCATCCGGCAGCAAATCAAAGAATTACTNCGCCGGAATTTCTTTGAGCAGTTCTGTCAGCAGTGTCCAGTAATGACCTACGCTTTGGATGTGAACTTGNTCATCCGGGGGGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATGTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGaa                                          <  3:48311/266‑1 (MQ=255)

                                                                                                                                                                                                                                               ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGttta                        >  4:463242/1‑273 (MQ=255)

                                                                                                                                                                                                                                                gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTc                                                        >  4:108376/1‑241 (MQ=255)

                                                                                                                                                                                                                                                  tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg                                                                                                                                                                                      <  6:580112/113‑1 (MQ=255)

                                                                                                                                                                                                                                                  tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg                                                                                                                                                                                      >  5:580112/1‑113 (MQ=255)

                                                                                                                                                                                                                                                  tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa                                                                                                                                                                                     <  7:969253/114‑1 (MQ=255)

                                                                                                                                                                                                                                                      aCTGCCGGAAGGCGATTAAACGCCATCCGGCATCAAATCAAATAATTCCTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAcc                 <  2:455140/274‑1 (MQ=255)

                                                                                                                                                                                                                                                        tGCCGGAAGGCGATTAAACTCCATCCGGCAGAAAATCAAATAATTGCTTCGCCGGAATTTCATTCAGAAGTTCAGTAAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGGGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGa              <  4:164540/275‑1 (MQ=255)

                                                                                                                                                                                                                                                         gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa                                                                                                                                                                               <  8:470931/113‑1 (MQ=255)

                                                                                                                                                                                                                                                         gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa                                                                                                                                                                               >  7:470931/1‑113 (MQ=255)

                                                                                                                                                                                                                                                          ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                                                                                                                                                                                                              >  4:31649/1‑65 (MQ=255)

                                                                                                                                                                                                                                                           cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg                                                                                                                                                                                                   <  7:542539/91‑1 (MQ=255)

                                                                                                                                                                                                                                                           cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg                                                                                                                                                                                                   >  8:542539/1‑91 (MQ=255)

                                                                                                                                                                                                                                                             gatggCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc                                                                                                                                                                                       >  6:330517/4‑101 (MQ=39)

                                                                                                                                                                                                                                                             gatggCGTTTAATCGCCTTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc                                                                                                                                                                                       <  5:330517/98‑1 (MQ=39)

                                                                                                                                                                                                                                                              aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                                                                                                                                                                                                   >  7:931845/1‑40 (MQ=255)

                                                                                                                                                                                                                                                              aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                                                                                                                                                                                                                   <  8:931845/40‑1 (MQ=255)

                                                                                                                                                                                                                                                              aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                                                                                                                                                                                                                          <  2:31650/65‑1 (MQ=255)

                                                                                                                                                                                                                                                                ggCGATTAAACGCCATCCGGCAGCAAATTAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTAGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAc      <  4:206007/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                                                                                                                                                                                                      <  7:486871/82‑1 (MQ=255)

                                                                                                                                                                                                                                                                 gCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa                                                                                                                                                                                                      >  8:486871/1‑82 (MQ=255)

                                                                                                                                                                                                                                                                    aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCCTCCGGAGAGTGTGGACCGGTGATCGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGCCCCGCGTGAGAAATCTGGATGTTTGGCGTCTTGTTTAACAGGCCCTGAAAGGTTTAAGCTACCAGATGCATccgcggc  >  4:388678/1‑269 (MQ=255)

                                                                                                                                                                                                                                                                    aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGAcca                                                                                                                                                     <  6:92815/128‑2 (MQ=255)

                                                                                                                                                                                                                                                                    aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCg                                                                                                                                                     >  5:92815/1‑128 (MQ=255)

                                                                                                                                                                                                                                                                    aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCCGCGTCTTGTTGAACAGGCGCTGATAGGTTTTACGTACCAGATGCATCACCgg   >  1:66655/1‑274 (MQ=255)

CGACTCGGCGGTACGCGCCACGCATTTGGCATCCGGTATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGA  >  NC_000913/253968‑254502

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