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%	496.7 / NA	140	Q47L (CAG→CTG)	phnK	ATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (55/85); total (55/85)

Only 100 of 146 total aligned reads displayed.

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGC  >  NC_000913/4316113‑4316290

 ggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCg                                                                                      >  1:396146/1‑93 (MQ=255)

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

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

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

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

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

   tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGCAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                    <  2:2184757/93‑1 (MQ=255)

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

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

     cgCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGt                                                                                  <  2:3589144/92‑1 (MQ=255)

      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCTGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                 >  2:2307665/1‑93 (MQ=255)

      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                 <  1:4675828/93‑1 (MQ=255)

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

      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                 <  2:2304206/93‑1 (MQ=255)

       ccTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCt                                                                                <  1:1693585/93‑1 (MQ=255)

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

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

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

          tgCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                             <  1:4182909/92‑1 (MQ=255)

          cGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                             >  1:4506503/1‑93 (MQ=255)

          cGCTCAGTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCTTGGTCTTcc                                                                             <  2:2519444/93‑1 (MQ=255)

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

           gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                            >  1:4158778/1‑93 (MQ=255)

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

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

             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                          <  1:3278299/93‑1 (MQ=255)

              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                         <  1:75418/93‑1 (MQ=255)

              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                         <  2:3598966/93‑1 (MQ=255)

              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                         >  1:3437120/1‑93 (MQ=255)

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

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

                 tGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGcc                                                                      <  2:990287/93‑1 (MQ=255)

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

                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                   <  1:507920/93‑1 (MQ=255)

                     tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                                  <  1:2374669/93‑1 (MQ=255)

                     tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                                  >  1:2415718/1‑93 (MQ=255)

                     tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                                  <  1:810156/93‑1 (MQ=255)

                       cAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTc                                                                <  1:3089136/93‑1 (MQ=255)

                        aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcc                                                               >  2:2842247/1‑93 (MQ=255)

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

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

                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                           >  1:3618990/1‑93 (MQ=255)

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

                              cGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACa                                                         <  2:2988265/93‑1 (MQ=255)

                               ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAcga                                                        >  1:2303470/1‑91 (MQ=255)

                               ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                        >  2:4689022/1‑93 (MQ=255)

                                gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaag                                                       >  1:20488/1‑92 (MQ=255)

                                gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                       <  2:3202855/93‑1 (MQ=255)

                                gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                       <  2:1722728/93‑1 (MQ=255)

                                 ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                                      <  2:1197683/93‑1 (MQ=255)

                                  tctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGc                                                     <  2:19710/93‑1 (MQ=255)

                                   ctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACCATGcc                                                    >  1:2863742/1‑93 (MQ=255)

                                    tcGTAGTGACTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                   <  2:2808571/93‑1 (MQ=255)

                                    tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                   <  2:4397806/93‑1 (MQ=255)

                                       tAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                                >  1:1696894/1‑93 (MQ=255)

                                       tAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                                >  2:2067297/1‑93 (MQ=255)

                                        aGTGAGTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                               <  2:4454517/93‑1 (MQ=255)

                                                tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                                       <  1:1059583/93‑1 (MQ=255)

                                                tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                                       <  2:781878/93‑1 (MQ=255)

                                                tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                                       >  1:3884165/1‑93 (MQ=255)

                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgg                                     >  2:872186/1‑93 (MQ=255)

                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgg                                     >  2:872185/1‑93 (MQ=255)

                                                   cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTGCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGc                                    <  2:1804861/93‑1 (MQ=255)

                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                   <  2:4155138/93‑1 (MQ=255)

                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGCCTTCAGCAGCGTGGTCTGCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                  <  1:638419/93‑1 (MQ=255)

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

                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTCGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                  <  1:17461/93‑1 (MQ=255)

                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                <  2:1659977/93‑1 (MQ=255)

                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                <  1:1244421/93‑1 (MQ=255)

                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                              <  2:2697132/93‑1 (MQ=255)

                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                              <  2:4247280/93‑1 (MQ=255)

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

                                                            gccgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                           <  1:40804/91‑1 (MQ=255)

                                                             gcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATcga                          >  1:1521361/1‑91 (MQ=255)

                                                             gcgTCAGGCGAGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaa                          >  1:3585198/1‑93 (MQ=255)

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

                                                               gTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                        <  1:3511796/93‑1 (MQ=255)

                                                               gTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                        <  1:2615156/93‑1 (MQ=255)

                                                               gTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                        <  1:3511797/93‑1 (MQ=255)

                                                                tCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                       >  2:794761/1‑93 (MQ=255)

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

                                                                tCAGGCCCGGGGAGATCGACTTCAGCAGCCTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                       <  2:3553312/93‑1 (MQ=255)

                                                                   tgcgcgcGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                    >  1:432020/2‑93 (MQ=255)

                                                                   ggCTCGCGGAGACCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                    <  2:564306/93‑1 (MQ=255)

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

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

                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                 >  1:3540107/1‑93 (MQ=255)

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

                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg               >  1:4349592/1‑93 (MQ=255)

                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg               >  2:4675126/1‑93 (MQ=255)

                                                                         cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc              <  1:3398609/93‑1 (MQ=255)

                                                                           gagaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTa            <  2:3001529/93‑1 (MQ=255)

                                                                             gaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa          <  2:3914979/93‑1 (MQ=255)

                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc        >  1:32970/1‑93 (MQ=255)

                                                                                cGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGcc       <  2:1828962/93‑1 (MQ=255)

                                                                                  gcTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt     <  2:1124778/92‑1 (MQ=255)

                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCTGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt     <  1:254714/93‑1 (MQ=255)

                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt     >  2:4710497/1‑93 (MQ=255)

                                                                                   cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt    <  1:1651317/93‑1 (MQ=255)

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

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGC  >  NC_000913/4316113‑4316290

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