BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	DH1	2,568,066	2 bp→CT	coding (522‑523/720 nt)	tonB ←	protein TonB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	2,568,066	0	G	C	97.9%	1459.0 / 12.0	376	Q175E (CAG→GAG)	tonB	protein TonB
Reads supporting (aligned to +/- strand): ref base G (6/2); new base C (190/178); total (196/180)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
*	DH1	2,568,067	0	C	T	97.9%	1413.9 / 9.6	376	G174G (GGG→GGA)	tonB	protein TonB
Reads supporting (aligned to +/- strand): ref base C (6/2); new base T (188/177); total (196/180)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 376 total aligned reads displayed.

CCACTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTATTTTC  >  DH1/2567921‑2568213

||

     gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATg                                                                                                                                           >  5:279976/1‑151 (MQ=255)

     gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATg                                                                                                                                           >  1:767310/1‑151 (MQ=255)

     gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATg                                                                                                                                           >  8:1134903/1‑151 (MQ=255)

      ccTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgc                                                                                                                                          >  5:1560874/1‑151 (MQ=255)

       cTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgcg                                                                                                                                         >  8:616523/1‑151 (MQ=255)

         ggCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgcgc                                                                                                                                        >  6:189997/1‑150 (MQ=255)

            ttACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATg                                                                                                                                    >  7:776855/1‑151 (MQ=255)

             tACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGc                                                                                                                                   >  7:552474/1‑151 (MQ=255)

              aCCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGcc                                                                                                                                  >  7:365097/1‑151 (MQ=255)

               cccGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCt                                                                                                                                 <  7:1324528/151‑1 (MQ=255)

                ccGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTACCTTTCACCTCTCCGTCAATGCGCAATGCCtg                                                                                                                                >  4:591509/1‑151 (MQ=255)

                ccGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCtg                                                                                                                                <  1:1002159/151‑1 (MQ=255)

                    cTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCt                                                                                                                            >  2:616166/1‑151 (MQ=255)

                    cTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCt                                                                                                                            >  1:365419/1‑151 (MQ=255)

                       aTAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCTAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGt                                                                                                                         >  6:100467/1‑151 (MQ=255)

                       aTAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGt                                                                                                                         <  6:975749/151‑1 (MQ=255)

                            gCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCgg                                                                                                                    >  1:822902/1‑151 (MQ=255)

                              cATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGAt                                                                                                                  >  6:181605/1‑151 (MQ=255)

                              cATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCCGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGAt                                                                                                                  <  3:1354873/151‑1 (MQ=255)

                                tCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATAc                                                                                                                >  3:764441/1‑151 (MQ=255)

                                 cTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGCCCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACt                                                                                                               >  5:813200/1‑151 (MQ=255)

                                   gcgcATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATa                                                                                                                 <  2:1432273/147‑1 (MQ=255)

                                   gcgcATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACt                                                                                                               <  4:446992/149‑1 (MQ=255)

                                   gcgcATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTg                                                                                                              <  8:177718/150‑1 (MQ=255)

                                    cgcatcgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCg                                                                                                            <  5:702274/151‑1 (MQ=255)

                                    cgcatcgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCg                                                                                                            >  5:201453/1‑151 (MQ=255)

                                         cgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTg                                                                                                        >  5:4291/1‑150 (MQ=255)

                                         cgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTg                                                                                                        <  1:1518071/150‑1 (MQ=255)

                                         cgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTg                                                                                                        <  2:536495/150‑1 (MQ=255)

                                           catTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGAtt                                                                                                     <  6:120878/151‑1 (MQ=255)

                                            atTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAAGCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGAtt                                                                                                     >  4:76739/1‑150 (MQ=255)

                                            atTTTTCACCTCACGCTCAAACATGTTCGAAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTa                                                                                                    <  4:1449965/151‑1 (MQ=255)

                                                 tCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTAGTGCCGGATACTGAGGCTGATTACGGCt                                                                                               >  7:149468/1‑151 (MQ=255)

                                                  cACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCtt                                                                                              >  7:1124542/1‑151 (MQ=255)

                                                  cACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCt                                                                                               <  8:852668/150‑1 (MQ=255)

                                                      tCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATg                                                                                          >  3:1039326/1‑151 (MQ=255)

                                                       cACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATgc                                                                                         >  3:493871/1‑151 (MQ=255)

                                                       cACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATgc                                                                                         <  8:472696/151‑1 (MQ=255)

                                                         cGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATgcgc                                                                                       <  5:1073648/151‑1 (MQ=255)

                                                            tCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTg                                                                                    >  2:955568/1‑151 (MQ=255)

                                                             cAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTgg                                                                                   <  2:758579/151‑1 (MQ=255)

                                                             cAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTg                                                                                    <  7:249017/150‑1 (MQ=255)

                                                                aCATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTcc                                                                                <  7:656871/151‑1 (MQ=255)

                                                                  aTGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTg                                                                              >  8:386278/1‑151 (MQ=255)

                                                                      tCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAg                                                                           <  3:1180335/150‑1 (MQ=255)

                                                                        gCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCca                                                                        <  1:1320079/151‑1 (MQ=255)

                                                                        gCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCca                                                                        <  5:628184/151‑1 (MQ=255)

                                                                         cAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCcac                                                                       >  8:1333795/1‑151 (MQ=255)

                                                                          aGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAAGCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCcac                                                                       >  2:1087709/1‑150 (MQ=255)

                                                                           ggCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCacac                                                                     <  7:597061/151‑1 (MQ=255)

                                                                             gttGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCGTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTATTGCGGGTGGTCCTGATGCCACACTg                                                                   <  3:1482313/150‑1 (MQ=255)

                                                                               tGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTgg                                                                  <  4:173911/150‑1 (MQ=255)

                                                                               tGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTg                                                                   <  2:1278559/149‑1 (MQ=255)

                                                                               tGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGt                                                                 >  6:29667/1‑151 (MQ=255)

                                                                                ggCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGt                                                                 >  5:744829/1‑150 (MQ=255)

                                                                                 gCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTaa                                                               <  2:1540358/151‑1 (MQ=255)

                                                                                 gCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTa                                                                >  5:1177546/1‑150 (MQ=255)

                                                                                     agagGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCgg                                                           <  1:382231/151‑1 (MQ=255)

                                                                                       atgATTTGTCCTTTATCCACGCGCCCATCCGGATTAACGTCAAATTTAACTTTAACCTCTCCTTCCATGCGCAATGCCTGTGCGCGTGCCGGATACTGCGGCGGATTACGGCTTAATGCGCGTGGGCCTGACGCCACACTGCTACCCGGc                                                          <  4:979819/148‑1 (MQ=255)

                                                                                        ggATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCtt                                                        >  1:1122821/1‑151 (MQ=255)

                                                                                           tttGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCTGATACTGGGGCTGATTACGGCTTAATGCGCGTGGTCCTGAATCCACATTGGTAACCGGCTTGCt                                                     >  3:849873/1‑151 (MQ=255)

                                                                                            ttGTACGTTATCGACGCGACCATCCGGGGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGGAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATCACGGCTTATTGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCt                                                     <  6:1294628/150‑1 (MQ=255)

                                                                                              gTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGt                                                  >  3:259649/1‑151 (MQ=255)

                                                                                                aCGTTATCCACGCGACCATCCGGCGTAACGTCTAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTg                                                <  8:736908/151‑1 (MQ=255)

                                                                                                aCGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTg                                                >  1:1418489/1‑151 (MQ=255)

                                                                                                 cGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGc                                               >  4:869750/1‑151 (MQ=255)

                                                                                                  gTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGc                                               >  7:461111/1‑150 (MQ=255)

                                                                                                  gTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGAGCAATGCCTGTGCTCGTGCCGGATACTGCGGCAGATTACGGCTTAATGCGCGTCGTCCTGAAGCCACACTGCTAACCGGCTTGCTCGttcc                                               >  6:780437/1‑148 (MQ=255)

                                                                                                     aTCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGcc                                           >  5:334769/1‑151 (MQ=255)

                                                                                                       ccACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCg                                          <  6:1480942/150‑1 (MQ=255)

                                                                                                       ccACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCATGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGt                                         <  1:581410/151‑1 (MQ=255)

                                                                                                          cgcgACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGt                                                  <  7:1381175/139‑1 (MQ=255)

                                                                                                          cgcgACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGc                                      <  3:1432362/151‑1 (MQ=255)

                                                                                                           gcgACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCt                                     >  4:326421/1‑151 (MQ=255)

                                                                                                             gACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGt                                   <  8:322212/151‑1 (MQ=255)

                                                                                                              aCCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTa                                  <  4:1185194/151‑1 (MQ=255)

                                                                                                               ccATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTAc                                 <  1:1183192/151‑1 (MQ=255)

                                                                                                                cATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACt                                >  5:1323946/1‑151 (MQ=255)

                                                                                                                 aTCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACtt                               <  7:975314/151‑1 (MQ=255)

                                                                                                                      gCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGt                          >  5:1351129/1‑151 (MQ=255)

                                                                                                                          aaCGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTCATGTCAgg                      >  4:1285792/1‑151 (MQ=255)

                                                                                                                            cGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCg                    >  7:1150714/1‑151 (MQ=255)

                                                                                                                            cGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCg                    <  1:991368/151‑1 (MQ=255)

                                                                                                                             gTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgt                   >  6:872442/1‑151 (MQ=255)

                                                                                                                               cAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgc                 >  6:663616/1‑151 (MQ=255)

                                                                                                                               cAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtg                  <  2:261366/150‑1 (MQ=255)

                                                                                                                                aaaTTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATAGTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgcc                >  6:1051793/1‑151 (MQ=255)

                                                                                                                                 aaTTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGAAGCCGTTGCTGTACTTGATGTCAAGCgtgccg               >  6:440591/1‑151 (MQ=255)

                                                                                                                                  aTTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccgg              <  3:570731/151‑1 (MQ=255)

                                                                                                                                   tttaactttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccggt             >  3:975006/1‑151 (MQ=255)

                                                                                                                                     taactttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccggtgc           >  3:1311606/1‑151 (MQ=255)

                                                                                                                                      aactttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccggtgc           <  3:923214/150‑1 (MQ=255)

                                                                                                                                       actttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccggtgccg         <  5:1458040/151‑1 (MQ=255)

                                                                                                                                       actttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgccggtgccg         >  1:54761/1‑151 (MQ=255)

                                                                                                                                         tttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTATTGCGCGTCGTCCTGAAGCCACACTGGTAACCGGCTTGAAGCTAGCACCCGTTGCTGTACTTGATATCAGGAgtgccggtgacga        >  8:1519542/1‑149 (MQ=255)

                                                                                                                                          ttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGt        <  4:1134607/149‑1 (MQ=255)

                                                                                                                                          ttaacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGt        <  8:1407699/149‑1 (MQ=255)

                                                                                                                                            aacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTAttt    >  7:831048/1‑151 (MQ=255)

                                                                                                                                            aacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTAtt     <  2:368498/150‑1 (MQ=255)

                                                                                                                                              ccTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGGTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTATTTTc  <  2:872097/151‑1 (MQ=255)

||

CCACTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTATTTTC  >  DH1/2567921‑2568213

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