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	24.2%	264.2 / 98.4	174	intergenic (+53/‑10)	ECDH1_RS05245/ECDH1_RS05250	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (70/43); new base G (33/3); total (104/70)
Fisher's exact test for biased strand distribution p-value = 7.13e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-02

Only 100 of 200 total aligned reads displayed.

AAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACT  >  DH1/1065472‑1065766

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

 aagtaagtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATATTACTAGATATGGCTGTAGTATGGTTTG‑G‑‑AAT‑‑ATGTGCAGTAATATAA‑Attagg                                                                                                                                                 >  1:249354/1‑148 (MQ=14)

      aacagccctttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATg                                                                                                                                            <  1:1365961/141‑1 (MQ=255)

          gccatttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATg                                                                                                                                            <  8:1474586/141‑1 (MQ=255)

          gccatttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATCGTACTTGATATGGGGTAAGTAGCGGTGT‑C‑‑AAT‑‑TTAAGCAGTAAAAAAA‑AAATCCCCATGcg                                                                                                                                          >  7:1474432/7‑148 (MQ=14)

             cGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACTAAATATGGATCTATTAGC‑GTATAG‑‑AAT‑‑ATCAACAGTAATAAAA‑ATAAGCCGATACAta                                                                                                                                        >  8:255344/1‑147 (MQ=2)

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

                          gCTCAGCTGGATAGAGTACTCGGCTACGAACCGCG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGAt                                                                                                                        <  7:1074744/151‑1 (MQ=255)

                                           aCTCGGCTACGAACCGTT‑GGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTAATTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAGAATAAGTTTTCCCGATGCATCCGTa                                                                                                                                    <  5:402910/123‑1 (MQ=11)

                                            cTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACTAGATCTGGCTGAAGTATCGGTAT‑C‑‑AAT‑‑ATCTGCAGGAATTAAA‑ATTTCGCCATTCCTCGGAACCGCAGCATGATAGACTACTCCGCTACGac                                                                                                      >  5:1481501/1‑150 (MQ=9)

                                             tCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAGAATAA‑GTTTtccccttgca                                                                                                                                         >  1:826596/1‑108 (MQ=17)

                                                   cactaCCGAGCCGCTCGGAGCTTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGc                                                                                                <  7:1278808/147‑1 (MQ=14)

                                                    cGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCgg                                                                                              >  8:325851/1‑151 (MQ=255)

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

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

                                                               ggTCGGAGGTTCGAATCCTCCCGGATGCACCATATCGTACTTCATATGGCTGTAGTAGCGGTAG‑C‑‑AAG‑‑ATCTGGAGGAAAATAA‑ATTACCCGATGCGTAGGTAGGACAGATCGATAGAGTAATCGGATACGAACCGAGCGGTCGGAGGTTc                                                                                    >  1:173108/1‑151 (MQ=17)

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

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

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

                                                                           caaTCCTCCGGGAATCACGATCTTTGAGTTGAGAGGGGTTTACTAGCGATCT‑C‑‑AAT‑‑ATCTGCAGTAAACTAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCgg                                                                        <  5:670616/150‑1 (MQ=11)

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

                                                                            aagcctccCGGATCCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑A‑‑AAT‑‑ATCAGCAGTAAACTAA‑CTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCg                                                                         <  7:1000131/146‑1 (MQ=255)

                                                                            aaTCCTCCCTGATGCAAGATCTCTTACTTTATATCGCTTTAGGAGCTGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCgg                                                                        <  7:1445061/150‑1 (MQ=14)

                                                                              ccctccCGGATGGACTATCTTTTTCTTGATATGTATTTAGTAGCTGTAT‑C‑‑AAT‑‑ATCTGCAGTAAAATAA‑TTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATg                                                                     <  3:1051690/150‑1 (MQ=11)

                                                                               cctccCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGc                                                                    >  7:1249250/1‑151 (MQ=39)

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

                                                                                ctccCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAA‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCAGGATGc                                                                    >  1:47439/1‑150 (MQ=255)

                                                                                    cGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCat                                                               >  8:1064088/1‑151 (MQ=39)

                                                                                          cACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTACCCATTCATACGTACCTGATCTTGATAGAGTAATACGCTACGAACCCACCCGTGCGATGTTAGAATCCTGCCCGaatcaacaaaatat                                                          >  8:310017/1‑137 (MQ=2)

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

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

                                                                                                  ctTACTTTATATCTCTTTAGTAGCTCTAT‑A‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGaaa                                                    <  1:359317/148‑1 (MQ=14)

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

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

                                                                                                     aCTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATGCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACGGAGCGGACCCCGGTTCAAATCATTCCCGAAGATCCATATtctaagccctcttct                                              >  1:851580/1‑139 (MQ=14)

                                                                                                       ttGATACGGCTTTAGTAGCGGTAT‑CAAAAA‑‑ATCTGCGGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTAGATCTTcc                                              >  1:1096613/1‑150 (MQ=14)

                                                                                                       ttGATACGGCTTTAGTAGCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCatataatactaaatattc                                               >  8:602556/1‑136 (MQ=17)

                                                                                                        tGATACGGCTTTAGTAGCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATACTACTAAATATTc                                               >  7:720576/1‑149 (MQ=14)

                                                                                                            atGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATg                                        >  7:691510/1‑150 (MQ=255)

                                                                                                              ggCTTTAGTAGCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACTATATTCTCCGAGATCTTCatctagt                                        >  7:1123407/1‑144 (MQ=11)

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

                                                                                                                 tttAGTAGCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTACGAAATATTCAGCAATTAAgg                                    >  6:87469/1‑151 (MQ=11)

                                                                                                                 tttAGTAGCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCAAATCCTCCCGGATGCACCATATTATACTAAATATTCATAAATTAAgg                                    >  4:975218/1‑151 (MQ=14)

                                                                                                                 tttAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAAACTTCAGCAATGAAGGtt                                  >  5:1335511/1‑150 (MQ=25)

                                                                                                                  ttAGTAGCTTTAT‑A‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTa                                  <  7:614797/150‑1 (MQ=255)

                                                                                                                  ttAGTACCGTTAT‑A‑‑AAT‑‑ATAAGCATTAAAATAT‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGTAACCTTCAGCAATGAAGGTa                                  <  6:1175978/150‑1 (MQ=11)

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

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

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

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

                                                                                                                       agCGGTAT‑CAAAAA‑‑ATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGATATCTTCATCTATGAAGGtttgg                               >  3:1326945/1‑146 (MQ=11)

                                                                                                                        gCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGa                              >  5:642906/1‑148 (MQ=255)

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

                                                                                                                        gCGGTAT‑C‑‑AAAAAATCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATCCTTAATCTGCATAAATGAAGGTaagga                              >  8:408371/1‑146 (MQ=11)

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                    at‑‑atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTAAGAAATCTTAAGCAATGAAGGTATCGAAGAGGTAgagt                   >  8:1193609/1‑147 (MQ=255)

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

                                                                                                                                     t‑‑atCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAACCTTCAGCAATGAAGGTATCGAAGAGATaagctac                 >  7:343083/1‑145 (MQ=255)

                                                                                                                                       aatCTGCAGTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAGATATTCAGCTAGTACCGGATCGAAGATGTCGCGGa                  >  1:285315/2‑150 (MQ=11)

                                                                                                                                       aatCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTc                                                                                                               <  5:1307963/56‑1 (MQ=255)

                                                                                                                                        atttccattAAATTAT‑ATTTTCCGATGTATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGaa                 <  6:294695/142‑1 (MQ=14)

                                                                                                                                        attagcataaaagTAT‑ATTTTCCGATGTATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGaa                 <  2:996570/141‑1 (MQ=14)

                                                                                                                                         tCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCt                                                                                                           <  1:528452/59‑1 (MQ=255)

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

                                                                                                                                         tCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAg                                                                                        >  6:167109/1‑78 (MQ=255)

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

                                                                                                                                          cagcagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCCGTCGGAGGTTCGAATAGTCGCGGATGCACAATATTCTCCGAAATCTGCAGGATGGAGGGTATCGATGAGGTTGCGGAtttg              >  3:225165/1‑150 (MQ=14)

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

                                                                                                                                             cagTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGa                                                                                         <  7:842019/73‑1 (MQ=38)

                                                                                                                                             cactatactaa‑ctTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAAcc           <  3:127623/143‑1 (MQ=255)

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

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

                                                                                                                                               gTAAAGTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGa                  <  6:237430/142‑1 (MQ=255)

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

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

                                                                                                                                                 aaaaTAA‑GTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCa        <  6:284556/150‑1 (MQ=255)

                                                                                                                                                  aattaa‑gtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCg          <  6:107636/144‑1 (MQ=255)

                                                                                                                                                  aagtaa‑gtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTACGAAATCTTCAGGAATGAAGGGTTCGAATAGCAAGCGGAATTAAcctagac      >  5:107868/4‑146 (MQ=11)

                                                                                                                                                  aaataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCgg                                                                                              >  8:639726/1‑63 (MQ=39)

                                                                                                                                                  aaataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCg                                                                                           >  4:934714/1‑66 (MQ=39)

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

                                                                                                                                                   attaa‑ttTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAgg      <  6:1209175/148‑1 (MQ=255)

                                                                                                                                                   agtaa‑gtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGATATCTTCAGGAATGACGGTATCGAAGAGGTAGAGTAATTAACCGctgc      >  5:1208217/3‑147 (MQ=14)

                                                                                                                                                   aataa‑ttTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGc     <  1:430408/151‑1 (MQ=255)

                                                                                                                                                   aataa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCGAAGAGATAGCGGAATTAACCGCAGGc     >  6:374924/1‑151 (MQ=25)

                                                                                                                                                     taa‑gtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATGCTCCGAAATATGCAGCAAGGATGGGATCGAAGAGCTAGCGGAAATAACAGCACGccc   >  3:335910/1‑149 (MQ=11)

                                                                                                                                                     taa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCAc   >  1:622723/1‑151 (MQ=40)

                                                                                                                                                      aa‑atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACt  <  1:949698/151‑1 (MQ=40)

                                                                                                                                                         atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAgg                                                                                       >  8:587790/1‑64 (MQ=31)

                                                                                                                                                         atTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAgg                                                                                       <  7:587798/64‑1 (MQ=31)

AAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAG‑CGGTCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTAT‑C‑‑AAT‑‑ATCAGCAGTAAAATAA‑ATTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATTAACCGCAGGCACT  >  DH1/1065472‑1065766

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