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%	1450.1 / 10.0	374	Q175E (CAG→GAG)	tonB	protein TonB
Reads supporting (aligned to +/- strand): ref base G (4/4); new base C (180/186); total (184/190)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	DH1	2,568,067	0	C	T	98.1%	1431.5 / 5.4	374	G174G (GGG→GGA)	tonB	protein TonB
Reads supporting (aligned to +/- strand): ref base C (3/4); new base T (181/186); total (184/190)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

AATCCCACTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTATTT  >  DH1/2567917‑2568211

||

      aCTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTc                                                                                                                                             >  7:599527/1‑150 (MQ=255)

      aCTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTc                                                                                                                                             >  3:507925/1‑150 (MQ=255)

      aCTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCa                                                                                                                                            <  5:1311411/151‑1 (MQ=255)

        tGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAAt                                                                                                                                          >  1:1422138/1‑151 (MQ=255)

         gCCTGGCTTACCCGGCTCATACAGCCATCTGCGCATCGCAGTTTTCACCTCACGCTCAAACATGTGCTCAGGCTTGTCTGAGAGGATTTGTACGGTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCATTCAAt                                                                                                                                          <  4:66546/150‑1 (MQ=255)

         gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCaa                                                                                                                                           <  3:353685/149‑1 (MQ=255)

         gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAAt                                                                                                                                          <  6:1154922/150‑1 (MQ=255)

         gCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAAt                                                                                                                                          <  7:473238/150‑1 (MQ=255)

          ccTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgc                                                                                                                                        >  7:1512751/1‑151 (MQ=255)

            tGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTAACCTCACGCTCAAACATGGTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgcg                                                                                                                                       >  5:399961/1‑150 (MQ=255)

             ggCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgcg                                                                                                                                       <  1:164707/149‑1 (MQ=255)

             ggCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATgcg                                                                                                                                       <  1:500224/149‑1 (MQ=255)

              gCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCa                                                                                                                                     >  2:658722/1‑150 (MQ=255)

                ttACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACTCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATg                                                                                                                                  >  5:1470396/1‑151 (MQ=255)

                 tACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATg                                                                                                                                  <  2:951588/150‑1 (MQ=255)

                 tACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGc                                                                                                                                 >  8:25475/1‑151 (MQ=255)

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

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

                     cGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCtgt                                                                                                                             >  7:94/1‑151 (MQ=255)

                        cTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGCTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCt                                                                                                                          >  5:527315/1‑151 (MQ=255)

                         tCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTc                                                                                                                         >  6:457693/1‑151 (MQ=255)

                          cATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCg                                                                                                                        <  6:209139/151‑1 (MQ=255)

                           aTAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGt                                                                                                                       <  8:1459855/151‑1 (MQ=255)

                                  cATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGAt                                                                                                                >  8:1257923/1‑151 (MQ=255)

                                   aTCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGAt                                                                                                                >  3:1326858/1‑150 (MQ=255)

                                   aTCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGAt                                                                                                                >  8:482312/1‑150 (MQ=255)

                                      tGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTg                                                                                                            >  7:840004/1‑151 (MQ=255)

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

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

                                       gcgcATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGc                                                                                                           <  2:871934/151‑1 (MQ=255)

                                           atcgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCt                                                                                                       <  7:289212/151‑1 (MQ=255)

                                           atcgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCt                                                                                                       >  5:951508/1‑151 (MQ=255)

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

                                             cgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTg                                                                                                      >  7:414204/1‑150 (MQ=255)

                                             cgcatTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGa                                                                                                     <  4:1075054/151‑1 (MQ=255)

                                                atTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGAtt                                                                                                   >  4:144834/1‑150 (MQ=255)

                                                atTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGAtt                                                                                                   >  7:987024/1‑150 (MQ=255)

                                                   tttCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACgg                                                                                               >  3:88522/1‑151 (MQ=255)

                                                    ttCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACgg                                                                                               <  6:535836/150‑1 (MQ=255)

                                                    ttCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGc                                                                                              >  8:306580/1‑151 (MQ=255)

                                                     tCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGc                                                                                              <  1:831771/150‑1 (MQ=255)

                                                     tCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCt                                                                                             >  2:1382440/1‑151 (MQ=255)

                                                      cACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCtt                                                                                            <  7:1148829/151‑1 (MQ=255)

                                                      cACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTAAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCt                                                                                             <  4:37844/150‑1 (MQ=255)

                                                       aCCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCtt                                                                                            >  2:46048/1‑150 (MQ=255)

                                                        ccTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTaa                                                                                          <  4:1968/151‑1 (MQ=255)

                                                           cACGCTCAAACATGTTCGCATGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATg                                                                                        <  5:599365/150‑1 (MQ=255)

                                                           cACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAAt                                                                                         <  5:1188000/149‑1 (MQ=255)

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

                                                              gCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGcgcg                                                                                    <  4:348561/151‑1 (MQ=255)

                                                              gCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGcgcg                                                                                    <  2:1232040/151‑1 (MQ=255)

                                                                 cAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTgg                                                                                 >  7:343317/1‑151 (MQ=255)

                                                                   aaCATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTc                                                                               >  7:696978/1‑151 (MQ=255)

                                                                   aaCATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTc                                                                               <  1:957364/151‑1 (MQ=255)

                                                                    aCATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTcc                                                                              >  2:1606800/1‑151 (MQ=255)

                                                                    aCATGTACTCAGGCTTGGCTGAGAGGATTTGTAAGTTATCCACGCGACGATCCGGCGTAAGGTCAAATTTAACTTTAACCTCTCCTTGAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGtac                                                                              <  4:1297823/151‑3 (MQ=255)

                                                                     cATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCt                                                                             <  3:1040102/151‑1 (MQ=255)

                                                                     cATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCt                                                                             <  1:71196/151‑1 (MQ=255)

                                                                         ttCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCCCGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCGTAATGCGCGTGGTCCTGacg                                                                         >  8:1632784/1‑149 (MQ=255)

                                                                           cGCAGGCTTGGCTGAGAGGATTTGTACGTTATCAACGCGACCATCCGGCGTAACGTCAAATTTAAGTTTACCCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGcc                                                                       >  1:1197741/1‑151 (MQ=255)

                                                                             cAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCcac                                                                     >  5:1662828/1‑151 (MQ=255)

                                                                               ggCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCacac                                                                   <  6:338854/151‑1 (MQ=255)

                                                                                gCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACt                                                                  >  7:522021/1‑151 (MQ=255)

                                                                                     gCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTaa                                                             <  4:1217653/151‑1 (MQ=255)

                                                                                       tGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAAcc                                                           >  2:1654974/1‑151 (MQ=255)

                                                                                       tGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACAATAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAAcc                                                           >  4:1598541/1‑151 (MQ=255)

                                                                                        gagaGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAAcc                                                           <  6:159621/150‑1 (MQ=255)

                                                                                        cagagGCTTTGTACGTTATCCACGCGACCATCCCGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAAcc                                                           <  2:1127148/149‑1 (MQ=255)

                                                                                           agGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCt                                                       <  8:557626/151‑1 (MQ=255)

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

                                                                                             gATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGTCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCtt                                                      <  6:890678/150‑1 (MQ=255)

                                                                                                   tACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGtt                                               >  8:824010/1‑151 (MQ=255)

                                                                                                   tACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGCGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGtt                                               >  5:750112/1‑151 (MQ=255)

                                                                                                    aCGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGtt                                               <  3:1170509/150‑1 (MQ=255)

                                                                                                    aCGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGtt                                               <  1:318419/150‑1 (MQ=255)

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

                                                                                                      gTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGc                                             >  6:1019320/1‑150 (MQ=255)

                                                                                                         aTCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGc                                          <  3:548056/150‑1 (MQ=255)

                                                                                                             aCGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTg                                     >  3:1635902/1‑151 (MQ=255)

                                                                                                             aCGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTg                                     <  3:431690/151‑1 (MQ=255)

                                                                                                                cgACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAATCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTg                                  >  5:1032597/1‑151 (MQ=255)

                                                                                                                cgACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCATGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTg                                  >  2:1470361/1‑151 (MQ=255)

                                                                                                                 gcccATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGt                                 >  6:698786/3‑151 (MQ=255)

                                                                                                                 gACCATCCGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGt                                 >  5:1579670/1‑151 (MQ=255)

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

                                                                                                                       ccGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTg                            <  5:416042/150‑1 (MQ=255)

                                                                                                                        cGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGAt                          <  3:948629/151‑1 (MQ=255)

                                                                                                                        cGGCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGAt                          >  6:1210537/1‑151 (MQ=255)

                                                                                                                          gCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTTGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGt                        <  3:1215219/151‑1 (MQ=255)

                                                                                                                          gCGTAACGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGt                        >  7:1378767/1‑151 (MQ=255)

                                                                                                                               aCGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGc                   >  4:7987/1‑151 (MQ=255)

                                                                                                                               aCGTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGc                   <  8:1059545/151‑1 (MQ=255)

                                                                                                                                 gTCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgt                 <  3:758832/151‑1 (MQ=255)

                                                                                                                                  tCAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtg                >  8:643651/1‑151 (MQ=255)

                                                                                                                                   cAAATTTAACTTTAACCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCgtgc               <  1:1117978/151‑1 (MQ=255)

                                                                                                                                    aaaTTTAACTTTAACCTCTCCTTCAATGCGCAATGACTGTGCTCGTGCCGTATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGAAGCCGTTGCTGTACTTGATGTCAGGCgtgcc              >  8:1564196/1‑151 (MQ=255)

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

                                                                                                                                             tttaacCTCTCATTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGAGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGt      >  6:961014/1‑150 (MQ=255)

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

                                                                                                                                                aacCTCTCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTAtt   <  7:182786/150‑1 (MQ=255)

||

AATCCCACTGCCTGGCTTACCCGGCTCATAACGCCATCTGCGCATCGCATTTTTCACCTCACGCTCAAACATGTTCGCAGGCTTGGCTGAGAGGATTTGTACGTTATCCACGCGACCATCCGGCGTAACGTCAAATTTAACTTTAACCTGCCCTTCAATGCGCAATGCCTGTGCTCGTGCCGGATACTGCGGCTGATTACGGCTTAATGCGCGTGGTCCTGAAGCCACACTGGTAACCGGCTTGCTGGTTGCAGCCGTTGCTGTACTTGATGTCAGGCGTGCCGGTGCCGTATTT  >  DH1/2567917‑2568211

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