BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	1,065,614	0	A	G	24.3%	251.1 / 93.7	168	intergenic (+49/‑14)	ECDH1_RS05245/ECDH1_RS05250	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (66/43); new base G (31/4); total (100/68)
Fisher's exact test for biased strand distribution p-value = 1.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 192 total aligned reads displayed.

GCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGC  >  DH1/1065467‑1065763

tatgtaatgtaaGATTCCCTATGCATCCGTAGCTGATCTGTATAGAGTACTCGGCTACGAAGCGATCGGTCTGAGGTTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAaaata                                                                                                                                                      <  8:64100/147‑1 (MQ=21)

gCAGAAAAGTAATTTGCCCTCTGCATCCGGAGCTCACCGGGAGAGAGGAATCGGCTACTAACGGTAGGCGCCGAGGTTCGAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAgaataa                                                                                                                                                     <  3:953812/151‑1 (MQ=9)

     aaagtaaGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAATACTAAATATAGCTATAAAAACGGTATA‑‑AAT‑‑ATCAAAAATAATATAA‑Attt                                                                                                                                                >  4:479414/1‑150 (MQ=14)

       agtaagtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTCATAGATATAGCGTAATGAACGGAAAC‑‑AAT‑‑ATGAACAGTAATTTAA‑ATCTccc                                                                                                                                             >  8:1157929/1‑151 (MQ=14)

         tttcctttcttttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAAGCGATCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTcc                                                                                                                                              <  6:192850/136‑1 (MQ=21)

               gccatttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATCGTACTTGATATGGGGTAAGTAGCGGTGTC‑‑AAT‑‑TTAAGCAGTAAAAAAA‑AAATCCCCATGcg                                                                                                                                       >  7:1474432/7‑148 (MQ=14)

                   gATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGGTACGAACCGAGCGGTAGGAGGTTCGAATCCTACTGGATGCACCATATCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTa                                                                                                                                 <  5:557156/151‑1 (MQ=255)

                              aGCTCATCTTGATTGAGTAATCGGCTACTAAGGGTTGGGTCTGAGGTTGGAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAgaataagt                                                                                                                                                   <  4:88402/123‑3 (MQ=11)

                                        gtttttctaatCGGCGACGATCCCACCGCGCGGCGCTTCTAATCCTCCGGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACt                                                                                                             <  7:1110951/140‑1 (MQ=14)

                                                aCTCGGCTACGAACCGTTGGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAGAATAAGTTTTCCCGATGCATCCGTa                                                                                                                                 <  5:402910/123‑1 (MQ=11)

                                                 cTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGGAGTAAAATAA‑ATATCCCGATGCATCCGTAGCTGAGCTGGATAGAGTACTCGGCTACGa                                                                                                    >  7:73415/1‑150 (MQ=255)

                                                  tCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAGAATAA‑GTTTtccccttgca                                                                                                                                      >  1:826596/1‑108 (MQ=17)

                                                   ggcacgtaacccGAGCGGCCGG‑GGGTGGAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAGAATAAGTTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAAcc                                                                                                 <  3:265661/141‑1 (MQ=11)

                                                          gatcgctccGGTCGGCGCTTCTAATCCTCCCTTATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGt                                                                                          <  3:977382/143‑1 (MQ=14)

                                                              gctccgctCGGAGAGTCTAATCCTCCCGGATTCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGa                                                                                      <  7:854225/144‑1 (MQ=14)

                                                               gtatgttctgAGGTTCTTATCCTCCGGTATTCAAGATCTTTTTTTTGATATGTCTTTTTTAGGGTTATT‑‑AAT‑‑ATCAGCAGTAAAATAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAg                                                                                     <  4:1318842/142‑1 (MQ=2)

                                                                  cGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGGAGCGGTATG‑‑AAG‑‑AGCAGGAGGAAAATTA‑CTTTCCCTATGCAGCCGTAGCACAGCTGGATAGAGTACTCGGCTGGGCCCGGCGCGGTGGGAGGtt                                                                                  >  7:371251/1‑151 (MQ=17)

                                                                   gctctgcgcgtctaaTCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTc                                                                                 <  7:181930/138‑1 (MQ=255)

                                                                       ggcgagtcttaTCCGCCGGGAGTCACCATCTCTTACTTGATATGGGTTTAGTAGCGCTATA‑‑AAT‑‑ATCAGCAGTAAACTAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAt                                                                             <  1:1027021/141‑1 (MQ=14)

                                                                       ggAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑ATC‑‑TGAAGCAGTAAAGTAT‑ATTTCCCGATCCAGCCGTTCCGGATCTGCTTAGAGATCACGGCTACGAACACAGCGGGCGGAGGttggccg                                                                             >  7:325540/1‑146 (MQ=9)

                                                                         ggcgtctaaTCCTCCCGGATGCACCATCTCTTACTTGATATGTCTTTACTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcc                                                                           <  1:256040/144‑1 (MQ=255)

                                                                           ggtctaaTCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGGAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcctc                                                                         <  1:481673/146‑1 (MQ=255)

                                                                              gtaaTCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCg                                                                      <  1:355662/149‑1 (MQ=38)

                                                                              ataaTCCTCCCTGACTCAATATCTTTTAATTTATATGTCTTTATTAGCTTTATA‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCg                                                                      <  2:675707/149‑1 (MQ=11)

                                                                               gAATCCTCCCGGATGCACCATATAGTACTTGATATTCATGTATTAGCGTTAAC‑‑AAT‑‑ATCAACAGTAAAAAAA‑ATTTCCCGATACATACGAAGAGCATCTCGATAGATTAATCGGCTACGAAACTAGATGTCGGAGGTTCGAATCCTCCCgg                                                                     >  4:987119/1‑151 (MQ=9)

                                                                               gAATCCACCCGGATGCACCATCTCTTACTTGAAAGCGCTATACTAGCGGTATA‑‑AAT‑‑ATCAGCAGTAAATTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGATCGGTCGGAGGTTCGAATCCTCCCgg                                                                     <  7:1007488/151‑1 (MQ=14)

                                                                                aatcctcggcgtcccACTATCTCTTTTTTTATATGTTTTTATTAGCTCTATA‑‑AAT‑‑ATCAGCAGTAAACTAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTccc                                                                       <  2:1155592/134‑1 (MQ=11)

                                                                                aaTCCTCGCTGATGCACCATCTCTTACTTGATATGTCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTccc                                                                       <  1:1247169/148‑1 (MQ=255)

                                                                                   catccCAGATCCACGATCTCTTACTTGATATGGCTATACCAGCGGTATA‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATg                                                                  <  3:371277/148‑1 (MQ=14)

                                                                                        cGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAGCCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCat                                                            >  7:116548/1‑151 (MQ=255)

                                                                                          gttccACGATCTATTACTTGATATCGCTTTAGTAGGGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCatat                                                          <  7:1306046/147‑1 (MQ=255)

                                                                                          gtctgaCCATCTTTTAGTTTATAGGTCTATAGTAGCTCTATA‑‑AAT‑‑ATCAGCAGTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCatat                                                          <  7:594736/146‑1 (MQ=11)

                                                                                               aCCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATtatc                                                      >  2:584460/1‑150 (MQ=255)

                                                                                                ccATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCg                                                    >  3:21283/1‑151 (MQ=39)

                                                                                                 cATCTCTTACTTGATCTGGCTATACGAGCTCTATA‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGa                                                   <  7:1250538/151‑1 (MQ=14)

                                                                                                  atatttttcTTTATATTGTTTTAGTAGAGGTATC‑‑TAT‑‑TTCTTTATTAAAGTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGaa                                                  <  8:427791/143‑1 (MQ=11)

                                                                                                   tctcTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGaa                                                  >  8:29018/1‑150 (MQ=255)

                                                                                                   tctcTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGaa                                                  >  2:952803/1‑150 (MQ=255)

                                                                                                      ctTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATg                                                                                                                                         <  3:495783/60‑1 (MQ=255)

                                                                                                         agttgatcgctcTTTACTAGGGTTAAC‑‑AAT‑‑ATCAGCAGTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCa                                           <  3:1383305/140‑1 (MQ=14)

                                                                                                           ttGATACGGCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTACTAAATATTc                                            >  8:595526/1‑150 (MQ=14)

                                                                                                             gATACGGCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATCCGAAATATTCatt                                         >  6:2542/1‑149 (MQ=14)

                                                                                                             gATACGGCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACGAAATATTCAgt                                         >  8:639253/1‑150 (MQ=14)

                                                                                                               aatgggtattctagCGCTAAA‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATg                                     <  5:705160/140‑1 (MQ=14)

                                                                                                                  gtctttattaggtttata‑‑AAT‑‑ATCTGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGaa                                   <  8:320/136‑1 (MQ=11)

                                                                                                                  ggCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCTACATATTCAGCAAtt                                     >  5:656339/1‑149 (MQ=14)

                                                                                                                  ggCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGATATAttgctctatt                                     >  1:938258/1‑142 (MQ=14)

                                                                                                                  ggCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTATATATTCataaatt                                     >  3:125374/1‑144 (MQ=14)

                                                                                                                  ggCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAAACTTCAGCAATGa                                    >  5:1174962/1‑149 (MQ=25)

                                                                                                                     tttAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTCCGAAATATTCCGAAATGAAgg                                 >  8:1231412/1‑151 (MQ=11)

                                                                                                                       tagtagCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGa                                                                                                    >  2:904922/1‑80 (MQ=255)

                                                                                                                        agtagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAGATATTCCTCAATTAGGGTAt                              >  3:680514/1‑151 (MQ=11)

                                                                                                                        agtagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACTAAATATTGAGCAATGAACGTAt                              >  5:28755/1‑151 (MQ=11)

                                                                                                                        agtagCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCg                            >  5:1348244/1‑151 (MQ=255)

                                                                                                                          ttgCGGTATC‑‑AAAAAATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGc                                         <  6:16225/136‑1 (MQ=11)

                                                                                                                          tagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTGGAGCagctcggatgagc                            >  5:16248/1‑139 (MQ=11)

                                                                                                                          tagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACGAAATCTTCATAAATGAACGTAtac                            >  8:254074/1‑149 (MQ=11)

                                                                                                                           aggtttata‑‑AAT‑‑ATCTGCTGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGaa                          <  2:1099706/145‑1 (MQ=11)

                                                                                                                            gCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCGa                           >  7:454279/1‑148 (MQ=255)

                                                                                                                             cGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtc                                                                            >  1:848295/1‑98 (MQ=255)

                                                                                                                             cGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtc                                                                            <  2:850018/98‑1 (MQ=255)

                                                                                                                             cGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAag                         >  5:639580/1‑149 (MQ=255)

                                                                                                                             cGGTATC‑‑AAT‑‑ATCAGAAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGt                                                                                          <  5:1318075/84‑1 (MQ=255)

                                                                                                                              ggTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAgg                      >  1:359632/1‑151 (MQ=255)

                                                                                                                              ggTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGATACGAACCGAGCGGTCGGAGGTTCGAATCCTCCAGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAgg                      >  6:963084/1‑151 (MQ=255)

                                                                                                                              ggTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGATACGAACCCAGCGGTCGGAGGTTCGAATCCTCCAGGATGCACCATATTCTCCGAAATCTTCAGGAATGAAGGTATCGAAGAgg                      >  2:1406643/1‑151 (MQ=255)

                                                                                                                              ggTATC‑‑AAT‑‑ATCAGCAGT‑AAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATGGAAGAGGt                     >  8:787962/1‑151 (MQ=255)

                                                                                                                              gaggtatcaaaaaatCTGCAGTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCt                                              <  3:1021632/118‑1 (MQ=1)

                                                                                                                               ttata‑‑AAT‑‑ATCTGCAGTAAATTAA‑TTTTCCCTATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAgg                      <  6:1047667/149‑1 (MQ=14)

                                                                                                                               ttata‑‑AAT‑‑ATCTGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAgg                      <  2:1470469/149‑1 (MQ=255)

                                                                                                                                tata‑‑AAT‑‑ATCAGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt                     <  8:1377061/150‑1 (MQ=255)

                                                                                                                                 aTC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt                     >  4:966721/1‑149 (MQ=255)

                                                                                                                                  tatcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATATTCAGCAATGAACGTATCGa                           >  2:1011320/10‑146 (MQ=11)

                                                                                                                                   atcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGACATATTCAGCTATTAAGGTAATGAAGAGGt                     >  3:860734/9‑151 (MQ=11)

                                                                                                                                   atcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGAGGTCTTAGGTTCGAATCCTCAAGGATTCACCATATTATACGAAATATTCAGAAATGAACGTATAGaatatc                      >  6:1168596/9‑146 (MQ=11)

                                                                                                                                   ataaaT‑‑ATCTGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGTCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAg                   <  4:357507/148‑1 (MQ=255)

                                                                                                                                    cgatattataagctttAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAACCTTCAGCAATGAAGGTATCGAAGAGATa                    <  8:343152/136‑1 (MQ=1)

                                                                                                                                      aaT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGACCCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATATTCAGCAATGAAGGAATCGAAGAGGTAGCgg                >  4:1052218/1‑151 (MQ=255)

                                                                                                                                      aaT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAAACTTCAGCAATGAAGGTATCGAagag                       >  6:323878/1‑144 (MQ=25)

                                                                                                                                        t‑‑atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGa               >  2:230551/1‑150 (MQ=255)

                                                                                                                                            tCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAg                                                                                     <  5:167347/78‑1 (MQ=255)

                                                                                                                                             cagcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAACCTTCAGCAATg                                     >  1:788498/1‑125 (MQ=18)

                                                                                                                                             cagcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGATCTCGGCTGCGAACCGCGCCCGCGGAGGTTCGAATACTCCCGGATGCACCATAGTGTCCGAAATCTACAGCAATGAAGGTTGCGAAGGGGGAGCggcgtgc           >  7:69185/1‑146 (MQ=14)

                                                                                                                                               gcagTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTaa          <  1:1008413/150‑1 (MQ=255)

                                                                                                                                                cagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCGAAGAGATAGCgg                >  6:1182213/1‑143 (MQ=255)

                                                                                                                                                 agTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCgg                                                                                       <  2:618972/71‑1 (MQ=38)

                                                                                                                                                 agTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCgg                                                                                       >  1:617609/1‑71 (MQ=38)

                                                                                                                                                 agTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAAccc       >  8:824174/1‑150 (MQ=255)

                                                                                                                                                  gtaaagtaa‑aTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGc      <  1:195145/151‑1 (MQ=255)

                                                                                                                                                  gTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGCAGGTTCGAATCCTCCCGGATGCACCATATTGTCCGAAATATTGCGGAATGAAGGTATCGAATAGGTAGCGGaagatagga       >  5:237830/1‑143 (MQ=14)

                                                                                                                                                  gTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAAACTTCAGCAATGAAGGTATCGAAGAGATAGCGGTATTAACCGc      >  6:284602/1‑151 (MQ=2)

                                                                                                                                                  gTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGc      >  1:152020/1‑151 (MQ=255)

                                                                                                                                                   tAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGc      >  3:1361073/1‑150 (MQ=255)

                                                                                                                                                   tAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGc      >  8:456386/1‑150 (MQ=255)

                                                                                                                                                   tAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAACCTTCAGCAATGAAGGTATCGAAGAGATAGCGGTATTAACCg       >  7:1450549/1‑149 (MQ=9)

                                                                                                                                                    aaagTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACTAGATATTGATGAATGAGCGTTTCGAAGAGCTAGAGTAATTaaaagcac    >  5:1249098/1‑145 (MQ=1)

                                                                                                                                                     aaataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCgg                                                                                           <  7:639429/63‑1 (MQ=39)

                                                                                                                                                     aaataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAAACTTCAGCAATGAAGGt                                >  6:1472801/1‑122 (MQ=18)

                                                                                                                                                     aaataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATACTCAGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGAGGAATTAACAGCAgg   >  6:867111/1‑151 (MQ=255)

                                                                                                                                                      aataa‑atTTCCCGATGCCTCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGCGCGGTCGGAGGTTCGAATACTCCCTGATGCACAATATTCTCCGAAATCTGCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTACCCGCAgcc  >  1:456516/1‑149 (MQ=14)

GCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGC  >  DH1/1065467‑1065763

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: `‑`