BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,316,203	T→A	Q47L (CAG→CTG)	phnK ←	ATP‑binding component of phosphonate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,316,203	0	T	A	100.0%	869.0 / NA	236	Q47L (CAG→CTG)	phnK	ATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (114/122); total (114/122)

Only 100 of 253 total aligned reads displayed.

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292

cggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCagcagc                                                                                         <  1:1109147/93‑1 (MQ=255)

 ggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCg                                                                                        <  1:297842/93‑1 (MQ=255)

  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCgg                                                                                       >  1:624571/1‑92 (MQ=255)

  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       >  1:5999518/1‑93 (MQ=255)

  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       <  2:1074353/93‑1 (MQ=255)

  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       <  1:1838499/93‑1 (MQ=255)

  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGAACGACTTCAGCAGCgg                                                                                       >  2:296226/1‑92 (MQ=255)

   tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      <  1:4789543/93‑1 (MQ=255)

    cGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTgg                                                                                     <  1:6093629/93‑1 (MQ=255)

    cGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTgg                                                                                     >  2:5494951/1‑93 (MQ=255)

     ggCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGt                                                                                    <  1:3889114/93‑1 (MQ=255)

      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   >  2:4768333/1‑93 (MQ=255)

      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   >  1:3748205/1‑93 (MQ=255)

        cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 <  2:2787091/93‑1 (MQ=255)

         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  2:4551751/1‑93 (MQ=255)

         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                <  1:3329298/93‑1 (MQ=255)

         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                <  1:4887792/93‑1 (MQ=255)

           gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              <  2:500183/93‑1 (MQ=255)

            cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCg                                                                             <  1:1756514/93‑1 (MQ=255)

             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  1:4458656/1‑93 (MQ=255)

               aTTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAg                                                                          >  2:2668049/1‑93 (MQ=255)

               aTTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAg                                                                          >  2:2143499/1‑93 (MQ=255)

                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  2:2219540/1‑93 (MQ=255)

                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  1:2746769/1‑93 (MQ=255)

                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         <  1:2497401/93‑1 (MQ=255)

                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  2:3990206/1‑93 (MQ=255)

                  gCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCg                                                                       >  2:3546853/1‑93 (MQ=255)

                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     >  2:4036522/1‑93 (MQ=255)

                    aTACAGCGAACGGTTCTCGTAGTGAATTTACCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     >  2:977278/1‑93 (MQ=255)

                          cGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcccc                                                               <  1:3407260/93‑1 (MQ=255)

                          cGAACGGTTCTCGTAGTGAATGTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcccc                                                               <  2:4233885/93‑1 (MQ=255)

                           gAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCg                                                              >  2:3714881/1‑93 (MQ=255)

                           gAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCg                                                              <  2:4716897/93‑1 (MQ=255)

                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             <  2:2388880/93‑1 (MQ=255)

                             aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCTGGAGCCGGATTCCCCGaa                                                            >  2:4191705/1‑92 (MQ=255)

                             aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                                            <  1:3508199/93‑1 (MQ=255)

                             aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                                            >  2:1921463/1‑93 (MQ=255)

                               ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  1:743694/93‑1 (MQ=255)

                               ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  2:1858369/93‑1 (MQ=255)

                                 ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                                        >  1:5318435/1‑93 (MQ=255)

                                 ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGACTTCCCGGAGCCGGATTCCCCGACACTg                                                        >  1:1497208/1‑93 (MQ=255)

                                  tctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGc                                                       >  1:315423/1‑93 (MQ=255)

                                   ctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGcc                                                      <  2:5192681/93‑1 (MQ=255)

                                    tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  1:469900/1‑93 (MQ=255)

                                     cGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCa                                                    <  2:1031504/93‑1 (MQ=255)

                                     cGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCa                                                    <  2:3621311/93‑1 (MQ=255)

                                      gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   >  2:5349698/1‑93 (MQ=255)

                                      gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   <  1:2465283/93‑1 (MQ=255)

                                        aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 >  2:6134203/1‑93 (MQ=255)

                                        aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 <  1:809216/93‑1 (MQ=255)

                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  1:5050203/93‑1 (MQ=255)

                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  1:4014925/93‑1 (MQ=255)

                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             >  2:759437/1‑93 (MQ=255)

                                              tttCCCCCTGCTGCGGCGTCAGGCGCGGGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccc                                           <  2:3291567/93‑1 (MQ=255)

                                               ttCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcccc                                          <  2:5505006/93‑1 (MQ=255)

                                                 cccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCg                                        >  1:5034626/1‑93 (MQ=255)

                                                 cccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCg                                        <  1:1580291/93‑1 (MQ=255)

                                                   cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTGCCCCGGc                                      >  1:1628706/1‑93 (MQ=255)

                                                   cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGc                                      >  2:5953100/1‑93 (MQ=255)

                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     >  2:5115247/1‑93 (MQ=255)

                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     >  1:4644027/1‑93 (MQ=255)

                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:5260857/1‑93 (MQ=255)

                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:2296304/1‑93 (MQ=255)

                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    <  2:5588297/93‑1 (MQ=255)

                                                      tgctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCAt                                   >  1:6232285/1‑93 (MQ=255)

                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGTTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                  <  2:1536965/93‑1 (MQ=255)

                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                  >  1:5480668/1‑93 (MQ=255)

                                                        ctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaa                                 <  2:5820476/93‑1 (MQ=255)

                                                        ctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaa                                 >  1:6120615/1‑93 (MQ=255)

                                                          gcggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAAt                               >  1:1965974/1‑93 (MQ=255)

                                                           cggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATc                              <  2:655826/93‑1 (MQ=255)

                                                             gcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaa                            >  1:4556933/1‑93 (MQ=255)

                                                             gcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaa                            <  1:5832045/93‑1 (MQ=255)

                                                             gcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaa                            >  1:2166950/1‑93 (MQ=255)

                                                              cgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaa                           <  2:4911954/93‑1 (MQ=255)

                                                               gTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                          >  2:4655024/1‑93 (MQ=255)

                                                                tCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         <  1:5332427/93‑1 (MQ=255)

                                                                 cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        >  2:3156153/1‑93 (MQ=255)

                                                                 cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        <  2:2549542/93‑1 (MQ=255)

                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      >  1:1420144/1‑93 (MQ=255)

                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:5599117/93‑1 (MQ=255)

                                                                    gcgcgcGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGAc                     <  2:4645889/93‑1 (MQ=255)

                                                                     cgcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  2:5325469/93‑1 (MQ=255)

                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  2:5681040/93‑1 (MQ=255)

                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  2:5850378/93‑1 (MQ=255)

                                                                       cgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATc                  <  1:1188201/93‑1 (MQ=255)

                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGGGCCGGATTCCCCGACAATGGCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg                 >  2:3869476/1‑93 (MQ=255)

                                                                          ggAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCt               >  2:3419822/1‑93 (MQ=255)

                                                                           gagaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTa              >  2:3061183/1‑93 (MQ=255)

                                                                             gaTCTACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTGCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            >  2:3483068/1‑93 (MQ=255)

                                                                             gaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            >  1:2058768/1‑93 (MQ=255)

                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          <  2:5034818/93‑1 (MQ=255)

                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          <  1:4162660/93‑1 (MQ=255)

                                                                                cGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGcc         >  2:3818155/1‑93 (MQ=255)

                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:3720963/93‑1 (MQ=255)

                                                                                    ttCAGCAGCGTGGTCTTCGCGGAGCCTGATACCCCGACAATGCCCAGCTCTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     <  2:1941050/93‑1 (MQ=255)

                                                                                    ttCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     >  2:4463839/1‑93 (MQ=255)

                                                                                     tCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGc    <  2:1312792/93‑1 (MQ=255)

                                                                                      cagcagCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGcc   <  1:3908007/93‑1 (MQ=255)

                                                                                       agcagcGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  <  2:2300998/93‑1 (MQ=255)

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292

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