BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	1,065,618	0	A	G	31.9%	414.3 / 239.4	314	intergenic (+53/‑10)	ECDH1_RS05245/ECDH1_RS05250	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (121/56); new base G (67/16); total (195/119)
Fisher's exact test for biased strand distribution p-value = 3.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 337 total aligned reads displayed.

AGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACTAG  >  DH1/1065469‑1065768

    aagtaagtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATTCACCATCTCTTACTTGATATTGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTcc                                                                                                                                                   <  3:1820182/151‑1 (MQ=255)

       taagtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTGGTACTAGATATGGAGTTAGTAGCGGTTGC‑‑AAT‑‑ATCAGGAGTAAATTAA‑ATTTcgcg                                                                                                                                                 >  8:57150/1‑150 (MQ=17)

         agtTTCCCTTTGCATACGTAGCTCATCTGGATAGAGTACTCGGCTACTA‑TCGGTTGGG‑TCTCAGGTTGTAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAGAATAAGTTTTCCCGAt                                                                                                                                               <  6:728392/151‑1 (MQ=11)

                atttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGt                                                                                                                                       <  3:479157/148‑1 (MQ=255)

                   tGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTTGTACGTGATATGGCTTTATGAGCGGTATG‑‑AAT‑‑ATCTGCAGTAAATTAA‑ATATCCCGCTgggtcggatgc                                                                                                                                    >  7:490298/1‑141 (MQ=17)

                   tGCAGCGGTAGGTCAGCTGGCTTGACTAATCGGCTAGGATGGCGAG‑GC‑TCGGCGCTTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAg                                                                                                                                     <  7:187822/150‑1 (MQ=14)

                       tCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATGGTACTAGATATGGCGTTAGGAGCGGTTAC‑‑AAT‑‑ATCTTCAGTAATTTAA‑AGTACGCGAttggtaagaacaag                                                                                                                                  >  3:2088528/1‑136 (MQ=9)

                         cGTAGCTCAGCTGGATTGAGTACTCGGCTACGA‑ACGGCGCGG‑TCGGAGGGTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCt                                                                                                                              <  7:1290866/151‑1 (MQ=255)

                            aGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACTAAATATGGCTTTAGGATCGGTAGC‑‑AAT‑‑ATCTCGAGTAATATAA‑AGTACCCGCTACATACGTAGATCAGcggtt                                                                                                                           >  7:407721/1‑147 (MQ=9)

                              gacAGCTGGATATACTACTCGGCGACG‑‑AGCGCGCGGCTCGGCGCGTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATa                                                                                                                         <  1:2216816/149‑1 (MQ=14)

                                     ggATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGTTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGTTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATagag                                                                                                                      <  7:2053295/147‑1 (MQ=255)

                                               cTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCAAGAAATAAATAAAATCGAGTTAGTAAAGGTATA‑‑ATT‑‑ATCAGCAGTAAAATAA‑ATTTCACGAAGCATCAGTAAAAAAGATGGATAGAGTACTCGGCTACGac                                                                                                        >  8:1781686/1‑150 (MQ=12)

                                                  ggCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑Gttttc                                                                                                                                                   >  1:697421/1‑103 (MQ=255)

                                                       cGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTGGTACTTGATATGCATGTAGGAGGGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCCTCCGAAGATCAGCTGTTATCAGTACTCGGCTACGAACCCAGccg                                                                                                >  1:592509/1‑149 (MQ=18)

                                                         a‑aCCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACGTGATATGCCTGTAGTAGCGGTAGG‑‑AAT‑‑ATGAGCAGGAATAAAA‑AGTACCCGATACATAGGAACCCCATCGGTATAGACAACTCCGATACGAACCCAGcactc                                                                                              >  5:2129863/1‑147 (MQ=9)

                                                            ccGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGAGGTAAC‑‑AAT‑‑ATCAGAAGTAAAAAAA‑AGATGCCGATGCATACGTAGCTCAGCTGTAAAGAGTAATCGGATACGAAACGAgcgggcgg                                                                                            >  2:383516/1‑151 (MQ=17)

                                                              gtccGG‑TCGGAGGGTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAg                                                                                          <  7:1275462/148‑1 (MQ=255)

                                                                gCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACTTGATATGGATGTAGTAGGGGTTGC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATATCCCGCTGCGTAAGTAGAACAGCAGGTTAGAGAACTCGGCTACGAGCAGAGCGTTCGGAGGt                                                                                        >  2:1956139/1‑151 (MQ=19)

                                                                   gctcggcgcgtctaaTCCTCGCGGATGCACGATCTCTTACTTTATATGTCTTTAGTAGCTGTATA‑‑AAT‑‑ATCAGCAGTAAATTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTc                                                                                      <  1:267756/138‑1 (MQ=11)

                                                                   gctccgcgcgtctaaTCCTCGCTGTTGCACCATCTCTTACTTGATATGGCTTTACTAGCTGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTc                                                                                      <  8:1442051/138‑1 (MQ=14)

                                                                              cGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTccc                                                                            <  1:1856518/150‑1 (MQ=255)

                                                                               taaTCCTCGCTGATTCACAATCTCTTTCTTTATCTGTTTTTTTGAGCTGTATA‑‑AAT‑‑ATCTGCAGTAAATTAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCg                                                                           <  4:164651/149‑1 (MQ=11)

                                                                                aaTCCTCCCTGATTCACGATCTCTTACTTGATATGGCTTTATTAGCTTTATA‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCgg                                                                          <  8:1298870/150‑1 (MQ=14)

                                                                                aaTCCTCCCGGATGCACCATCTCTTACTTGATATTGCTATAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATACATCCGTAGCTCAGCAGGATAGAGTACTATGCTAAGAACCGAGAGGTGGGAGGTTCGATTCCTCCCgg                                                                          >  4:1510282/1‑150 (MQ=14)

                                                                                       ccGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAaaataa‑at                                                                                                                                                       >  8:2091767/1‑66 (MQ=255)

                                                                                             gCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATtct                                                            >  8:1947509/1‑151 (MQ=39)

                                                                                                 cATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCg                                                         >  5:1459013/1‑150 (MQ=40)

                                                                                                           ttGATACGGCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACTAAATATTcc                                                >  3:2179808/1‑150 (MQ=14)

                                                                                                                acggCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCTACATATTCCTCAAt                                           >  1:929524/3‑151 (MQ=14)

                                                                                                                  gtctttaagcgctatata‑‑AAT‑‑ATCTGCAGTAAAGTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGaa                                        <  7:1466358/136‑1 (MQ=11)

                                                                                                                  ggCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGa                                         >  7:273647/1‑149 (MQ=255)

                                                                                                                   gCTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATACTCCCGGATGCACCATATTGTCCTAAATATTCAGCAATGa                                         >  4:846465/1‑150 (MQ=11)

                                                                                                                    cTTTAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTACATATTCAGATATGaac                                       >  5:984236/1‑150 (MQ=14)

                                                                                                                     tttAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGACATATTCAGCTATGACgg                                      >  5:1773232/1‑151 (MQ=11)

                                                                                                                     tttAGTAGCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATATTCATAAATTAAgg                                      >  4:937797/1‑151 (MQ=11)

                                                                                                                     tttAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGt                                     >  7:1666307/1‑150 (MQ=255)

                                                                                                                     tttAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGt                                     >  8:525933/1‑150 (MQ=255)

                                                                                                                      ttAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTAt                                   >  5:248601/1‑151 (MQ=255)

                                                                                                                      attatagggatata‑‑AAA‑‑ATCAGCAGTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTa                                    <  7:986231/140‑1 (MQ=14)

                                                                                                                       tagtagCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATc                                  >  4:345093/1‑151 (MQ=255)

                                                                                                                        agtagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATACTCCGAAATATTCAGATATGAACGTa                                    >  2:1135523/1‑150 (MQ=11)

                                                                                                                         gtagCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCg                                 >  8:415951/1‑150 (MQ=255)

                                                                                                                          tagCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTAAATCTTCataaa                                            >  4:1176705/1‑136 (MQ=11)

                                                                                                                           cgctctata‑‑AAT‑‑ATCACCAGTAAAATAA‑ATTTCCTGATGCATCCGTAGCTCATCTGGATATACTAATCGGCTACTATGCGTTGGCGCTTTGGGTGGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGa                                <  8:1663821/144‑1 (MQ=9)

                                                                                                                           agCGGTATCAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTACGAAATATTCAGCAATTATGGTAt                                   >  5:317863/1‑148 (MQ=11)

                                                                                                                           agCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCGa                                >  7:1327746/1‑149 (MQ=255)

                                                                                                                            gCGGTATC‑‑AAAAAATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGAGATCTTCAGGAATGAAGGTATGGa                                >  1:2007623/1‑150 (MQ=11)

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

                                                                                                                             cGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAACCTTCAGCAATGAAGGTATCGAag                              >  2:839542/1‑149 (MQ=38)

                                                                                                                              ggaaTC‑‑AAA‑‑ATCCGCCGTAAAATAC‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcctc                                                                              <  1:346987/97‑1 (MQ=11)

                                                                                                                              gcggtgtctaaaaatCAGCAGTAAAGTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGaa                               <  3:390551/138‑1 (MQ=14)

                                                                                                                               cggtatcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGa                                         <  2:187648/128‑1 (MQ=11)

                                                                                                                               ataac‑‑AAT‑‑ATCAGCAGTAAAATAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt                          <  7:2177861/147‑1 (MQ=255)

                                                                                                                                tATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCCTCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGACCCGCGCGGTCGGAGGTACGAATCCTCCCGGCTGCACCATATTGTCCGAAATCTTCAGCTATGAAGGGTTCGAAGAGGtg                         >  6:1798738/1‑150 (MQ=14)

                                                                                                                                  tatcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGACATATTCAGATATGATGGTTTCGAagatg                           >  5:55258/10‑149 (MQ=11)

                                                                                                                                  tatcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATACTCCGAAATATTCAGCAATTAAGGTATACAagatg                           >  6:1368140/10‑149 (MQ=11)

                                                                                                                                  gtcttcatatcttgcagTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAgg                           <  2:984755/138‑1 (MQ=14)

                                                                                                                                   gcttcatatcttgcagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt                          <  1:1297706/139‑1 (MQ=14)

                                                                                                                                   atcaaaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTTCTACTAAATATTCAGGTATGATCGTATACAATAGGt                          >  3:163061/9‑151 (MQ=14)

                                                                                                                                    tcaaaaaatCTGCAGTAAAATAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt                          <  6:55180/143‑1 (MQ=255)

                                                                                                                                     cccaaacgcagcagTAAATTTA‑TTTTCCCGATGCATCCGTAGCTCAGCTGCATACTTCCCTAGTGTACTTAAGTAGTGGGCTGATGCTATAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAg                        <  1:1511112/143‑1 (MQ=9)

                                                                                                                                     ataaaaatCTGCAGTAAAATAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATACTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTa                         <  3:846828/144‑1 (MQ=255)

                                                                                                                                     aaaaaaatCTGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAg                        <  6:209305/145‑1 (MQ=255)

                                                                                                                                     aaaaaaatCTGCAGTAAAGTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTa                         <  3:1172598/144‑1 (MQ=255)

                                                                                                                                     aaaT‑‑ATCTGCAGTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGc                       <  4:1069546/149‑1 (MQ=255)

                                                                                                                                      aaaaaatCTGCAGTAAATTAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGc                       <  6:875796/146‑1 (MQ=255)

                                                                                                                                      aaT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGa                                                                                    >  7:1015032/1‑88 (MQ=255)

                                                                                                                                      aaT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCgg                     >  1:1950651/1‑151 (MQ=255)

                                                                                                                                      aaT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCG‑ATCCTCCCGGATGCACCATATTCTCCGAAATATTCAGAAATGAAGGTATCGAAGAGGTAGCgg                     >  2:541069/1‑150 (MQ=255)

                                                                                                                                       aaaaatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGAAATCTTCAGCAATGAACGTATagaagatctaccg                      >  6:1137853/5‑144 (MQ=11)

                                                                                                                                        aaaatCAGCAGTAAACTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCgg                     <  5:1138636/148‑1 (MQ=255)

                                                                                                                                        aaaatCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCgg                     <  3:911527/148‑1 (MQ=255)

                                                                                                                                          aatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTTCATATTCCGCTATGATCGTATCGAATATGTacagga                    >  1:1364913/2‑145 (MQ=1)

                                                                                                                                          aatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGGATGGAAGAGGTacagga                    >  1:1098478/2‑145 (MQ=14)

                                                                                                                                           atCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATCCGAAATATTCAGATATTAACGTATGGAAGATGTAGAGTAAt                  >  2:2058513/1‑151 (MQ=11)

                                                                                                                                           atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcc                                                                                >  1:1443484/1‑89 (MQ=40)

                                                                                                                                           atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcc                                                                                <  2:1443051/89‑1 (MQ=40)

                                                                                                                                           atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAAt                  >  7:1730292/1‑151 (MQ=255)

                                                                                                                                            tctgcagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGAAATATGCAGGTATGAAGGGATGCAAGAGCAAGAGGAAt                  >  1:484738/4‑150 (MQ=11)

                                                                                                                                             ctgcagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATGTTACGAAATATTCCGAAATTAACGTATCGAAAAggacgagtaata                 >  1:178677/3‑140 (MQ=11)

                                                                                                                                             cagcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAAtt                 >  5:1598583/1‑150 (MQ=255)

                                                                                                                                             cagcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCGAAGAGATAg                        >  8:2183040/1‑143 (MQ=255)

                                                                                                                                              tgcagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATCCGAAATCTTCAGAAATGAAGGTATAGAAGAGGt                          >  2:1296554/2‑140 (MQ=14)

                                                                                                                                              tgcagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATCCGAAATATTCAGCAATTAACGTATAGaatataaatagtaac                  >  6:2153685/2‑135 (MQ=11)

                                                                                                                                              agcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTa                >  4:1805991/1‑150 (MQ=255)

                                                                                                                                               gcagTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTaa               <  2:2278493/150‑1 (MQ=255)

                                                                                                                                               gcagTAAATTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAAc              <  2:1307367/151‑1 (MQ=255)

                                                                                                                                               gcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAAc              >  6:359639/1‑151 (MQ=255)

                                                                                                                                                cagTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCGAGATCTGGAGGTAGTAAGGGATGGAAGAGGTAGAGGAATTAAcc             >  5:57063/1‑151 (MQ=14)

                                                                                                                                                 agTAAATTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCg            <  8:1111606/151‑1 (MQ=255)

                                                                                                                                                 agTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCg            >  3:416996/1‑151 (MQ=255)

                                                                                                                                                  gTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGATGGTATCGAAGAGGTAGCGGAATTAACCGc           <  4:549998/151‑1 (MQ=255)

                                                                                                                                                   tAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCGAGATCTGCAGGTATGAAGGGTTGGAAGAGGTAGCGGAATTTACCgtt          >  5:1369055/1‑149 (MQ=14)

                                                                                                                                                    aaatTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAg         <  4:204048/151‑1 (MQ=255)

                                                                                                                                                    aaagTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAg         <  4:45886/151‑1 (MQ=255)

                                                                                                                                                    aaagTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATGTTACTTAATATTCAGGTATTATGGTTTGCAATAGGTAGCGGAATaaaaagaac         >  1:2273918/1‑145 (MQ=11)

                                                                                                                                                    aaaaTAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGtaacc                                                    >  4:297031/1‑106 (MQ=17)

                                                                                                                                                       ataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGc                                                                                                  <  3:1383896/59‑1 (MQ=37)

                                                                                                                                                         aa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACt    >  3:1562532/1‑151 (MQ=40)

                                                                                                                                                            atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACTAg  <  3:1241320/151‑1 (MQ=40)

AGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGA‑ACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACTAG  >  DH1/1065469‑1065768

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