BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	2,595,851	0	A	G	52.5%	55.7 / inf	425	pseudogene (2764/3744 nt)	narZ	nitrate reductase
Reads supporting (aligned to +/- strand): ref base A (98/103); new base G (105/117); total (204/221)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-02

Only 100 of 430 total aligned reads displayed.

AGGCTTTCACGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCGCCTTT  >  DH1/2595701‑2596001

aGGCTTTCACGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTa                                                                                                                                                        <  7:325521/151‑1 (MQ=255)

     ttCACGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGCCCTTCCGCTGTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGgtgt                                                                                                                                                   >  3:692650/1‑151 (MQ=255)

     ttCACGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGgtg                                                                                                                                                    <  8:330624/150‑1 (MQ=255)

        aCGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTcc                                                                                                                                                <  1:1462825/151‑1 (MQ=255)

        aCGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTcc                                                                                                                                                <  7:415785/151‑1 (MQ=255)

          ggCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAg                                                                                                                                              <  2:1525301/151‑1 (MQ=255)

          ggCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAg                                                                                                                                              <  7:182150/151‑1 (MQ=255)

          ggCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGGGTAGTTGAGCTTACGCAGCAGCTCCATCTCGCTCTGGGTGTTcc                                                                                                                                                >  8:215117/1‑149 (MQ=255)

            cTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGc                                                                                                                                            >  5:284810/1‑151 (MQ=255)

                cTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATc                                                                                                                                        >  1:46863/1‑151 (MQ=255)

                  gACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATccc                                                                                                                                      <  5:1185663/151‑1 (MQ=255)

                         ggTTTCCGGTGCCAGTGTCAGGAGCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTAcc                                                                                                                                <  1:536960/150‑1 (MQ=255)

                             tCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCa                                                                                                                            <  5:625449/150‑1 (MQ=255)

                             tCCGCTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCATTTCGCCGGACCTTCCGCATTGGTGTAGTTGAGCATACGCAGCAGATCCATCTCCCTCTGGGTGTTCCAGGCAATCCCATTACCGCCAt                                                                                                                           >  4:1094695/1‑151 (MQ=255)

                              cgggTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCAtt                                                                                                                          <  2:1594797/149‑1 (MQ=255)

                              ccGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCAtt                                                                                                                          <  3:607104/151‑1 (MQ=255)

                               cGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTGCCGCCATTa                                                                                                                         <  1:1267332/151‑1 (MQ=255)

                                ggTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGTAATCCCTTTACCGCCATTAc                                                                                                                        >  1:1705754/1‑151 (MQ=255)

                                   gCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCg                                                                                                                      <  6:1713976/150‑1 (MQ=255)

                                   gCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGa                                                                                                                     <  5:1398310/151‑1 (MQ=255)

                                     cAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGAtt                                                                                                                   <  1:1130777/151‑1 (MQ=255)

                                     cAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGCCCTTTCGCCGGACCAACCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCGCTGGGTGTTCCAGGAACTCCCTTTACCGCCATTAGCGagt                                                                                                                   >  4:734546/1‑149 (MQ=255)

                                      aGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGAttt                                                                                                                  >  3:891175/1‑151 (MQ=255)

                                             ggATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCAt                                                                                                           <  1:892022/151‑1 (MQ=255)

                                             ggATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCAt                                                                                                           <  2:55275/151‑1 (MQ=255)

                                                     tctGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGc                                                                                                   >  4:901599/1‑151 (MQ=255)

                                                      ctGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCgg                                                                                                     <  3:1843821/148‑1 (MQ=255)

                                                      ctGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGACTTTCTCCATCAGCggg                                                                                                    <  4:643666/149‑1 (MQ=255)

                                                       tGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGccc                                                                                                 >  5:804616/1‑150 (MQ=255)

                                                        gCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCg                                                                                                 >  8:899041/1‑150 (MQ=255)

                                                        gCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCg                                                                                                 >  1:322294/1‑150 (MQ=255)

                                                          cGCATCAATGTCGGTGATCAGCCTCGTCTGGCCTTTCGCCTGACCTTACTCTTTGTTGTAGGTGAGCTTACGCAGCATATCCATCTCGCTCTAGGTGTTCCAGTCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATa                                                                                              <  6:947569/151‑1 (MQ=255)

                                                              tCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAgg                                                                                          >  6:1153553/1‑151 (MQ=255)

                                                                aaTTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGt                                                                                         >  1:157669/1‑150 (MQ=255)

                                                                 aTTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTaa                                                                                       <  7:505666/151‑1 (MQ=255)

                                                                 aTTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTaa                                                                                       <  2:1858960/151‑1 (MQ=255)

                                                                   tGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTaaa                                                                                      <  5:730991/150‑1 (MQ=255)

                                                                    gCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGc                                                                                    <  2:365351/151‑1 (MQ=255)

                                                                     cGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCg                                                                                   <  4:580948/151‑1 (MQ=255)

                                                                         gtTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCg                                                                               >  3:987782/1‑151 (MQ=255)

                                                                          ttCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGt                                                                              >  5:1759566/1‑151 (MQ=255)

                                                                            cAGCATCGGCTGGCGTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTaa                                                                            <  4:1233390/151‑1 (MQ=255)

                                                                            cAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTa                                                                             <  1:217208/150‑1 (MQ=255)

                                                                              gCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGt                                                                          <  7:1321771/151‑1 (MQ=255)

                                                                              gCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGt                                                                          <  2:1697391/151‑1 (MQ=255)

                                                                                aTCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCg                                                                        <  1:1767104/151‑1 (MQ=255)

                                                                                   ggCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCg                                                                     <  7:91144/151‑1 (MQ=255)

                                                                                    gCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCgg                                                                    <  7:170267/151‑1 (MQ=255)

                                                                                        gCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATa                                                                 >  8:575842/1‑150 (MQ=255)

                                                                                         ccTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCC‑TTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATc                                                              <  4:1522705/151‑1 (MQ=255)

                                                                                          cTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATc                                                              >  1:1612744/1‑151 (MQ=255)

                                                                                           tttCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTGTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGCTAATcg                                                             >  2:1076473/1‑151 (MQ=255)

                                                                                            ttCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATcgc                                                            >  1:1175551/1‑151 (MQ=255)

                                                                                            ttCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATcg                                                             <  6:260275/150‑1 (MQ=255)

                                                                                            ttCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATc                                                              <  7:668794/149‑1 (MQ=255)

                                                                                             tCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATcgcg                                                           >  1:40926/1‑151 (MQ=255)

                                                                                                   gACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTAc                                                     <  5:318637/151‑1 (MQ=255)

                                                                                                      cTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGAc                                                  >  3:946064/1‑151 (MQ=255)

                                                                                                       ttCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGAcc                                                 >  7:975984/1‑151 (MQ=255)

                                                                                                        tCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGAcc                                                 >  2:1292222/1‑150 (MQ=255)

                                                                                                        tCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGAcc                                                 <  8:1807221/150‑1 (MQ=255)

                                                                                                        tCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGAc                                                  <  1:345107/149‑1 (MQ=255)

                                                                                                         ccGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCAt                                               >  5:1443811/1‑151 (MQ=255)

                                                                                                         ccGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCAt                                               <  4:962060/151‑1 (MQ=255)

                                                                                                          cGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATa                                              <  6:1029195/151‑1 (MQ=255)

                                                                                                           gCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATaa                                             <  2:843635/151‑1 (MQ=255)

                                                                                                            cTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAAt                                            >  7:1530796/1‑151 (MQ=255)

                                                                                                               tGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAAtgtg                                         >  8:1245738/1‑151 (MQ=255)

                                                                                                                  tgtAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTgg                                      >  2:265060/1‑151 (MQ=255)

                                                                                                                    tAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGgcg                                    >  1:679243/1‑151 (MQ=255)

                                                                                                                    tAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGgcg                                    >  6:1564929/1‑151 (MQ=255)

                                                                                                                     aGTTGAGCTTACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTTgcgc                                   >  3:1444167/1‑151 (MQ=255)

                                                                                                                      gTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGcgcg                                  <  8:1634327/151‑1 (MQ=255)

                                                                                                                        tGAGCTTACGCAGCACATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCgg                                 <  2:118464/150‑1 (MQ=255)

                                                                                                                         gAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGt                                <  7:109887/150‑1 (MQ=255)

                                                                                                                           gCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGtttt                             <  3:113641/151‑1 (MQ=255)

                                                                                                                             ttACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGcgcg                                  <  3:1420961/144‑1 (MQ=255)

                                                                                                                             ttACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTAc                           >  3:1080818/1‑151 (MQ=255)

                                                                                                                              tACGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTAcg                          >  6:1895878/1‑150 (MQ=255)

                                                                                                                               aCGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCt                         >  8:858741/1‑151 (MQ=255)

                                                                                                                                cGCAGCAGATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGGTCTACGACCATAATGTGTGGAGCGGTTTTACCTg                        >  1:1014574/1‑151 (MQ=255)

                                                                                                                                cGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTg                        >  5:1114985/1‑151 (MQ=255)

                                                                                                                                   agcagATCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTggg                      >  8:1195916/1‑150 (MQ=255)

                                                                                                                                   agcagACCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGa                     <  6:674677/151‑1 (MQ=255)

                                                                                                                                     cagATCCATCTCGCTATGGGTGTTCCAGGCAATACCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATc                   <  7:251/151‑1 (MQ=255)

                                                                                                                                        aTCCATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGCCCATAATGTGTGGCGCGGt                                >  5:1897423/1‑135 (MQ=255)

                                                                                                                                        aTCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGCCCATAATGTGTGGCGCGGTTTTACCTGGGATCAgg                >  5:470835/1‑151 (MQ=255)

                                                                                                                                        aTCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAgg                >  3:730029/1‑151 (MQ=255)

                                                                                                                                          ccATCTCGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTc              <  4:1424145/151‑1 (MQ=255)

                                                                                                                                          ccATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCGATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGt               <  8:1392507/150‑1 (MQ=255)

                                                                                                                                          ccATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTc              <  5:1689817/151‑1 (MQ=255)

                                                                                                                                               tcGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTACTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCAAAATGTCTGGCGCGGTTTTACCTGGGATCAGGTCGCACt         >  6:43154/1‑151 (MQ=255)

                                                                                                                                                cgcgCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCa           <  4:427372/148‑1 (MQ=255)

                                                                                                                                                cGCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGc            <  4:209260/147‑1 (MQ=255)

                                                                                                                                                cGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCa           <  8:755616/148‑1 (MQ=255)

                                                                                                                                                cGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCAc          <  5:545104/149‑1 (MQ=255)

                                                                                                                                                 gCTCTGGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCg       <  6:1447940/151‑1 (MQ=255)

                                                                                                                                                  ctctAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCGc      >  8:419649/1‑151 (MQ=255)

                                                                                                                                                   tctAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCGcc     >  5:1597307/1‑151 (MQ=255)

                                                                                                                                                     tAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCGCCtt   >  8:1734942/1‑151 (MQ=255)

AGGCTTTCACGGCTACCTGACCGTTGGTTTCCGGTGCCAGTGTCAGGATCATCTCTGCCGCATCAATTGCGGTGTTCAGCATCGGCTGGCCTTTCGCCGGACCTTCCGCTTTGGTGTAGTTGAGCTTACGCAGCAGATCCATCTCGCTCTAGGTGTTCCAGGCAATCCCTTTACCGCCATTACCGATTTTCTCCATCAGCGGGCCGATAGAGGTAAAGCGTTCGTAAGTCGCCGGATAATCGCGCTCTACGACCATAATGTGTGGCGCGGTTTTACCTGGGATCAGGTCGCACTCGCCTTT  >  DH1/2595701‑2596001

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