BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	3,970,143	G→T	E6* (GAG→TAG)	ECDH1_RS19620 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,970,143	0	G	T	98.2%	1683.5 / 13.7	449	E6* (GAG→TAG)	ECDH1_RS19620	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (7/1); new base T (245/196); total (252/197)
Fisher's exact test for biased strand distribution p-value = 8.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 449 total aligned reads displayed.

GCCAACACTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATAT  >  DH1/3969996‑3970290

  cAACACTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAggtg                                                                                                                                                >  7:719163/1‑151 (MQ=255)

  cAACACTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAggtg                                                                                                                                                >  8:400189/1‑151 (MQ=255)

    acacTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAggtggt                                                                                                                                              <  1:17289/151‑1 (MQ=255)

      acTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGttt                                                                                                                                            >  2:620451/1‑151 (MQ=255)

      acTCCACAGTTTTCCCCCAAAAGATGAAACATTTGGCGACATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTATCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGCTCGttt                                                                                                                                            >  7:655393/1‑151 (MQ=255)

         gCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTg                                                                                                                                         >  1:1423360/1‑151 (MQ=255)

             aGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGGATTAGGTGGttttttg                                                                                                                                        >  4:290107/1‑146 (MQ=255)

              gTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGgagag                                                                                                                                    <  3:1355168/151‑1 (MQ=255)

                 ttCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGTCAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGagaa                                                                                                                                  <  1:1306840/150‑1 (MQ=255)

                  tCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGagaag                                                                                                                                 <  3:93488/150‑1 (MQ=255)

                    ccccAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAATCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAAt                                                                                                                              <  7:1472867/151‑1 (MQ=255)

                     cccAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGAATTAGGTGGTTTTTGGAGAGAAGAAtt                                                                                                                             >  8:1076915/1‑150 (MQ=255)

                      ccAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGa                                                                                                                            <  2:154391/151‑1 (MQ=255)

                        aaaaGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAgg                                                                                                                          >  3:951070/1‑151 (MQ=255)

                             aTGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGAt                                                                                                                     <  5:229630/151‑1 (MQ=255)

                             aTGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGAt                                                                                                                     <  1:73840/151‑1 (MQ=255)

                                aaaCATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCg                                                                                                                  <  2:470355/151‑1 (MQ=255)

                                aaaCATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCg                                                                                                                  >  1:239555/1‑151 (MQ=255)

                                  aCATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTc                                                                                                                >  2:931897/1‑151 (MQ=255)

                                  aCATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTc                                                                                                                <  2:1021968/151‑1 (MQ=255)

                                  aCATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTc                                                                                                                >  7:22653/1‑151 (MQ=255)

                                   cATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCg                                                                                                               >  4:513714/1‑151 (MQ=255)

                                        ggggCCATTTTGACTCATAGAGGAAAGCATCGCGGACATAGTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCca                                                                                                          <  6:1217291/151‑1 (MQ=255)

                                        ggggCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCca                                                                                                          <  6:433578/151‑1 (MQ=255)

                                                 ttGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGt                                                                                                 >  1:882154/1‑151 (MQ=255)

                                                   gACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGtt                                                                                                >  7:527439/1‑150 (MQ=255)

                                                     cTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGc                                                                                             >  7:1147343/1‑151 (MQ=255)

                                                      tCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCt                                                                                            >  5:713369/1‑151 (MQ=255)

                                                       cATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCt                                                                                            <  1:1174340/150‑1 (MQ=255)

                                                                 aGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGc                                                                                 >  7:1284423/1‑151 (MQ=255)

                                                                   cATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCgg                                                                               <  3:1089501/151‑1 (MQ=255)

                                                                    aTCGCGCACAAACATTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTATAATTATAGATGTATATATTCCATTACCTCCTTTTTCGAGACAAGAATAAGGAAGATGCGACGAAACACATAAACGATAACTTAACAtataaacc                                                                               >  6:931769/1‑146 (MQ=255)

                                                                      cgcgGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAgg                                                                            >  5:608834/1‑151 (MQ=255)

                                                                      cgcgGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAgg                                                                            <  1:60643/151‑1 (MQ=255)

                                                                      cgcgGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAgg                                                                            >  6:1432763/1‑151 (MQ=255)

                                                                      cgcgGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAgg                                                                            >  4:987511/1‑151 (MQ=255)

                                                                         ggACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGa                                                                         <  6:96368/151‑1 (MQ=255)

                                                                           aCAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGAtg                                                                       <  1:420623/151‑1 (MQ=255)

                                                                             aaaCTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGAtgtg                                                                     <  6:1150604/151‑1 (MQ=255)

                                                                             aaaCTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGAtgtg                                                                     <  3:1199480/151‑1 (MQ=255)

                                                                              aaCTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGa                                                                    <  4:677298/151‑1 (MQ=255)

                                                                               aCTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTaaa                                                                   <  8:650692/151‑1 (MQ=255)

                                                                               aCTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGa                                                                    >  2:520761/1‑150 (MQ=255)

                                                                                  ttttCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGaaa                                                                  <  4:540134/149‑1 (MQ=255)

                                                                                  ttttCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGAATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAaatt                                                                >  4:18546/1‑151 (MQ=255)

                                                                                    ttCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAaatta                                                               <  1:329826/150‑1 (MQ=255)

                                                                                    ttCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAaatta                                                               >  3:417763/1‑150 (MQ=255)

                                                                                     tCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAaattaat                                                             >  7:1210947/1‑151 (MQ=255)

                                                                                       attttAATTTTTTGCTTAATAATCTATTTTTATGTTTATTTTTAGCTGTTTAAATTGGATTATGTGGTTTTTTGAGATCACCATGCCGAACATAAGTCGCTCCCCACTAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAAttt                                                           <  8:1311925/149‑1 (MQ=255)

                                                                                        gTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAAttt                                                           >  8:306868/1‑150 (MQ=255)

                                                                                        gTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAAttt                                                           >  2:82176/1‑150 (MQ=255)

                                                                                        gTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAAttt                                                           >  7:93706/1‑150 (MQ=255)

                                                                                        gTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTa                                                          <  7:399871/151‑1 (MQ=255)

                                                                                         tttATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTAc                                                         >  5:1196198/1‑151 (MQ=255)

                                                                                           tATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACa                                                        <  3:647701/150‑1 (MQ=255)

                                                                                             tttGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAAt                                                      >  7:1077500/1‑150 (MQ=255)

                                                                                             tttGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAAt                                                      >  8:191108/1‑150 (MQ=255)

                                                                                               tgttgGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAAt                                                      >  8:939621/1‑148 (MQ=255)

                                                                                                gttgGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTCCATTAGCTGGTTTTTGGAGAGAAGAATGAGGAAGATCCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGTTGATGTGAAATTAATTTACAATAGa                                                   >  2:987212/1‑150 (MQ=255)

                                                                                                 ttgGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGaa                                                  >  2:308150/1‑150 (MQ=255)

                                                                                                   ggCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAaata                                               <  3:787735/151‑1 (MQ=255)

                                                                                                     cTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAaataat                                             >  6:644580/1‑151 (MQ=255)

                                                                                                      ttAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAAtt                                            >  4:226621/1‑151 (MQ=255)

                                                                                                       tAATCTTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAaataat                                             <  8:1157978/149‑1 (MQ=255)

                                                                                                       tAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAaataat                                             <  4:109339/149‑1 (MQ=255)

                                                                                                        aaTATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAAttt                                           <  6:241528/150‑1 (MQ=255)

                                                                                                         atatTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTAc                                         >  5:930520/1‑151 (MQ=255)

                                                                                                           atTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACat                                       >  2:207677/1‑151 (MQ=255)

                                                                                                            ttCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACat                                       >  7:456959/1‑150 (MQ=255)

                                                                                                             tCTATTGTTATGTTTATTTATAGATGTTTATATTGGATTAGGTGGTTTTTGTAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGGTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACata                                      <  5:1410395/150‑1 (MQ=255)

                                                                                                             tCTATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGACGAAGATGAGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACatat                                     >  2:1260605/1‑151 (MQ=255)

                                                                                                               tATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATc                                    >  5:1296373/1‑150 (MQ=255)

                                                                                                               tATTGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATc                                    <  2:73281/150‑1 (MQ=255)

                                                                                                                 ttGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGGTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCaa                                  >  6:1274111/1‑150 (MQ=255)

                                                                                                                  tGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCaaa                                 >  2:369931/1‑150 (MQ=255)

                                                                                                                  tGTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCaaa                                 <  5:259011/150‑1 (MQ=255)

                                                                                                                   gTTATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAAc                                <  3:539728/150‑1 (MQ=255)

                                                                                                                     tATCTTTATTTATAGATGTTTATATTTCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAg                              <  3:381330/150‑1 (MQ=255)

                                                                                                                     tATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGt                             >  5:477788/1‑151 (MQ=255)

                                                                                                                     tATCTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGt                             >  7:137262/1‑151 (MQ=255)

                                                                                                                        cTTTATTTATAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGtt                            >  8:76093/1‑149 (MQ=255)

                                                                                                                          ttatttatAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGtt                            <  1:1054396/147‑1 (MQ=255)

                                                                                                                           tatttatAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGAt                        <  5:368829/150‑1 (MQ=255)

                                                                                                                           tatttatAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATg                       >  3:808766/1‑151 (MQ=255)

                                                                                                                           tatttatAGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATg                       >  2:768874/1‑151 (MQ=255)

                                                                                                                               tataGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCtt                    <  4:886896/150‑1 (MQ=255)

                                                                                                                                ataGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCtttt                  >  4:763449/1‑151 (MQ=255)

                                                                                                                                ataGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCtttt                  >  3:19043/1‑151 (MQ=255)

                                                                                                                                ataGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttt                   <  8:1440409/150‑1 (MQ=255)

                                                                                                                                 taGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttttt                 >  8:1134029/1‑151 (MQ=255)

                                                                                                                                 taGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttttt                 >  1:526162/1‑151 (MQ=255)

                                                                                                                                 taGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttttt                 >  1:480350/1‑151 (MQ=255)

                                                                                                                                  aGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttttt                 <  7:320127/150‑1 (MQ=255)

                                                                                                                                  aGATGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCttttt                 <  3:250128/150‑1 (MQ=255)

                                                                                                                                    aTGTTTATATTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTc              <  4:1404424/151‑1 (MQ=255)

                                                                                                                                       tttATATTGCATTAGGTTGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGt            <  5:768681/150‑1 (MQ=255)

                                                                                                                                         tataTTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGttt          <  5:556909/150‑1 (MQ=255)

                                                                                                                                           tatTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGtttttt       >  6:861525/1‑151 (MQ=255)

                                                                                                                                            atTGCATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTa      <  4:119539/151‑1 (MQ=255)

                                                                                                                                                cATTAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAatat  <  3:1262054/151‑1 (MQ=255)

GCCAACACTGCACAGTTTTCCCCCAAAAGATGAAACATTTGGGGCCATTTTGACTCATAGAGGAAAGCATCGCGGACAAACTTTTTCAGTTTATTTGTTGGCTTAATATTCTATTGTTATCTTTATTTATAGATGTTTATATTGCATGAGGTGGTTTTTGGAGAGAAGAATGAGGAAGATGCGTCGAGCCACAGAAACGTTAGCTTTACATATAGCGGAGGTGATGTGAAATTAATTTACAATAGAAATAATTTACATATCAAACAGTTAGATGCTTTTTGTCGTTTTTTAATAT  >  DH1/3969996‑3970290

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