BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,522,824	A→G	10.7%	S424S (TCT→TCC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,522,824	0	A	G	10.7%	345.1 / 3.8	136	S424S (TCT→TCC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (59/55); new base G (16/3); total (78/58)
Fisher's exact test for biased strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TCCAACGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCGGGCCATTAGGCAGTAC  >  NC_000913/4522587‑4523060

cccAACGCCAAAAATGACCAGAAGTCCGCCCGCCTTGACTTTGTCTGGTTTGAAATCGTAAAGCCCACTTTTCCCGGAAAGATAAACCTTCGCCACACCTTGATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGga                                                                                                                                                                                                                                              <  2:477241/237‑1 (MQ=255)

          ataatGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGCAGGGG‑ATCGTTGAGGTTAATCCGCGCGC‑CTTT‑GCTTG‑CCGCGGAGACATGTTc                                                                                                                                                                                                                                     >  1:941408/1‑237 (MQ=255)

               gccccaatatcccccccccTTGATTTTTTCATGTTTTATTTCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGt                                                                                                                                                                                                                              <  2:854495/226‑1 (MQ=255)

                aCCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGGTTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGGAGGGG‑GTAGTTGAGGGTAATCCGCGCGC‑TTTT‑GCTTG‑CTTGGGAGAAATGAACACTTCGta                                                                                                                                                                                                                             >  2:799276/1‑238 (MQ=255)

                         tcccccggcTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGAt                                                                                                                                                                                                                    <  1:982170/236‑1 (MQ=255)

                         tCGCCCGCCTTGGTCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTTCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCCCGCGC‑TTTT‑GCTTG‑CTGCGGGGCCATGATCACCTCCTCAGccgttat                                                                                                                                                                                                                    >  1:2427/1‑233 (MQ=255)

                         cccccctgcTTTATTTTTTCATGTTTTATTTCGTAAAACGCATTTTTCGTGGACAGATAAACCTTCGCCACGCCTTTATGGTAATAGATTCCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTTGCGCAATATTGCCCATCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGAt                                                                                                                                                                                                                    <  2:2427/231‑1 (MQ=255)

                               gccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑GTAGTTGAGGTTAATCCGCGCGC‑TTTT‑CCTTG‑CTGCGGGGACCTTATACCTTCGTCCACGTTGGTTTGtt                                                                                                                                                                                                               >  1:765720/1‑238 (MQ=255)

                               gccTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGtt                                                                                                                                                                                                               >  1:793077/1‑238 (MQ=255)

                                      tCTTGTCATGTTTGATAACGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATGCCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCg                                                                                                                                                                                                       <  1:396318/239‑1 (MQ=255)

                                         tGTCTTGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCGTCGCCACACCTTTATGGTAATAGATCCCTGGCCCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGc                                                                                                                                                                                                    <  2:93564/239‑1 (MQ=255)

                                          gTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑GTAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGACATGATCAATTCGTCAGCGTTGATTTGTTCCTGATCGAGCa                                                                                                                                                                                                   >  1:690133/1‑239 (MQ=255)

                                           tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATTACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑AAAGGTGAGGGTAATCCGCGCGC‑TTTT‑GCTTT‑CTCCGGAGACATGATCACTTCGTCAGCCGTGTTTTTTTCCCGTTCGGGCAg                                                                                                                                                                                                  >  1:383698/1‑239 (MQ=255)

                                           tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTTGCGCAATATTGCCCACCGGGAAGGCG‑ATAGTTGAGGTTAATCCCCCCGC‑TTTT‑GCTTG‑CTGCGGGGACATGTTCACTTTGTCAGCGTTGATTTGTTCCTGCTCGCGCAg                                                                                                                                                                                                  >  1:192310/1‑239 (MQ=255)

                                           tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCCCTTCGTCAGCGTTGAATTGTTCCtttttgagccg                                                                                                                                                                                                  >  1:903767/1‑229 (MQ=255)

                                           tCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGGGCAg                                                                                                                                                                                                  >  1:635799/1‑239 (MQ=255)

                                             attgttttattctgaaaactccgttttCCCGGAACGATAAACCATAGCCCCACCTTTAGGGTAATAGAATCCTTGCTCATCAATCAGCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGa                                                                                                                                                                                                <  2:78188/216‑1 (MQ=255)

                                               gTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCCATATTGCCCACCGGGAAGGGG‑GTAGTTGTGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGCCATGATCACTTCGTCAGCGTTGCTTTGTTCCTGATCGAGCAGGGGt                                                                                                                                                                                              >  1:132952/1‑239 (MQ=255)

                                                     tatCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGCACCTTCTGGGGCCATATTGCCCCCCGGGGCAGTG‑ATAGCTGGGGGTAATACGCGCGC‑TTTT‑TCTTT‑CTGCGGGGACCTGTTCCCTTCGTCAGCCTTCATTTTTTCCTTCTTGCGCAGGAGGTTCTTg                                                                                                                                                                                        >  2:168964/1‑239 (MQ=255)

                                                          tAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTTAGGTTAATCCGCGCGG‑TTTT‑GCTTG‑CTGCGGCGACCTGATCACTTCGTCACCGTTGATTTTTTCCTGATCGAGCAGGAGGTGCTTGAAccc                                                                                                                                                                                   >  1:688653/1‑238 (MQ=255)

                                                             aTCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCCCAATATTGCCCACCGGGAAGGTG‑GTAGTTGAGGTTAATCCGCGGGC‑TTTT‑GCTTT‑CTGCGGGGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCt                                                                                                                                                                                 >  1:771906/1‑238 (MQ=255)

                                                               cccccTTTTCCCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGAc                                                                                                                                                                              <  1:712422/235‑1 (MQ=255)

                                                                gtgaataagagacagaGAAAAACCTTCGCCACACCTTTATGGTAATAGATTCCTTGCTCAACAAGCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGAc                                                                                                                                                                              <  1:648232/223‑1 (MQ=255)

                                                                         tcGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTTGCGCAATATTGCCCACCGGGGAGGGG‑GTAGTTGAGGGTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGACATTATCCCTTCGTCACCGTTGATTTGTTCCTGATCGGGCCGGGGTTGCTTGAACCCCTGACGCCGTTCggg                                                                                                                                                                    >  2:1006858/1‑238 (MQ=255)

                                                                           gggacgATAACCCTTCCCCACACCTTGTGGGTAATAGATACCTTGCTCCTCAATCAACGGGGGGT‑CAATGGCGGTTGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGt                                                                                                                                                                  <  2:481526/234‑1 (MQ=255)

                                                                           ggACAGATAACCCTACGCCACATCTTTATGGTAATAGATTCCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGt                                                                                                                                                                  <  1:320935/239‑1 (MQ=255)

                                                                                         tCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAACGTTGATTTTTTCCTGATCGAGCAGGAGTTGCTTGAACCCCTGACGCCGTTCGGAGGGGTCCCACCAGTCg                                                                                                                                                    >  1:867544/1‑239 (MQ=255)

                                                                                         tCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGGTGATTTTTTCCTGATCGGGCAGGAGTTGCTTGAAACGCTGACGCCGTTCGGGGGGCTCccaccaagcg                                                                                                                                                    >  1:668002/1‑235 (MQ=255)

                                                                                           gCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGACATGATCACTTCGTCAGCGGTGGTTTGTTCCTGATTGAGCAGGGGTTGCTTGAACCGCCGTCGCCGTTTGGAGGGGTCCCACCCGTCGAg                                                                                                                                                  >  1:970319/1‑239 (MQ=255)

                                                                                           gCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGTTCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAg                                                                                                                                                  >  1:225272/1‑239 (MQ=255)

                                                                                              cccccTTTATGGTAATTGATACCTTGCTCATCAATCAGCGAGGGGG‑CAATGGCGGTTGATTTGATCCCCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGtt                                                                                                                                               <  2:132952/236‑1 (MQ=255)

                                                                                                acCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATAACCGAACCTTCTGGCGCAATATTGCCCACCGGGGAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGACCTGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCCCCAGTCGAGGtttt                                                                                                                                             >  1:460733/1‑239 (MQ=255)

                                                                                                           aaTAGATCCCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtg                                                                                                                                  <  1:58555/239‑1 (MQ=255)

                                                                                                            atagaatccctgcatatcaatcaTCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGT‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtg                                                                                                                                  <  1:828785/223‑1 (MQ=255)

                                                                                                            aTAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTTAGGTTAATCCCCGCGC‑TTTT‑GCTTT‑CTGCGGGGACATTATCACTTCGTCAGCGTTGATTTTTTCCCTGTCGGGCGGAGGTTTCTTGAACCCCTGGCGCCGTTTGGGGTGTCCCCCCCCGTACAGGGttttttttcccgtgc                                                                                                                                 >  2:155187/1‑231 (MQ=255)

                                                                                                            aTAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGga                                                                                                                                                                                                                                              <  2:373588/130‑1 (MQ=255)

                                                                                                            aTAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGga                                                                                                                                                                                                                                              >  1:373588/1‑130 (MQ=255)

                                                                                                              aagtgtataagagacagagcagcGAGGGGT‑CAATGGCGGTTGATTTGAACCCCGAACCTTCTGGCCCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgct                                                                                                                                <  1:814103/216‑1 (MQ=255)

                                                                                                                    ctgtcTCATCAACCAGCGGGGGGT‑GAAGGGGGGTGGATTTGATCCCCGAACCTTCTGGCGCAATTTGGCCCACGGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCCCTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTa                                                                                                                         <  2:832868/235‑1 (MQ=255)

                                                                                                                    cTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGGAGGTG‑ATAGTTGAGGGTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGGGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTTGGGGTGGTCCCACCAGCCGAGGGTTTTTTTCAGCGTGCTCCCGGTa                                                                                                                         >  2:96576/1‑239 (MQ=255)

                                                                                                                    cTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTa                                                                                                                         >  1:295405/1‑239 (MQ=255)

                                                                                                                      tGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTTTCGCCATATTTCCCACCGGGGAGGTG‑ATTGTTGCGGTTAATCCGCGCGC‑TTTT‑TCTTT‑CTCCGGGGACGTTATCACTTCGTCAGCGGTGTTTTGTTCCTGGCCGAGCAGGGTGTGCTTTAACCCCTGACGCCGTTTTGAGGGCCCCCCCCCTGTGGGGGTTTCTTTCatggggcctccgggta                                                                                                                        >  2:828785/1‑223 (MQ=255)

                                                                                                                      ctgttcttcattctcGGGGGGTCCAATGGGGGTGTATTTGATCACCGACCCTTCTGGCCCAATATTGCCCACCGGTAAGGTG‑ATAGTTGACGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGGGGAGACATGATCACTTCGGCAGCGGTGATTTGTTCCTGATCGGGCAGGAGTTGCTTTAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGGTCCCGGTaa                                                                                                                        <  2:611736/226‑1 (MQ=255)

                                                                                                                       gcccacccagcatcGGGGGGT‑CAATGGGGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGT‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAAcc                                                                                                                      <  2:48915/229‑1 (MQ=255)

                                                                                                                        cTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACCCCGTTCGGCGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgcgtcccggaacc                                                                                                                      >  1:674325/1‑227 (MQ=255)

                                                                                                                              aatcaTCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGGGTGCTCCCCCCCGTGCCGGGTTTCTTTCAGCGTGCTGCCGGTAACCGTCCTa                                                                                                                >  1:48915/1‑238 (MQ=255)

                                                                                                                                       ggggT‑CATTGGCGGTGGTTTTGATCTCCGAACCTTCTGGCGCAATTTTGCCCACCGGGAAGGTG‑ATGGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGAGCCCTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCa                                                                                                       <  2:941408/238‑1 (MQ=255)

                                                                                                                                             catgggCGGTGGATTTGGTCACCGAACCTTCTGGCGCAATTTTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaa                                                                                                  <  2:793077/234‑1 (MQ=255)

                                                                                                                                             acagggtggTGGATTTGGGCCCCGAACCTTCTGTGGAAATATTGCCCACCGGG‑AGGTGTATAGGTGAGGTTAACCCGCGGCT‑TTTT‑GCTTG‑CTGCGGAGACATGACCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGCACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaa                                                                                                  <  2:397538/231‑1 (MQ=255)

                                                                                                                                              taTGGGGTTGGGTTTGTTCCCCGTACTTTTTGGCGCAATTTTGCCCCCCGGGAAGGTG‑TTAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACTTGATCCCTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaaca                                                                                                <  2:383698/238‑1 (MQ=255)

                                                                                                                                              aaTGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGGGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGGGCGCCCGGTAACCGTCATAACCGCTTCATGccaca                                                                                                >  1:1077274/1‑239 (MQ=255)

                                                                                                                                                     gTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGGG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGCGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGTGCAGGAGTTGCTTGAACCGCTGCCGCCGTTCGGAGTGCTCCCACCAACCGCGGGTTTCTTTCAGCGTGCTGCCGGTAACCGCCCTAACGCCTTCATGCAACAAACCGAg                                                                                         >  2:291503/1‑239 (MQ=255)

                                                                                                                                                          tttGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTgc                                                                                    >  1:808066/1‑239 (MQ=255)

                                                                                                                                                           ttgtcccccGAACCTTTTGGCGCAATATTGCCCACCGGGAAGGTT‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTgcg                                                                                   <  2:342435/232‑1 (MQ=255)

                                                                                                                                                           ttGAGCACAGAACCTTCTGGTCAAATATTGCCCTCCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGGCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTgcg                                                                                   <  1:1079333/239‑1 (MQ=255)

                                                                                                                                                              aTCACCGAACCTTCCGGCGCAAAATTGCCCACCGGGAAGGTT‑ATGGTTGGGGTTAATCCGAGCGC‑TTTT‑ACTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATTGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTg                                                                                     <  1:548296/234‑1 (MQ=255)

                                                                                                                                                               tccccGACCCTTCTGGCGCAATTTTGCCCCCCGGGAAGGGG‑ATAGTTGAGGTTAGACCGGGCGT‑TTTT‑GCTTG‑CTGGGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                                                               <  2:771906/236‑1 (MQ=255)

                                                                                                                                                               tCACCGAACCTTCTGGCGCAATATTGCCCCGCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGG‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                                                               <  1:1006858/239‑1 (MQ=255)

                                                                                                                                                               cccccGAACCTTCTGGCGAAATTTGCCCCCCCGGGAAGGTG‑AAAGGTGAGGTTAACCCGC‑CGCTTTTT‑GCCTG‑CTGGGGGAAAAGGATCCCTTCGTCAGCGGTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTGTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                                                               <  2:751373/236‑1 (MQ=255)

                                                                                                                                                                   cGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTaa                                                                                                                        <  2:600396/194‑1 (MQ=255)

                                                                                                                                                                   cGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTaa                                                                                                                        >  1:600396/1‑194 (MQ=255)

                                                                                                                                                                    gAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                                                          >  2:849534/1‑239 (MQ=255)

                                                                                                                                                                        cTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑CCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGATTGCTCCCACCAGCCGCGGTTTTCTTTCGCCGGGCTGCCGGGCACCGTCATAACGTCTTCTTGCCACAACCCGAGGCCGCGCaaaagaaactgg                                                                      >  2:548296/1‑228 (MQ=255)

                                                                                                                                                                             ggCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACtt                                                                  >  1:517152/1‑238 (MQ=255)

                                                                                                                                                                             ggCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTc                                                                 <  1:735207/239‑1 (MQ=255)

                                                                                                                                                                                   atatTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGCCGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGAAACCGTCATAACGTCTTCAGTCAACAATCCGGGGGTCGGCAGATGCAACAAGACTACCGCca                                                            >  2:240054/1‑238 (MQ=255)

                                                                                                                                                                                    tatTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGAGGCAACATGACTTCCGGcac                                                           >  1:756695/1‑238 (MQ=255)

                                                                                                                                                                                         cccccggggAGGGTG‑ATAGTTGGGTTTATGCCGCGCGC‑TTTG‑TCTTG‑TTGGGGGGACAGGGTCTTTGCGTCAGGGTTGATTGTGTGCGGACGGAGCAGTAGTGGCTTGACCCGCTGCCGCCGTTCGGAGTGCTCCCACGAGGCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAAATGACTTCCGGCACaccccc                                                     <  2:556924/233‑1 (MQ=255)

                                                                                                                                                                                         ccccccGGAAAGGTG‑TTTGTTGGGGTTAA‑CCCCGCGC‑TTTT‑TGTTGTCTGCGGGGACATGTTCCCTTCGTCGGCGTTGATTTGGTCCTGATCGAGCAGGGGTTTCTTGACCCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACaccccc                                                     <  2:688653/239‑1 (MQ=255)

                                                                                                                                                                                          cccgcggggaggggtttAGTTGAGTTTAATCCG‑GGGC‑TTTTGGTTTG‑CGGCGGAGAGTTGACCCTTACGCGGGGGTTGTTTTGTCCCTGGTCGAGCAGGAGTTGCTGGAACCGTGGGCGCCGTGCGGGTTGCTCCCACCAGTGGAGGTTTTCTTACAGCGGGCGGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAAAGGCAACATGACTTCCGGCACaccacc                                                     <  2:988327/223‑1 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGa                                                                                                                                                                                                          <  2:122414/85‑1 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCGGATGCAACATGACTTCCGGCACACCACCtccaa                                                >  2:666518/1‑235 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGcct                                                >  1:140348/1‑238 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCa                                                >  1:797738/1‑239 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGCCCTAACGTCTTCATGCAAAAATCCGAGGCTGGCCAGATGCACCCTGACTTCCGGCCAACCACCTGCCa                                                >  2:374463/1‑239 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCAACCAGTCGAGGTTTTCTTTCACCGTGCTGCCGGTCACCGTCATAACGTCTTCATGCAACAATCCCAAGCTGCGCAGATGCAAACTGACCTCTGGCACACACCCTGccc                                                >  2:1152/1‑238 (MQ=255)

                                                                                                                                                                                              cGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGAACACTACGTCAGCGTTGATTTGTTCCTGATCGAGCCGGAGTTGCTTGAACCGCGGCGGCCTTTCGGAGTGCTCCAACCAGTCGAGGTTTTCTTTCATCGTGCTGCCGGGAACGGTCATAACGCTTTCATGCAACAATCCAAGGTGGCCGAGATGCCAACTAACTTCCGCCACACAAACGGCCa                                                >  2:845794/1‑239 (MQ=255)

                                                                                                                                                                                                          gggTTGGGTTTAATCCGGGCGT‑TTTT‑GCTTG‑GTTGGGAGACAGTTTCATTTCGGCAGGTTTGATGTGTCCCGAATCGAGCAGGAGTTGCTTGAACCGCTGCAGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCACCAATCCGAGGCTGCGAAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAtt                                      <  2:648837/236‑1 (MQ=255)

                                                                                                                                                                                                           aGTTGAGTTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTAAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGTCTGCGCAGATGCACCATGATTTCCGGCACACCACTTGCCATAAAGGCATTg                                     >  1:506711/1‑238 (MQ=255)

                                                                                                                                                                                                            tttGGGGGTAACCCGGGCTC‑TTTT‑GTTTG‑CGGGGGAGAAATGATCCCTTCGTCGGGGTGGATTTGTTCCTGATCGGGAAGGGTTTGCTTGAACCGCGGCCGCCGTTCGGATGTCTCCCACCAGTCGGGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACAGGACTTCCGGCACACCACCTGCCATAAAGGCATTGAc                                   <  2:756695/238‑1 (MQ=255)

                                                                                                                                                                                                              tGAGGTTAATCCGCGCGC‑TTTG‑GCTGG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGATTTGCTGGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCg                                 <  1:848046/239‑1 (MQ=255)

                                                                                                                                                                                                                aGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCCTAACGTCTTCATGCAACAATCCGAGGCTGCGCCGATGCAACATGACTTCCGGCACACCCCCTGCCATAAAGGCATTGACAGtt                               >  1:586434/1‑239 (MQ=255)

                                                                                                                                                                                                                aGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGtt                               >  1:1010211/1‑239 (MQ=255)

                                                                                                                                                                                                                aGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCAAAAAGGAATTgcacgtg                               >  2:259369/1‑233 (MQ=255)

                                                                                                                                                                                                                 ggTTAATCCGGGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGtt                               <  1:1012190/238‑1 (MQ=255)

                                                                                                                                                                                                                 ggTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGCCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGACCACCACCTCCCCAAAAGGCATTAACCGtgg                              >  2:940070/1‑237 (MQ=255)

                                                                                                                                                                                                                       tCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtg                                                                                                                                  >  2:810999/1‑133 (MQ=255)

                                                                                                                                                                                                                       tCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtg                                                                                                                                  <  1:810999/133‑1 (MQ=255)

                                                                                                                                                                                                                           cgcgc‑TTTT‑GCTTG‑CTGCGGAGGCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                                                               <  2:244304/180‑1 (MQ=255)

                                                                                                                                                                                                                           cgcgc‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCGGATGCACCATGACTTCCGGCCACCCCCCGGCCAAAAAGGCATTGACCGTTGgagaaaaaac                    >  2:600016/1‑231 (MQ=255)

                                                                                                                                                                                                                           cgcgc‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGa                                                                               >  1:244304/1‑180 (MQ=255)

                                                                                                                                                                                                                              gc‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGATTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  1:155187/238‑1 (MQ=255)

                                                                                                                                                                                                                              gc‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGACCCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  1:96576/238‑1 (MQ=255)

                                                                                                                                                                                                                              gc‑TTTT‑GCTGG‑CTGGGGAGCCTTGATCCCTTCTTCACCGTTTATGTGTTCCTGATCGAGCAGGAGTTGCGTGAACCGCTGACCCCGTTCGGAGTGCTCCCACCATTCGAGGTTTTTTTTCAGCGTGCGGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  1:817127/238‑1 (MQ=255)

                                                                                                                                                                                                                               gttttt‑GCTTG‑GTGGGGAGACATGATCCCTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  1:993704/236‑1 (MQ=255)

                                                                                                                                                                                                                               ggtttc‑GCTTG‑CTGGGGAGACATGAGCACTTGGTCAGGTTTGATTTGTTCCTGATCGAGCGGGATTTGCTTGAACCGCGGCCCCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTCCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAACCCGAGGCGGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  2:687671/236‑1 (MQ=255)

                                                                                                                                                                                                                                gtttt‑GCTTG‑CTGGGGAAACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCggg                <  1:259369/238‑1 (MQ=255)

                                                                                                                                                                                                                                ggctttgttctggcggggAGATCTGATCCCTTCGGCAGCTTTGATTTGTGCCTGATCCAGCAGGGGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGGCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCg                  <  1:862953/223‑1 (MQ=255)

                                                                                                                                                                                                                                  gttttttttggggcGGAGAAATGATCCTTTGGTCGGGTTTTGTTTGTTCCGTATCGGGCAGGATTTGGTTGATCCGGTGACGCCGTGCGGAGTCCTCCCACCAGTCGAGGTTTTCTTTCAGCGGGCGGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCggg                <  2:690133/227‑1 (MQ=255)

                                                                                                                                                                                                                                             tgcGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCGGGCCATTAGGCAg     >  1:653965/1‑239 (MQ=255)

TCCAACGCCAATAATGACCAGAATATCGCCCGCCTTGATCTTGTCATGTTTGATATCGTAAATCGCACTTTTCTCGGACAGATAAACCTTCGCCACACCTTTATGGTAATAGATACCTTGCTCATCAATCATCGAGGGGT‑CAATGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTG‑ATAGTTGAGGTTAATCCGCGCGC‑TTTT‑GCTTG‑CTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCATTGACCGTTGGATGATAAACCGGGCCATTAGGCAGTAC  >  NC_000913/4522587‑4523060

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