BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,156,590	C→A	6.6%	intergenic (+72/‑33)	yqhD → / → dkgA	aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,156,590	0	C	A	6.6%	673.9 / 27.7	197	intergenic (+72/‑33)	yqhD/dkgA	aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand): ref base C (84/100); new base A (9/4); total (93/104)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-02

Only 100 of 199 total aligned reads displayed.

ATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGA  >  NC_000913/3156321‑3156861

aTTGACGCCGCGATTGCCGCAACCCGCAATGTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACGGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCt                                                                                                                                                                                                                                                                             <  1:121010/274‑1 (MQ=255)

    aCGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCGATTAGGCGTGCCGACCCACCTCTCCGGCTACGGTCTGGACGGCAGCGCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                                                                                                                                                                                                         <  2:256026/274‑1 (MQ=255)

     tgctgggATTGCGGCAACTCGCAAGTTCTTTGGGGAATGAGGCGTGCTGACCCGCCTATCCGCCGCCGGTCTGGACGGCAGTGCCATCCCGGCTGTGCGGAAAGAACTGGAATGGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTGTGATGTCAGCCGCCGTAGATACGAAGCCTGCCGCTCAGCGTTTTACGGCTCAAACTTTCGTTTTCGGGCATTGCGTCCAGACTTAAGTTGACAACACCTCTCCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        <  4:226591/268‑1 (MQ=255)

                                   tGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCACGACTTAAGTTCACAACACCTCACCGGACCCTGCTCCGGTAAGTTCATACAAAGGGGGAAAGTATGGCCaa                                                                                                                                                                                                                                         >  4:591872/1‑275 (MQ=255)

                                         aTTAGGCGGGCCGACCCACCCTGCCGACTACGGGCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACGGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                                                                                                                                                   <  4:532345/275‑1 (MQ=255)

                                            aGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGt                                                                                                                                                                                                                                <  4:263848/275‑1 (MQ=255)

                                                                         gTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGACTGCTCCGGTGAGTTCATATAAaggagg                                                                                                                                                                                                                                                      >  2:32389/1‑224 (MQ=255)

                                                                          tCTGGACGGCAGCGCCATCCCTGCTTTGCTGAAAAAACTGGAAGAGCACGACATGACCACACTGAGCGAAAATCATGACATTACGTTGGATGTCAGACGCCGTATATACGAAGCCGCCCGCTAAGCTGTTTACGCCTCAAACTTTCGTTTTCGGGCACTTCGTCCAGACTTAAGTTCACAGCACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                                                                                                                                                                                                  <  2:463212/275‑1 (MQ=255)

                                                                          gttgggcgcgaGCTCCATCCCGGCTTTGCTGAAAACACTGGAAGAGCACGGCATTACCCAACTGGGCGAAAATCATGACATTACATTGGATGTCAGCCGCCGTATATACGAGGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATACAGGAGGAACGTGTGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                                                                                                                                                                                                  <  4:31026/265‑1 (MQ=255)

                                                                                   cAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCt                                                                                                                                                                                                                       >  4:290871/1‑245 (MQ=255)

                                                                                       tCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                                                                                                                                                                                                   <  2:290871/245‑1 (MQ=255)

                                                                                          aacccGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                                  <  1:67649/273‑1 (MQ=255)

                                                                                                   ttGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                                                                                                                                         >  4:93823/1‑275 (MQ=255)

                                                                                                   ttGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCCCAAACCTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCACCGGTGAGCTCATATAAAGGAGGAACGCATGGCTAATACAACGGTTAGTAAGCTACAGGATGGCAAGGGCATGCACCAGCTGGTACTGGACGGCTGgcaa                                                                                                                                                                         >  4:62913/1‑275 (MQ=255)

                                                                                                                        cACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        >  2:335766/1‑159 (MQ=255)

                                                                                                                            gCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                    <  4:335766/159‑1 (MQ=255)

                                                                                                                                gACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        >  2:688056/1‑151 (MQ=255)

                                                                                                                                 aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                                                                                                                                                                                                         >  8:524565/1‑149 (MQ=255)

                                                                                                                                  cccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        <  5:456027/149‑1 (MQ=255)

                                                                                                                                  cccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTGCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        <  8:72647/149‑1 (MQ=255)

                                                                                                                                      aCTGAGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGTTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGGTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTGATGCCCCATCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaa                                                                                                                                      <  1:117398/275‑1 (MQ=255)

                                                                                                                                         gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                                                                                                                                                                                                 <  7:453582/149‑1 (MQ=255)

                                                                                                                                            cGAAAATCAGGACATTACGTTGGATGTCGGCCGCCGTATATACGAAGCCGCCCGCTAATCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                                 <  2:434908/274‑1 (MQ=255)

                                                                                                                                                   cATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                                                                                                                                                       <  5:483576/149‑1 (MQ=255)

                                                                                                                                                   cATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                                                                                                                                                       <  8:848353/149‑1 (MQ=255)

                                                                                                                                                    aTGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                                                                                                                                                                                                      >  7:146416/1‑149 (MQ=255)

                                                                                                                                                     tGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                                                                                                                                                        <  5:896221/146‑1 (MQ=255)

                                                                                                                                                             cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                                                                                                                                                             <  5:529378/149‑1 (MQ=255)

                                                                                                                                                              gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                                                                                                                                                                                            >  8:66186/1‑149 (MQ=255)

                                                                                                                                                                 ggATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                                                                                                                                                         <  8:633800/149‑1 (MQ=255)

                                                                                                                                                                  gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTc                                                                                                                                          >  2:636820/1‑243 (MQ=255)

                                                                                                                                                                  gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTgg                                                                                                                        >  2:169538/1‑261 (MQ=255)

                                                                                                                                                                       cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATAc                                                                                                     >  2:214023/1‑275 (MQ=255)

                                                                                                                                                                       cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAACGCCATTCAAAAAGCg                                                                                                                                  >  2:88838/1‑246 (MQ=255)

                                                                                                                                                                        aGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                                                                                                                                                                                                 >  4:527686/1‑150 (MQ=255)

                                                                                                                                                                            gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTAACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                                                                                                                                >  4:126634/1‑163 (MQ=255)

                                                                                                                                                                                tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTAACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                                                                                                                                            <  2:126634/163‑1 (MQ=255)

                                                                                                                                                                                   ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCt                                                                                                                                                                                                                       >  6:10332/1‑149 (MQ=255)

                                                                                                                                                                                     aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                                                                                                                                                                                                     >  5:782997/1‑149 (MQ=255)

                                                                                                                                                                                     aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACATAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTAc                                                                                        <  2:80184/274‑1 (MQ=255)

                                                                                                                                                                                       gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                                                                                                                                                        <  6:807812/112‑1 (MQ=255)

                                                                                                                                                                                       gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                                                                                                                                                        >  5:807812/1‑112 (MQ=255)

                                                                                                                                                                                       gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                                                                                                                                                                                                      <  7:95435/98‑1 (MQ=255)

                                                                                                                                                                                        aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                                                                                                                                                                                                        >  4:697859/1‑95 (MQ=255)

                                                                                                                                                                                          gccgccCGCTAAGCTTTTTACGCCTCAACCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGTTGCGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCCTGCCCCAGCTGGGACTGGGCGTCTGGCAAGCCATTAATGAGGAAGTACACACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCCCGATTTATACCACCGCGGCCTACAAGaa                                                                                   >  4:136764/1‑274 (MQ=255)

                                                                                                                                                                                          gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                                                                                                                                                                                           >  4:435982/1‑106 (MQ=255)

                                                                                                                                                                                          gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                                                                                                                                <  6:390307/149‑1 (MQ=255)

                                                                                                                                                                                            cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                                                                                                                                                                                              <  7:390148/149‑1 (MQ=255)

                                                                                                                                                                                              cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                                                                                                                                                                                                           <  2:532753/134‑1 (MQ=255)

                                                                                                                                                                                              cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGCTTGATACCGCCGCGGCCTACAAGAACgaag                                                                              >  2:143161/1‑275 (MQ=255)

                                                                                                                                                                                              cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATAcc                                                                                                    >  4:171590/1‑253 (MQ=255)

                                                                                                                                                                                                  cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                                                                                                                                                   <  4:656836/122‑1 (MQ=255)

                                                                                                                                                                                                   tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                                                                                                                                                                                                       <  7:753020/149‑1 (MQ=255)

                                                                                                                                                                                                    aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                                                                                                                                                                                                      <  7:719947/149‑1 (MQ=255)

                                                                                                                                                                                                    aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                                                                                                                                                                                                      >  6:342277/1‑149 (MQ=255)

                                                                                                                                                                                                      gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                                                                                                                                                                                                    >  8:554617/1‑149 (MQ=255)

                                                                                                                                                                                                       cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                                                                                                                                                                                                   >  7:573968/1‑149 (MQ=255)

                                                                                                                                                                                                             aCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                                                                                                                                                                                             >  5:1063861/1‑149 (MQ=255)

                                                                                                                                                                                                              cGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTa                                                                                                                    <  2:618952/221‑1 (MQ=255)

                                                                                                                                                                                                                ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg                                                                                                                                                                                          <  5:235025/149‑1 (MQ=255)

                                                                                                                                                                                                                 cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                                                                                                                                                         >  6:669379/1‑149 (MQ=255)

                                                                                                                                                                                                                 cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                                                                                                                                                         >  5:13816/1‑149 (MQ=255)

                                                                                                                                                                                                                  tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTa                                                                                                                                                                                                                      >  6:490577/1‑119 (MQ=255)

                                                                                                                                                                                                                  tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                                                                                                                                                        <  7:61849/149‑1 (MQ=255)

                                                                                                                                                                                                                  tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                                                                                                                                                        <  8:811443/149‑1 (MQ=255)

                                                                                                                                                                                                                      aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                                    <  7:865646/149‑1 (MQ=255)

                                                                                                                                                                                                                      aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg                                                                                                                                                                                    <  8:918424/149‑1 (MQ=255)

                                                                                                                                                                                                                        tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc                                                                                                                                                                                  >  5:669797/1‑149 (MQ=255)

                                                                                                                                                                                                                        tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaa                                                                                                                                      <  2:401692/193‑1 (MQ=255)

                                                                                                                                                                                                                             ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                                                                                                                                             <  6:1059465/149‑1 (MQ=255)

                                                                                                                                                                                                                               ttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgc                                                                                                                                                                           <  6:922175/149‑1 (MQ=255)

                                                                                                                                                                                                                                 cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                                                                                                                                         <  6:201205/149‑1 (MQ=255)

                                                                                                                                                                                                                                 cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                                                                                                                                         <  5:408940/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGCCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGa                                                                                                                                                              >  7:352514/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAg                                                                                                                                                             <  8:122181/149‑1 (MQ=255)

                                                                                                                                                                                                                                             ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAg                                                                                                                                                             >  7:983867/1‑149 (MQ=255)

                                                                                                                                                                                                                                             ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAg                                                                                                                                                             <  7:678550/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                                          <  7:935073/149‑1 (MQ=255)

                                                                                                                                                                                                                                                caCTTAAGGGCACAACACCACTCCGGAGCCTGCACCGGTGAATTCATATAAGGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATCGCAATTTCATGCTCTAGCTGGGACTGGGCGTCTGGCAAGCAAGTAACGAGGaa                                                                                                                                                          <  6:611415/148‑1 (MQ=255)

                                                                                                                                                                                                                                                   ttAAGTTCACAACACGTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAGGGCTAATCCAACCGTTATTAAGCTACAGGCTGGCAGTGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                                                                                                                                                       >  6:13616/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                                                                                                                                                       >  5:33089/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aaGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                                                                                                                                                        >  7:256731/1‑66 (MQ=255)

                                                                                                                                                                                                                                                     aaGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                                                                                                                                                        <  8:256731/66‑1 (MQ=255)

                                                                                                                                                                                                                                                          cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTAc                                                                                        <  2:179052/205‑1 (MQ=255)

                                                                                                                                                                                                                                                           acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGc                                                                                                                                               >  6:555533/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaagacCTCACCGGAGCCTGCTGCGGTGAGTTGAAATAAAGGAGCAACGTGTGGGTAATCCAACCGTTATTAATCGACAGAGTGGCAAGTGCATGCCCAGGCTGTGACTGGGCGTATGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTAAAAAAGCGTTAGAAGTGGGTTATCGCTAGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTAAACAGAGAAGAACTGTTCATCACctct                <  4:516680/271‑1 (MQ=255)

                                                                                                                                                                                                                                                             aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa                                                                                                                                             <  6:762137/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTAATTAAGCTACAGGATGGCAATGTCATGCCCCGCCTGGGACTGGGCGTTTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                                       >  6:986146/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                               <  7:227179/125‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                               >  8:227179/1‑125 (MQ=255)

                                                                                                                                                                                                                                                                   tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                                       >  8:1015226/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa                                                                                                                                       >  6:895572/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTc                                                                      >  4:188031/1‑212 (MQ=255)

                                                                                                                                                                                                                                                                     aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa                                                                                                                                     >  7:306277/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa                                                                                                                                     >  6:41089/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aCCGAAGCCTGCACCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACGGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaaa                                                                                                                                     <  8:150436/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc                                                                                                                                   <  8:621138/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                                 <  8:398444/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            ccTGCTCCGGTGAGTCCATATAAAGGAGGAACGTATGGCTCACCCAACCGCTATTAAGCTACAGGATGGAAATGTCTTGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAGGATGAACAAACAACCGCCATTCANAAAGCGTTAGAATTGGGTTATCGTCCGACTCATACTCCCCGGGCCTACAAGAACGAAGAAGGTGCAGGCAAGGCNCTGCAAAACGCCACAGTCATCAAAGAAAAANNCTGCAACACCCCTAAGCAGTGGAaccc  >  3:86826/1‑271 (MQ=255)

                                                                                                                                                                                                                                                                            ccTGCTACGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAg                                                                                                                              >  8:480976/1‑149 (MQ=255)

ATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAAGTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGA  >  NC_000913/3156321‑3156861

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