BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,225,547	G→A	G972D (GGC→GAC)	ebgA →	evolved beta‑D‑galactosidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,225,547	0	G	A	100.0%	1220.4 / NA	299	G972D (GGC→GAC)	ebgA	evolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (176/123); total (176/123)

Only 100 of 299 total aligned reads displayed.

CTATACCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGGCTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGGAGAA  >  NC_000913/3225449‑3225644

 tataCCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTc                                                                                                >  2:5095296/1‑101 (MQ=255)

  ataCCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTcc                                                                                               >  1:5095297/1‑101 (MQ=255)

   taCCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCa                                                                                              >  2:5095302/1‑101 (MQ=255)

   taCCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCa                                                                                              >  1:5095303/1‑101 (MQ=255)

     cccAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAAc                                                                                            <  1:5094811/101‑1 (MQ=255)

      ccAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACt                                                                                           >  1:5095306/1‑101 (MQ=255)

       cAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTc                                                                                          >  2:5095310/1‑101 (MQ=255)

        aGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTcc                                                                                         >  2:5095314/1‑101 (MQ=255)

        aGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTcc                                                                                         >  1:5095315/1‑101 (MQ=255)

        aGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTcc                                                                                         <  1:5094933/101‑1 (MQ=255)

             aaCATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTgggg                                                                                    >  2:5095320/1‑101 (MQ=255)

             aaCATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTgggg                                                                                    >  1:5095321/1‑101 (MQ=255)

              aCATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGc                                                                                   >  2:5095326/1‑101 (MQ=255)

               cATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCa                                                                                  >  1:5095328/1‑101 (MQ=255)

                 tCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGc                                                                                >  1:5095332/1‑101 (MQ=255)

                 tCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGc                                                                                <  1:5094457/101‑1 (MQ=255)

                   cACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGa                                                                              >  2:5095334/1‑101 (MQ=255)

                    aCGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAg                                                                             <  1:5094922/101‑1 (MQ=255)

                     cGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAgg                                                                            <  1:5094825/101‑1 (MQ=255)

                     cGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAgg                                                                            <  1:5094901/101‑1 (MQ=255)

                      gCTGCCCAGCGCTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGt                                                                           <  1:5094833/101‑1 (MQ=255)

                      gCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGt                                                                           >  1:5095338/1‑101 (MQ=255)

                       cTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTg                                                                          <  1:5095018/101‑1 (MQ=255)

                       cTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTg                                                                          >  2:5095340/1‑101 (MQ=255)

                        tGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGACTCGACTCCAACTCCTGGGGCAGCGAGGTGc                                                                         >  2:5095344/1‑101 (MQ=255)

                         gCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCt                                                                        >  1:5095350/1‑101 (MQ=255)

                         gCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACCCCAACTCCTGGGGCAGCGAGGTGCt                                                                        >  2:5095347/1‑101 (MQ=255)

                          cccAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTg                                                                       >  2:5095351/1‑101 (MQ=255)

                              gCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACt                                                                   <  2:5094636/101‑1 (MQ=255)

                                aCTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTcc                                                                 <  1:5095110/101‑1 (MQ=255)

                                 cTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCt                                                                >  2:5095357/1‑101 (MQ=255)

                                 cTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCt                                                                >  1:5095356/1‑101 (MQ=255)

                                  tGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTg                                                               <  1:5094873/101‑1 (MQ=255)

                                  tGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTg                                                               >  2:5095361/1‑101 (MQ=255)

                                   gTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTgg                                                              <  1:5094883/101‑1 (MQ=255)

                                    tAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGgc                                                             >  1:5095365/1‑101 (MQ=255)

                                     aaCGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGgcg                                                            >  2:5095367/1‑101 (MQ=255)

                                      aCGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGgcgc                                                           >  1:5095369/1‑101 (MQ=255)

                                         aGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTg                                                                          >  2:5095375/1‑83 (MQ=255)

                                         aGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTc                                                        >  1:5095374/1‑101 (MQ=255)

                                           cTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTg                                                      >  2:5095378/1‑101 (MQ=255)

                                           cTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTg                                                      <  2:5094898/101‑1 (MQ=255)

                                            tGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTgg                                                     <  1:5095074/101‑1 (MQ=255)

                                            tGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTgg                                                     >  2:5095381/1‑101 (MQ=255)

                                              cAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGtt                                                   >  1:5095386/1‑101 (MQ=255)

                                              cAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGtt                                                   >  1:5095389/1‑101 (MQ=255)

                                                gcgcAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTcc                                                 >  2:5095391/1‑101 (MQ=255)

                                                gcgcAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTcc                                                 <  1:5094758/101‑1 (MQ=255)

                                                  gcAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGt                                               >  1:5095394/1‑101 (MQ=255)

                                                  gcAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGt                                               <  2:5094860/101‑1 (MQ=255)

                                                   cAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTg                                              <  2:5094937/101‑1 (MQ=255)

                                                   cAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTg                                              <  1:5095084/101‑1 (MQ=255)

                                                    aGTGATGACATCCCCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTg                                                      <  2:5095035/92‑1 (MQ=255)

                                                    aGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCt                                                                >  1:5095398/1‑82 (MQ=255)

                                                     gTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGAc                                            >  2:5095399/1‑101 (MQ=255)

                                                      tgatgaCATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACt                                           <  1:5094682/101‑1 (MQ=255)

                                                        atgaCATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTc                                         >  1:5095404/1‑101 (MQ=255)

                                                        atgaCATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTc                                         <  2:5094969/101‑1 (MQ=255)

                                                        atgaCATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTc                                         <  2:5094836/101‑1 (MQ=255)

                                                          gaCATCACCCTGAACCTTGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAg                                       <  1:5095242/101‑1 (MQ=255)

                                                          gaCATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAg                                       <  2:5095244/101‑1 (MQ=255)

                                                            cATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCt                                     >  1:5095407/1‑101 (MQ=255)

                                                              tCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCt                                                                                        <  1:5094780/48‑1 (MQ=255)

                                                              tCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTAc                                   >  2:5095409/1‑101 (MQ=255)

                                                                aCCCTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACgg                                 >  1:5095412/1‑101 (MQ=255)

                                                                 cccTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGc                                <  2:5094988/101‑1 (MQ=255)

                                                                  ccTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCt                               <  1:5094949/101‑1 (MQ=255)

                                                                   cTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCtt                              >  1:5095416/1‑101 (MQ=255)

                                                                   cTGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCtt                              <  1:5094850/101‑1 (MQ=255)

                                                                    tGAACCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCttt                             >  2:5095418/1‑101 (MQ=255)

                                                                       aCCTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACg                          >  2:5095422/1‑101 (MQ=255)

                                                                        ccTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGt                         >  1:5095423/1‑101 (MQ=255)

                                                                         cTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGtt                        >  1:5095429/1‑101 (MQ=255)

                                                                         cTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGtt                        <  2:5094874/101‑1 (MQ=255)

                                                                         cTCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGtt                        >  2:5095428/1‑101 (MQ=255)

                                                                          tCGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCt                                                       >  1:5095434/1‑69 (MQ=255)

                                                                           cGATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgc                      >  2:5095436/1‑101 (MQ=255)

                                                                            gATCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgct                     <  2:5095085/101‑1 (MQ=255)

                                                                             aTCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgctg                    <  1:5094976/101‑1 (MQ=255)

                                                                             aTCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgctg                    >  2:5095440/1‑101 (MQ=255)

                                                                              tCACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgctgc                   >  1:5095441/1‑101 (MQ=255)

                                                                               cACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTtgct                     >  2:5095447/1‑98 (MQ=255)

                                                                               cACCAGCTGCTTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGcc                  >  1:5095445/1‑101 (MQ=255)

                                                                                     ctgctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTc            >  2:5095452/1‑101 (MQ=255)

                                                                                     ctgctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTc            <  2:5095179/101‑1 (MQ=255)

                                                                                      tgctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCt           <  1:5095043/101‑1 (MQ=255)

                                                                                      tgctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCt           >  1:5095453/1‑101 (MQ=255)

                                                                                       gctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTg          <  1:5094964/101‑1 (MQ=255)

                                                                                        ctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTgg         <  1:5094617/101‑1 (MQ=255)

                                                                                        ctTGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTgg         <  2:5094788/101‑1 (MQ=255)

                                                                                         ttGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCg                 >  2:5095457/1‑92 (MQ=255)

                                                                                          tGGCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCg       <  1:5095165/101‑1 (MQ=255)

                                                                                            gCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGga     <  2:5094890/101‑1 (MQ=255)

                                                                                            gCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGga     >  1:5095461/1‑101 (MQ=255)

                                                                                            gCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGga     <  1:5095181/101‑1 (MQ=255)

                                                                                            gCCTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGga     >  2:5095462/1‑101 (MQ=255)

                                                                                             ccTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGgag    <  2:5095016/101‑1 (MQ=255)

                                                                                             ccTCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGgag    <  2:5094954/101‑1 (MQ=255)

                                                                                               tCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGGAGaa  >  1:5095467/1‑101 (MQ=255)

                                                                                               tCGACTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGGAGaa  <  1:5095141/101‑1 (MQ=255)

CTATACCCAGGAAAACATCCACGCTGCCCAGCACTGTAACGAGCTGCAGCGCAGTGATGACATCACCCTGAACCTCGATCACCAGCTGCTTGGCCTCGGCTCCAACTCCTGGGGCAGCGAGGTGCTGGACTCCTGGCGCGTCTGGTTCCGTGACTTCAGCTACGGCTTTACGTTGCTGCCGGTTTCTGGCGGAGAA  >  NC_000913/3225449‑3225644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`