BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	2,688,095	0	A	G	52.9%	76.6 / inf	438	K115E (AAA→GAA)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand): ref base A (88/118); new base G (102/129); total (190/248)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

Only 100 of 445 total aligned reads displayed.

GGTTACCTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATCGCGATATC  >  DH1/2687945‑2688244

ggTTACCTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTa                                                                                                                                                       <  5:56324/151‑1 (MQ=255)

   tACCTTGCTGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGTAACTGCACGCGGTTTCTgaa                                                                                                                                                     <  6:122970/150‑1 (MQ=255)

   tACCTTGATGCATTCGTCGCCAGTTTCAGGGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTgaa                                                                                                                                                     <  5:986255/150‑1 (MQ=255)

   tACCTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAAt                                                                                                                                                    >  1:704398/1‑151 (MQ=255)

    aCCTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAAtt                                                                                                                                                   <  1:1800927/151‑1 (MQ=255)

     ccTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTc                                                                                                                                                  >  3:520879/1‑151 (MQ=255)

     ccTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTc                                                                                                                                                  <  2:344869/151‑1 (MQ=255)

     ccTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGGTTCTGAATTc                                                                                                                                                  >  8:1639303/1‑151 (MQ=255)

        tGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATc                                                                                                                                               >  6:505276/1‑151 (MQ=255)

         gATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATcc                                                                                                                                              <  3:1465136/151‑1 (MQ=255)

            gcattgctcgCCACGTTCAGCGAAGAACTCGGGGCGCGCTTGAAAACTATGCTGGGAACCAGTAACCACCACCTTAAAGCGGTTCAGTATTTCCTGAAAGATAATGTGGCGGAGATCCCGGAACTGCAGGCGGGTTGTGAAATCATCCag                                                                                                                                            <  2:1131680/143‑2 (MQ=255)

             cAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAGGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACtt                                                                                                                                          <  8:408782/151‑1 (MQ=255)

                 cgtcgCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGc                                                                                                                                       <  5:134854/150‑1 (MQ=255)

                 cgtcgCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGcc                                                                                                                                      >  7:995981/1‑151 (MQ=255)

                  gtcgCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCt                                                                                                                                     >  7:163930/1‑151 (MQ=255)

                    cgCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGt                                                                                                                                   <  4:272210/151‑1 (MQ=255)

                    cgCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGt                                                                                                                                   >  6:1497570/1‑151 (MQ=255)

                        aGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACtt                                                                                                                               >  1:1542846/1‑151 (MQ=255)

                         gTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTc                                                                                                                              <  6:1690109/151‑1 (MQ=255)

                            tCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGGCTGTACTTCgg                                                                                                                            >  8:1314936/1‑150 (MQ=255)

                              aGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGa                                                                                                                           >  3:299700/1‑149 (MQ=255)

                                   agaagaTGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATc                                                                                                                    <  5:385001/151‑1 (MQ=255)

                                    gaagaTGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCa                                                                                                                   <  2:669882/151‑1 (MQ=255)

                                       gaTGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATa                                                                                                                <  7:547969/151‑1 (MQ=255)

                                        aTGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATaa                                                                                                               <  6:1372785/151‑1 (MQ=255)

                                         tGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCGTGTACTTCGGAAGATATCAATaa                                                                                                               <  4:412141/150‑1 (MQ=255)

                                          gcttcgctaaTCAAAACTATCGAGCGTACCAATAACCACGAAGTTAAAGCGGTTGAGTATTTCCAGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAAcc                                                                                                             <  8:582696/142‑1 (MQ=255)

                                          gcggcgCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAAcc                                                                                                             <  4:191772/151‑1 (MQ=255)

                                             gcgcgcATCAAAACTATCGAGCGTAACACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACctct                                                                                                          <  6:418548/151‑1 (MQ=255)

                                                cgcATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTccc                                                                                                       >  3:1140007/1‑151 (MQ=255)

                                                 gcATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCa                                                                                                      <  8:302210/151‑1 (MQ=255)

                                                  cATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCAc                                                                                                     <  3:449152/151‑1 (MQ=255)

                                                  cATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCAc                                                                                                     <  7:1017032/151‑1 (MQ=255)

                                                   aTCAAAACTATCGCGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACg                                                                                                    <  2:774376/151‑1 (MQ=255)

                                                   aTCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGGAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACg                                                                                                    <  4:1732435/151‑1 (MQ=255)

                                                   aTCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACg                                                                                                    <  1:342849/151‑1 (MQ=255)

                                                     cAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGc                                                                                                   >  3:826872/1‑150 (MQ=255)

                                                      aaaaCTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCAt                                                                                                 <  5:442347/151‑1 (MQ=255)

                                                        aaCTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTa                                                                                               <  8:1064821/151‑1 (MQ=255)

                                                           tATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATg                                                                                            >  3:1665445/1‑151 (MQ=255)

                                                            aTCCAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACTCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATg                                                                                            <  8:524111/150‑1 (MQ=255)

                                                              cGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTACACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTg                                                                                         <  1:791247/151‑1 (MQ=255)

                                                                 gCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGaa                                                                                       >  3:1876214/1‑150 (MQ=255)

                                                                  cGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGaaaa                                                                                     <  1:258179/151‑1 (MQ=255)

                                                                  cGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGaaa                                                                                      >  4:907615/1‑150 (MQ=255)

                                                                  cGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGGGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGaaaa                                                                                     <  7:1220104/151‑1 (MQ=255)

                                                                   gTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAAc                                                                                    <  6:1792100/151‑1 (MQ=255)

                                                                   gTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAAc                                                                                    <  1:1419849/151‑1 (MQ=255)

                                                                    tACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAccc                                                                                   >  8:178739/1‑151 (MQ=255)

                                                                     aCCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGaaaaac                                                                                   <  1:1059340/150‑3 (MQ=255)

                                                                     aCCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAAcc                                                                                   <  7:1922490/150‑1 (MQ=255)

                                                                     aCCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACcg                                                                                  <  5:1439988/151‑1 (MQ=255)

                                                                      ccACTAACCACGACGTTAAAGCGGTTGAGTAATTCCTGAAAGAAAAACTGCCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACcgc                                                                                 <  5:1537457/151‑1 (MQ=255)

                                                                      ccAATAACCACGACGTTAAATCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACcgc                                                                                 <  7:664967/151‑1 (MQ=255)

                                                                         atAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACcgcgcg                                                                              <  4:1028931/150‑1 (MQ=255)

                                                                           aaCCACGACGGAAAAGCCGGTGACCACTTCCCGAAACAAAAAGCGCCGGACCTCCCGGCACCGCACGCGGTTTCAAAATTCATCCACATTGCCTGTACTTCCGAACAAATCAAAAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGtg                                                                            <  6:1864800/151‑1 (MQ=255)

                                                                             ccACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGtgat                                                                          >  4:1431906/1‑151 (MQ=255)

                                                                              cACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGtgatg                                                                         <  2:542645/151‑1 (MQ=255)

                                                                               acgacgTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGtgatga                                                                        <  1:1353665/151‑1 (MQ=255)

                                                                                cgacgTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGaa                                                                       <  7:1512891/151‑1 (MQ=255)

                                                                                 gacgTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGtgatga                                                                        >  6:1066924/1‑149 (MQ=255)

                                                                                   cgTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAAATCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTg                                                                    <  5:1090209/151‑1 (MQ=255)

                                                                                   cgTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTg                                                                    <  2:1500044/151‑1 (MQ=255)

                                                                                     ttAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGAt                                                                  >  1:1140815/1‑151 (MQ=255)

                                                                                      tAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATc                                                                 >  5:793861/1‑151 (MQ=255)

                                                                                       aaaGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATc                                                                 <  4:540718/150‑1 (MQ=255)

                                                                                         aGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTg                                                              >  8:83652/1‑151 (MQ=255)

                                                                                          gCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGc                                                             >  3:943595/1‑151 (MQ=255)

                                                                                            ggTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGcc                                                            >  8:779520/1‑150 (MQ=255)

                                                                                                gAGTATTTCCTGAAAGAAAAAGTGGCGGAGATGCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACt                                                       <  5:92815/151‑1 (MQ=255)

                                                                                                gAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACt                                                       <  1:438093/151‑1 (MQ=255)

                                                                                                  gTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTgg                                                     <  1:1902530/151‑1 (MQ=255)

                                                                                                     tttCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGt                                                  >  7:771792/1‑151 (MQ=255)

                                                                                                       tCCTGAAAGAAAAAGTGGCGGAGACCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGTGTGATGAAGTGATCCTGCCATACTGGCGTc                                                 <  3:1800595/150‑1 (MQ=255)

                                                                                                         cTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAAc                                              >  8:170041/1‑151 (MQ=255)

                                                                                                              agaaaAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCCCTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAACCCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCACCtgatt                                         >  8:1588948/1‑151 (MQ=255)

                                                                                                                aaaaaGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGGTGAAGTGATCCTGCCATACTGGCGTCAACtgattga                                       <  1:12249/151‑1 (MQ=255)

                                                                                                                aaaaaGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCGATACTGGCGTCAACtgattg                                        >  3:779570/1‑150 (MQ=255)

                                                                                                                    aGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGc                                   <  4:1514881/151‑1 (MQ=255)

                                                                                                                       ggCGGAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCAt                                 <  5:1799112/150‑1 (MQ=255)

                                                                                                                        gCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTa                               <  6:241896/151‑1 (MQ=255)

                                                                                                                        gCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTa                               <  4:1217112/151‑1 (MQ=255)

                                                                                                                          ggAGATCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTaa                              >  6:129522/1‑150 (MQ=255)

                                                                                                                          ggAGAGCCCGGAACTGCACGCGGTTTCTGAATTCCTCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCGGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTaaa                             <  7:138659/151‑1 (MQ=255)

                                                                                                                             gaTCCCGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTaa                              >  1:164228/1‑147 (MQ=255)

                                                                                                                              aTCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTaa                              >  2:1443761/1‑146 (MQ=255)

                                                                                                                                cccGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGAtctc                       <  7:426731/151‑1 (MQ=255)

                                                                                                                                cccGGAACTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGAtctc                       >  5:1605508/1‑151 (MQ=255)

                                                                                                                                 ccGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCg                      >  7:889500/1‑151 (MQ=255)

                                                                                                                                   ggAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGcc                    >  4:1252929/1‑151 (MQ=255)

                                                                                                                                     aaCTGCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGt                  >  7:243448/1‑151 (MQ=255)

                                                                                                                                     aaCTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGt                  >  7:1254731/1‑151 (MQ=255)

                                                                                                                                         gCACTCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCa               <  2:984581/150‑1 (MQ=255)

                                                                                                                                         gCACGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCa               <  6:732805/150‑1 (MQ=255)

                                                                                                                                           aCGCGGTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGt             >  8:1444529/1‑150 (MQ=255)

                                                                                                                                               ggTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATcgc        >  2:1796034/1‑151 (MQ=255)

                                                                                                                                                gTTTCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGCCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATcgcc       >  2:1282811/1‑150 (MQ=255)

                                                                                                                                                 tttCTAAATTCATCCACTTAGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATcgcg       >  7:1855391/1‑150 (MQ=255)

                                                                                                                                                   tCTGAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAACCCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATCGCGat     >  8:257910/1‑150 (MQ=255)

                                                                                                                                                   tCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATCGCGat     >  6:951303/1‑150 (MQ=255)

GGTTACCTTGATGCAATCGTCGCCAGTTTCAGCGAAGAAGATGCGGCGCGCATCAAAACTATCGAGCGTACCACTAACCACGACGTTAAAGCGGTTGAGTATTTCCTGAAAGAAAAAGTGGCGGAGATCCCGGAACTGCACGCGGTTTCTAAATTCATCCACTTTGCCTGTACTTCGGAAGATATCAATAACCTCTCCCACGCATTAATGCTGAAAACCGCGCGTGATGAAGTGATCCTGCCATACTGGCGTCAACTGATTGATGGCATTAAAGATCTCGCCGTTCAGTATCGCGATATC  >  DH1/2687945‑2688244

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