BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	1,647,682	0	T	G	22.8%	218.1 / 21.7	148	intergenic (+261/+112)	dusC/ECDH1_RS08050	tRNA dihydrouridine synthase DusC/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (90/8); new base G (15/14); total (119/29)
Fisher's exact test for biased strand distribution p-value = 5.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCACCACCAGT  >  DH1/1647532‑1647830

ccctGGCTGGGGGCTGCGGGCGCGTTTTACCGGCCGGTTCCCTCGCGCCCGACCCGGTCGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGGGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCt                                                                                                                                                        <  5:408820/148‑1 (MQ=11)

 caTGGATGGGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTg                                                                                                                                                       <  4:188132/151‑1 (MQ=14)

  aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGCGGCGGGGCTTGG‑G‑CGCCCCCGg‑‑c                                                                                                                                                    >  8:392143/1‑151 (MQ=11)

  aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCTg‑‑c                                                                                                                                                    >  1:418747/1‑151 (MQ=14)

  aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCCGCGGCGCTTGG‑C‑CCCCCCCTg‑‑c                                                                                                                                                    >  6:213671/1‑151 (MQ=9)

    ggagggggggTGAGGGGGGGTTTTACAGGGATGTTACCTCGCGCCCGACCGGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTg‑‑cg                                                                                                                                                   <  7:967180/144‑1 (MQ=11)

        gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTc                                                                                                                                              >  1:272099/1‑151 (MQ=31)

        gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGCGGGTGCACGGGGCCCCGGCGGTTGG‑G‑CCCCCCCAg‑‑cgcgcgc                                                                                                                                              >  2:1065467/1‑149 (MQ=11)

             cTGAGGCCGCCTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGTCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCGATTGG‑C‑CGACCCCAG‑‑CCCCCTCCTTgc                                                                                                                                         >  6:623117/1‑151 (MQ=11)

              tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGGGATTGG‑C‑CGCCCCCGG‑‑CCCGCTCCGTgcg                                                                                                                                        >  4:1039051/1‑151 (MQ=255)

               gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTgcgc                                                                                                                                       >  4:674785/1‑151 (MQ=31)

                   gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGt                                                                                                                                   >  6:119409/1‑151 (MQ=31)

                     gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTg                                                                                                                                  >  2:267554/1‑150 (MQ=31)

                       gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGAtt                                                                                                                               >  3:746830/1‑151 (MQ=18)

                        ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGCGCCCGGGGGTGCACGGGGCGGCGGCGCTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGAtt                                                                                                                               >  2:373621/1‑150 (MQ=255)

                        ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCGTGG‑C‑CGCCCCCGG‑‑CGCGCGCCGTGCGGCAGTGAttt                                                                                                                              >  3:151758/1‑150 (MQ=11)

                         tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATtat                                                                                                                             >  2:577885/1‑151 (MQ=35)

                          ttACACGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCAGAGCCCGGGGGTGCAGGGCGCCCCGGGGAATGG‑G‑GGCCCCCCG‑‑CCCGCACCATGCGCCCGTgtaaaat                                                                                                                            >  8:655553/1‑145 (MQ=11)

                           tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTatat                                                                                                                           >  1:59451/1‑151 (MQ=37)

                           tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGATTatat                                                                                                                           >  6:305481/1‑151 (MQ=37)

                           tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCCG‑‑CCCGCTCCTTTCGCCAGTGATTatat                                                                                                                           >  1:484903/1‑151 (MQ=255)

                             cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATg                                                                                                                          >  5:684157/1‑150 (MQ=255)

                              aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGG‑C‑CGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGATTATATGa                                                                                                                         >  6:355213/1‑150 (MQ=255)

                                   tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAAt                                                                                                                   >  1:743398/1‑151 (MQ=255)

                                   tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGATTTTATGAAGCaa                                                                                                                    >  8:811088/1‑150 (MQ=255)

                                      tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGc                                                                                                                 >  2:483957/1‑150 (MQ=255)

                                      tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGc                                                                                                                 >  5:1079976/1‑150 (MQ=255)

                                      tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGGGTGG‑C‑CGCCCCCGG‑‑CGCCCGCCCTGGGCCAGGGCTTATATGAAGCAATGcc                                                                                                                >  1:888701/1‑150 (MQ=14)

                                        ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCACCTTGCGCCAGTGATTATATGAAGCAATGcac                                                                                                               >  4:286420/1‑150 (MQ=255)

                                          tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑G‑CGCCCCCTG‑‑CGCGCTCCGTGGGCCAGTGATTATATGAAGCAATGCACCTg                                                                                                            >  3:280434/1‑151 (MQ=255)

                                           cgcgCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGATTATATGAAGCAATGCCCATgc                                                                                                           >  3:851762/1‑150 (MQ=255)

                                           cgcgCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCCGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTTCGCCACTGATTATATTAAGCAATGCACATg                                                                                                            >  6:831915/1‑150 (MQ=255)

                                            gcgcCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCGTGCGCCACCGATTATATGAATCAATTcacaata                                                                                                           >  7:1074133/1‑147 (MQ=14)

                                              gcCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGAc                                                                                                          >  2:865071/1‑149 (MQ=255)

                                                  gACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCGG‑‑CGCGCGCCGTGCGCCAGTGGATATATGAACCAAAGCACCTgcacggcgc                                                                                                    >  1:319969/1‑143 (MQ=11)

                                                  gACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAgc                                                                                                    >  4:242266/1‑151 (MQ=255)

                                                      cGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCa                                                                                                 >  2:187762/1‑150 (MQ=255)

                                                         tAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGGTAGG‑C‑CCCCCCCAG‑‑CGCGCTGCTTGGGCCACTGATTTTATGAACCAATGGCCAACACACGAGCGCAACac                                                                                             >  7:97482/1‑150 (MQ=14)

                                                          aGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCGG‑‑CGCGCTCCTTGCGCCAGTGATTATATTAAGCAATGCACATGACAGGAGCGCAACAgg                                                                                            >  1:816782/1‑151 (MQ=255)

                                                           gCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCCCTCCTTGCGCCAGTGGTTATATGAAGCAATGCACATGACAGGCGCGCCACAgg                                                                                            >  8:274340/1‑150 (MQ=255)

                                                             cGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATg                                                                                         >  6:512703/1‑151 (MQ=255)

                                                               tAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCGG‑‑CGCGCGCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGCGCGCAACAGGATg                                                                                         >  3:292023/1‑149 (MQ=255)

                                                                    gATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTTGG‑C‑CGCCCCCGG‑‑CGCGCGCCGTGCGCCACGGTTTATATGAAGCAAAGGCCAAGCCAGGCGCGGAACAGGATGAGAGAcg                                                                                  >  7:622088/1‑150 (MQ=14)

                                                                     aTAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGGGATTGG‑C‑CGCCCCCTG‑‑CGCGCTGCTTGCGCCAGTGATTTTATGAAGCAATGCACCTGACAGGAGCGCAACAGGATGAGAGACTg                                                                                 >  7:428618/1‑151 (MQ=255)

                                                                     aTAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCGCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATgagagtcgg                                                                                 >  7:692896/1‑147 (MQ=255)

                                                                       aaGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCCGGGGCGGGGGGGCGTTG‑G‑CGCCCCCAG‑‑CGCGCGCCATGGGGCCGGGGTTTTATTTAGCAATGGGCATGCCAGGGGCGCCACCGGCTgggagccggg                                                                                >  5:959759/1‑145 (MQ=2)

                                                                       aaGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCCGCGGCGCTTGG‑G‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCCCTGATTATATGAAGGAAAGCACAAGACAGGCGCGCCACCCGATGCGAGCCTGa                                                                                >  3:704427/1‑150 (MQ=17)

                                                                           cGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCAG‑‑CGCGCTCCATGCGCCAGTGATTATATTAAGCAATGCACATGACAGGAGCGCAACACGATGAGAGACTGAAatat                                                                           >  2:875344/1‑151 (MQ=255)

                                                                                gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACAGAattattgg                                                                        >  1:544831/1‑147 (MQ=255)

                                                                                 cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACAGGCCAGGAGCGCCACAGGCTGAGCGACTGAAATATtctggg                                                                     >  3:31223/1‑148 (MQ=18)

                                                                                 cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGCTTGG‑C‑CGCCCCCTG‑‑CCCGCTCCGTGCGCCCGTGATTTTATTAAGCAAAGCACATGCCAGGCGCGCAACAGGCTGACAGACTGAAATATTCTGGt                                                                     >  3:998628/1‑151 (MQ=18)

                                                                                    cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGGGATTGG‑C‑CGCCCCCTG‑‑CGCCCTCCTTGCGCCAGTGATTATATGAAGCAATGCACAAGACACGAGGGCAACAGGTTGAGAGACTGAAATATTTTCGTTAt                                                                  >  5:1055167/1‑151 (MQ=21)

                                                                                         aGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTATAGTTTTTTTCg                                                             >  6:490037/1‑151 (MQ=255)

                                                                                         aGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCCGCGGCGGTTTG‑C‑CGCCCCCGG‑‑CGCGCTGCGTGCGGCACTGATTATATGAAGCAAAGCCCATGACAGGGGGCCCAAAGGGGGAGAGACTGAAAATTTTGGGTTATTTTCg                                                             >  8:576280/1‑151 (MQ=12)

                                                                                          gcggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGATTATATGAAGCAATGCACCTGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCGCGGTATTTTCgg                                                            >  4:235190/1‑151 (MQ=255)

                                                                                           cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACACGATGAGAGACTGAAATATTCTCGTTATTTTCggc                                                           >  3:395527/1‑150 (MQ=255)

                                                                                              cgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTTTATTAAGCAATGCACATGACAGGAGCGCAACAGGATGAGCGACTGTAATATTCTCGGTTTTTTCGGACCt                                                        >  5:422/1‑151 (MQ=21)

                                                                                               gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCCTTATTTTAGGACCtg                                                       >  1:615408/1‑151 (MQ=255)

                                                                                                tttttttgcccctgcgCCGGGGGGGGCTGGGGGGGGGGGCGTTGGC‑C‑CCCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCtg                                                       <  2:52034/135‑1 (MQ=25)

                                                                                                aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGG‑C‑CGCCCCCGG‑‑CGCGCGGCGTGCGCCACTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGCTGGGAGCCTGAAATATTCTCGTTATTTTCGGGCCtgg                                                      >  7:1004795/1‑150 (MQ=21)

                                                                                                aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACACGACCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCtg                                                       >  3:565063/1‑150 (MQ=255)

                                                                                                         gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGCGCGCAACAGGATGAGAGACTGAAATATTATCGTTATTTTAGGACCTGATGACaata                                              >  8:122465/1‑150 (MQ=255)

                                                                                                              agCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTATCGTTATTTTAGGACCTGATGACATtaattaat                                        >  4:985235/1‑151 (MQ=255)

                                                                                                               ggggggggtgggctggggggggggggTTTGG‑CGCGCCCCCGC‑‑CGCGCTCCTTGCGCGAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATtaattaat                                        <  2:593538/137‑1 (MQ=255)

                                                                                                                gggggggggGCGGGGGGGGGTGGGGTTTGG‑G‑CCCCCGCGGC‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTaattaatt                                       <  2:803881/148‑1 (MQ=255)

                                                                                                                 ccGGGGGTGCACGGGCCGGCGGCCATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACAAGACAGGAGCGCAACAGGATGAGAGACAGAAATAATCTCGGTATTTTCGGACCGGATGACATTAATTAAtttt                                     >  4:917935/1‑151 (MQ=37)

                                                                                                                  gggggggGCTGGGGG‑GGGGGTGTTTGGCC‑CCCCCCCCC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTa                                    <  6:149427/150‑1 (MQ=255)

                                                                                                                  gggggggGCGGGGGGGGGGGGCTAGGGC‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTa                                    <  6:154636/150‑1 (MQ=38)

                                                                                                                  ggggggTGGGGGGGGGGGTGG‑GTTTGG‑G‑CCCCCCCGGC‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTa                                    <  6:722653/150‑1 (MQ=255)

                                                                                                                  cGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCCATGCACATGACAGGAGCGCCACAGGATGAGAGACTGAAATATTCTCGTTATTTTAGGACCTGATGCCATTAATTAAtttt                                     >  5:1000301/1‑150 (MQ=255)

                                                                                                                  cGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGCGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGCCATTAATTAAtttt                                     >  3:856637/1‑150 (MQ=255)

                                                                                                                     ggggggggggggCGGCGGCTAGCGC‑C‑CGCCCCGGG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTaaaa                                 <  2:982417/144‑1 (MQ=255)

                                                                                                                     ggggcggggggggggCGG‑GGTTGG‑C‑CCCCCCCCGCCCGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTaaa                                  <  5:318160/144‑1 (MQ=255)

                                                                                                                     gggTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTaaa                                  >  7:227126/1‑150 (MQ=255)

                                                                                                                          cAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATatc                             >  2:406966/1‑150 (MQ=255)

                                                                                                                                gggggggttgtgcgcc‑cgcccgcCC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTa                       <  8:526690/136‑1 (MQ=255)

                                                                                                                                  ggggggtatggg‑C‑CGCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAAt                    <  2:652463/141‑1 (MQ=255)

                                                                                                                                  ggcggcgATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAAt                    >  2:669324/1‑151 (MQ=255)

                                                                                                                                   gcggcgATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAAt                    >  5:278952/1‑150 (MQ=255)

                                                                                                                                    gggcttgggg‑C‑CCCCCCGCG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATAc                  <  4:771126/143‑1 (MQ=255)

                                                                                                                                    gcggctctGG‑C‑CCCCCCCCGC‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATa                   <  6:693651/144‑1 (MQ=255)

                                                                                                                                       cgATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGcc               >  4:812851/1‑151 (MQ=255)

                                                                                                                                        ggtgcggccccccccccc‑‑cGCGCTCCTTGCGCGAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGc                <  2:525190/140‑1 (MQ=255)

                                                                                                                                        gctccggccccccccccc‑‑cGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGc                <  1:267778/140‑1 (MQ=255)

                                                                                                                                        gATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGcc               >  5:1622/1‑150 (MQ=255)

                                                                                                                                         aTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGcc               >  5:711457/1‑149 (MQ=255)

                                                                                                                                         aTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGcc               >  5:146915/1‑149 (MQ=255)

                                                                                                                                         aTTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCat             >  8:731059/1‑150 (MQ=255)

                                                                                                                                          gtggccccccccccGC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCa              <  1:902086/142‑1 (MQ=255)

                                                                                                                                           tGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTc           >  3:1046372/1‑151 (MQ=255)

                                                                                                                                            gg‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACAGGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTcc          >  5:203468/1‑150 (MQ=255)

                                                                                                                                             ggccccccccccc‑‑cGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTc           <  5:490671/145‑1 (MQ=255)

                                                                                                                                             ggcccccccccGC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTc           <  6:542181/145‑1 (MQ=255)

                                                                                                                                             gcccccccccccc‑‑cGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTc           <  3:1027928/145‑1 (MQ=255)

                                                                                                                                              gccgccccccTC‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTca          <  3:44482/146‑1 (MQ=255)

                                                                                                                                                  gccCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTcaccac      >  7:86761/1‑151 (MQ=255)

                                                                                                                                                   ccccccc‑‑cGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCaccacc     <  1:535199/151‑1 (MQ=255)

                                                                                                                                                    ccccTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCaccacc     >  3:75736/1‑150 (MQ=255)

                                                                                                                                                      ccTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCACCACCAGt  >  2:977282/1‑151 (MQ=255)

                                                                                                                                                      ccTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCACCACCAGt  >  4:983153/1‑151 (MQ=255)

ACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑C‑CGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGATTATATGAAGCAATGCACATGACAGGAGCGCAACAGGATGAGAGACTGAAATATTCTCGTTATTTTCGGACCTGATGACATTAATTAATTTTAAAATATCATCTTAAATACGCCAGTCACCACCAGT  >  DH1/1647532‑1647830

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