BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,566,875	A→G	R29R (CGT→CGC)	glgC ←	glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,875	0	A	G	100.0%	762.1 / NA	246	R29R (CGT→CGC)	glgC	glucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (162/84); total (162/84)

Only 100 of 338 total aligned reads displayed.

AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964

 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGt                                                                                        <  2:5207269/93‑1 (MQ=255)

  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtc                                                                                       <  1:4946850/93‑1 (MQ=255)

  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtc                                                                                       >  2:3901384/1‑93 (MQ=255)

   cAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcc                                                                                      <  1:5210927/93‑1 (MQ=255)

   cAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcc                                                                                      <  2:2345897/93‑1 (MQ=255)

    aGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcct                                                                                     <  1:3853034/93‑1 (MQ=255)

     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  2:602568/1‑93 (MQ=255)

     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  2:2562040/1‑93 (MQ=255)

      cgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctcc                                                                                   >  2:7478004/1‑93 (MQ=255)

       gcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTccc                                                                                  >  1:1358101/1‑93 (MQ=255)

        cAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCg                                                                                 <  1:6464686/93‑1 (MQ=255)

          aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               >  1:6806337/1‑93 (MQ=255)

           aGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCa                                                                              <  2:1940105/93‑1 (MQ=255)

           aGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCa                                                                              <  1:7703893/93‑1 (MQ=255)

           aGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCa                                                                              <  2:1428607/93‑1 (MQ=255)

            gTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCACGCGGCCTCCCGCCac                                                                             >  1:233037/1‑92 (MQ=255)

             tCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGt                                                                            >  1:613228/1‑93 (MQ=255)

              cGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           <  1:1130665/93‑1 (MQ=255)

                aTAATGCGGAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATc                                                                         <  2:5355534/93‑1 (MQ=255)

                 tAATGCGGAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCa                                                                        <  2:5141873/93‑1 (MQ=255)

                  aaTGCGGAACTTCCCGCGGAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAg                                                                       <  1:2586576/93‑1 (MQ=255)

                   aTGCGGACCTTCCCGCCGAGGTGTACGGCCGTTTTTGCTCGTTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      <  1:4371395/93‑1 (MQ=255)

                    tGCGTAACTTCCCGCCGAGTTGTACGGCCGGTTTGCCTCGTTTATTGGTTAAATCCTTCGGGCGGGTCCCGCGTCCTCCCGCCAGTATCAggg                                                                     <  2:6235615/93‑1 (MQ=255)

                       ggAACTTACCGCCGAGGTGTCCGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCaa                                                                  <  2:2603983/93‑1 (MQ=255)

                        gAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGTTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 <  2:5716984/93‑1 (MQ=255)

                        gAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 >  2:3197360/1‑93 (MQ=255)

                         aaCTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACa                                                                <  1:7901430/93‑1 (MQ=255)

                          aCTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               <  1:1660388/93‑1 (MQ=255)

                          aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               >  2:6634084/1‑93 (MQ=255)

                          aCTTACCGCCGAAGTGTACGCCCGGTTTTGCTCGCTTATTGGTTAATTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               <  2:1126688/93‑1 (MQ=255)

                            ttACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAt                                                             >  2:673160/1‑93 (MQ=255)

                             tACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAtt                                                            >  1:6285974/1‑93 (MQ=255)

                                 gccgCAGTGTACGGCCGGGTTTGCTCGCTTATTGGTTAAATACTTCAGGCGGGTACCGCGTCCGCCCGCCAGTATCAGGGCAACAGAGTTGaa                                                        >  2:7105023/1‑93 (MQ=255)

                                 gccgAAGTGTAGGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        <  1:2861600/93‑1 (MQ=255)

                                 gccgAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTGCGGGGGGGTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        <  2:1503228/93‑1 (MQ=255)

                                 gccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTAAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        <  2:312434/93‑1 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTTGTTAAATCCTTCAGGCGGGAACCGCGTCCTCCCGCCAGTATCAGGGCAAAAGGTTTCAAt                                                       >  2:5871732/1‑93 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  2:166485/93‑1 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  1:5020613/93‑1 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGAGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       >  1:4893728/1‑93 (MQ=255)

                                   cgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATg                                                      <  1:3623205/93‑1 (MQ=255)

                                    gAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATgg                                                     >  2:2083222/1‑93 (MQ=255)

                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    <  1:6731396/93‑1 (MQ=255)

                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    >  1:869823/1‑93 (MQ=255)

                                       gtgtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAg                                                  >  2:4204910/1‑93 (MQ=255)

                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                >  2:7373397/1‑93 (MQ=255)

                                           gcGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTgg                                              <  1:4631431/92‑1 (MQ=255)

                                           aCGGCGGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTgg                                              <  1:6552619/93‑1 (MQ=255)

                                             ggCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcg                                            >  2:6053905/1‑93 (MQ=255)

                                              gCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcgc                                           >  1:5800201/1‑93 (MQ=255)

                                                cGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcgcgc                                         >  2:5569153/1‑93 (MQ=255)

                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        >  2:7258865/1‑93 (MQ=255)

                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGGACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAAGGCAGGTGgcgcggc                                        >  2:2287554/1‑91 (MQ=255)

                                                  gTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCa                                       >  2:4925331/1‑93 (MQ=255)

                                                   ttttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCaa                                      >  2:239353/1‑93 (MQ=255)

                                                   ttttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCaa                                      >  2:1801835/1‑93 (MQ=255)

                                                    tttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAAc                                     >  1:7064028/1‑93 (MQ=255)

                                                     ttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGGACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACa                                    <  1:7651652/93‑1 (MQ=255)

                                                       gCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCCACAGATTTCAATGGCAGCTGGCGCGCCAACAtt                                  >  1:3916259/1‑93 (MQ=255)

                                                         tctcTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  1:1988070/90‑1 (MQ=255)

                                                         ccGCTTATTGTTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  1:5299940/92‑1 (MQ=255)

                                                           ccTTATTGGTTAAATCCTCCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              <  2:7775288/92‑1 (MQ=255)

                                                             ttATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGa                            >  1:4214873/1‑93 (MQ=255)

                                                             ttATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGa                            >  2:3623097/1‑93 (MQ=255)

                                                             ttATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGa                            <  2:3827069/93‑1 (MQ=255)

                                                              tATTGTTTAACTCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           <  2:6798556/93‑1 (MQ=255)

                                                              tATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           >  2:5647593/1‑93 (MQ=255)

                                                              tATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           >  2:3817417/1‑93 (MQ=255)

                                                               tttGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          <  1:2460205/92‑1 (MQ=255)

                                                               aTTGGTTAAACCCTTCAGGCGGGCACCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          <  1:1476266/93‑1 (MQ=255)

                                                                 tGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGt                        <  2:2464983/93‑1 (MQ=255)

                                                                  ggTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGtt                       <  2:4091027/93‑1 (MQ=255)

                                                                   gTTAATTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      <  1:837316/93‑1 (MQ=255)

                                                                   gTTAAATCCTTCGGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      <  1:4771456/93‑1 (MQ=255)

                                                                   gTTAAATCCTTCAGGCGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      <  1:2235810/93‑1 (MQ=255)

                                                                    ttAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     >  1:6667870/1‑93 (MQ=255)

                                                                      aaaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   >  2:7666594/1‑93 (MQ=255)

                                                                      aaaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   <  1:4116705/93‑1 (MQ=255)

                                                                        aTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTtctc                 <  1:7438894/93‑1 (MQ=255)

                                                                        aTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGCGATCGTTCTtctc                 >  1:4291155/1‑93 (MQ=255)

                                                                         tCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTctct                >  1:2081639/1‑93 (MQ=255)

                                                                           cTTCGGGCGGGTCCGCCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaa              <  1:1431857/93‑1 (MQ=255)

                                                                           cTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaa              >  2:5819989/1‑93 (MQ=255)

                                                                           cTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaa              >  1:6878521/1‑93 (MQ=255)

                                                                            ttCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaaa             >  1:3458450/1‑93 (MQ=255)

                                                                             tcgggCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAAc            <  2:458580/90‑1 (MQ=255)

                                                                               aGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          >  2:479319/1‑93 (MQ=255)

                                                                               aGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          >  1:5631443/1‑93 (MQ=255)

                                                                                ggCGGGTCCCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         <  1:5906687/93‑1 (MQ=255)

                                                                                ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  1:3045455/1‑93 (MQ=255)

                                                                                ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  1:3342590/1‑93 (MQ=255)

                                                                                 gggggTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAc        <  2:2853983/91‑1 (MQ=255)

                                                                                 gCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAc        <  1:5406093/93‑1 (MQ=255)

                                                                                  ggggTACCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAcc       <  2:7915264/92‑1 (MQ=255)

                                                                                   gggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      >  1:5149464/1‑93 (MQ=255)

                                                                                   gggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      >  2:4185416/1‑93 (MQ=255)

                                                                                    gttACCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     <  1:4673523/91‑1 (MQ=255)

                                                                                    ggTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     <  1:7246269/93‑1 (MQ=255)

                                                                                    ggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     <  1:2686348/93‑1 (MQ=255)

                                                                                       aCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAc  >  1:4468840/1‑93 (MQ=255)

AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964

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