BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,586,558	A→G	Q315R (CAG→CGG)	acrD →	AcrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,586,558	0	A	G	98.5%	429.0 / ‑2.4	136	Q315R (CAG→CGG)	acrD	AcrD
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (62/71); total (64/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 142 total aligned reads displayed.

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCT  >  NC_000913/2586468‑2586648

cagccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGt                                                                                          <  2:3364726/93‑1 (MQ=255)

cagccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGt                                                                                          <  1:2563658/93‑1 (MQ=255)

 cgccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                         >  2:1501875/2‑93 (MQ=255)

 agccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                         >  2:381545/1‑93 (MQ=255)

 agccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                         >  1:695084/1‑93 (MQ=255)

 agccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                         <  1:252548/93‑1 (MQ=255)

 agccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                         >  1:2412256/1‑93 (MQ=255)

  gccTCCGGGCTGGGCGTAAAACTGGCCGCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                        <  1:3346946/93‑1 (MQ=255)

   ccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAtt                                                                                       <  1:1089525/93‑1 (MQ=255)

    cTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAttt                                                                                      >  1:3840094/1‑93 (MQ=255)

      ccGGGCTGGGGGTAAAACTGGCCTCCGGCTATAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTcc                                                                                    <  2:3759495/93‑1 (MQ=255)

      ccGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTcc                                                                                    >  2:2645358/1‑93 (MQ=255)

       cGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTccc                                                                                   <  1:230434/93‑1 (MQ=255)

        gggCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCg                                                                                  >  1:443305/1‑93 (MQ=255)

         ggCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGc                                                                                 >  2:898392/1‑93 (MQ=255)

            tGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATg                                                                              <  1:3445231/93‑1 (MQ=255)

             gggggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATgg                                                                             >  2:200265/1‑93 (MQ=255)

              ggggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGa                                                                            >  1:2944597/1‑93 (MQ=255)

               gggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGAc                                                                           <  1:381347/93‑1 (MQ=255)

               gggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGAc                                                                           <  1:337484/93‑1 (MQ=255)

                ggTAAACCTGGCCTCCGGCGCTAACGCAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                          <  2:4221018/93‑1 (MQ=255)

                 gTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTg                                                                         <  2:4606748/93‑1 (MQ=255)

                   aaaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGa                                                                       >  2:2365012/1‑93 (MQ=255)

                      aCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATa                                                                    >  2:1405999/1‑93 (MQ=255)

                      aCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATa                                                                    <  2:854524/93‑1 (MQ=255)

                      aCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATa                                                                    >  2:3303079/1‑93 (MQ=255)

                       cTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATAc                                                                   <  2:4515387/93‑1 (MQ=255)

                        tGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACa                                                                  >  2:3909433/1‑93 (MQ=255)

                        tGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACa                                                                  <  1:3818978/93‑1 (MQ=255)

                         ggCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACaa                                                                 >  1:4082542/1‑93 (MQ=255)

                          gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                                >  1:3260415/1‑93 (MQ=255)

                          gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                                >  1:3260416/1‑93 (MQ=255)

                           ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCGTGGACTGGAATACAAgg                                                               <  2:693324/93‑1 (MQ=255)

                           ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAgg                                                               <  2:693325/93‑1 (MQ=255)

                           ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAgg                                                               >  2:1927364/1‑93 (MQ=255)

                            cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                              <  1:720078/93‑1 (MQ=255)

                            cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                              <  2:2685102/93‑1 (MQ=255)

                            cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                              <  1:510560/93‑1 (MQ=255)

                            cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTTCCGAATGGACTGGAATACAAAGt                                                              <  2:3288400/93‑1 (MQ=255)

                             tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCTGTATTTCCCGCATGGACTGGAATACAAGGTg                                                             >  2:1847777/1‑93 (MQ=255)

                             tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                             <  2:2381047/93‑1 (MQ=255)

                               cGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGc                                                           <  1:4068852/93‑1 (MQ=255)

                                ggCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGCATTTCCCGCATGGACTGGAATACAAGGTGGCg                                                          >  1:2869035/1‑93 (MQ=255)

                                 gcgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGt                                                         <  2:3741959/93‑1 (MQ=255)

                                   gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                       >  2:2872474/1‑93 (MQ=255)

                                   gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                       <  1:2668286/93‑1 (MQ=255)

                                    cTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATg                                                      <  1:3902507/93‑1 (MQ=255)

                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                     <  1:1493542/93‑1 (MQ=255)

                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaa                                                    >  2:3653060/1‑93 (MQ=255)

                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaa                                                    >  1:3873643/1‑93 (MQ=255)

                                       aCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaaa                                                   <  2:216457/93‑1 (MQ=255)

                                        cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGTATACAAGGTGGCGTATGAAac                                                  <  1:2090736/93‑1 (MQ=255)

                                        cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAac                                                  >  1:481969/1‑93 (MQ=255)

                                         gAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacc                                                 <  2:4203808/93‑1 (MQ=255)

                                           aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                               >  1:2827502/1‑93 (MQ=255)

                                           aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                               <  1:440954/93‑1 (MQ=255)

                                            aTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccacc                                              <  1:3408563/93‑1 (MQ=255)

                                             tGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCt                                             >  1:2439351/1‑93 (MQ=255)

                                                  aCAGCGGAGCTGGTGCTCAATCGTCTCCACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGttt                                        <  2:3819923/93‑1 (MQ=255)

                                                   cAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTg                                       <  1:1855868/93‑1 (MQ=255)

                                                   cAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTg                                       <  2:2583171/93‑1 (MQ=255)

                                                     gCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGtt                                     >  1:1858931/1‑93 (MQ=255)

                                                     gCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGtt                                     >  2:1595784/1‑93 (MQ=255)

                                                        gAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaaa                                  <  2:565285/93‑1 (MQ=255)

                                                        gAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaaa                                  <  2:3291917/93‑1 (MQ=255)

                                                        gAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaaa                                  >  2:2856063/1‑93 (MQ=255)

                                                        gAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaaa                                  <  2:3412478/93‑1 (MQ=255)

                                                         aGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTTTTAAAg                                 <  1:4498542/93‑1 (MQ=255)

                                                           cTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGcc                               <  2:2637798/93‑1 (MQ=255)

                                                           cTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGcc                               <  2:3680484/93‑1 (MQ=255)

                                                           cTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGcc                               <  2:160202/93‑1 (MQ=255)

                                                            tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                              >  1:4567478/1‑93 (MQ=255)

                                                            tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                              >  2:1835105/1‑93 (MQ=255)

                                                             ggTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTc                             >  2:1203280/1‑93 (MQ=255)

                                                             ggTGCTCAATCGTCTCGACGAGCTGGCGCGGGATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTc                             >  1:2602020/1‑93 (MQ=255)

                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                            >  1:3082153/1‑93 (MQ=255)

                                                               tGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCa                           <  2:2831543/93‑1 (MQ=255)

                                                                tcTCAAGCGTCTCGATGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACAACCTCGTTTGTTAAAGCCTCCAt                          <  1:739175/92‑1 (MQ=255)

                                                                  tCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTg                        >  1:305922/1‑93 (MQ=255)

                                                                     aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                     <  1:2872983/93‑1 (MQ=255)

                                                                      tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                    <  1:1155125/93‑1 (MQ=255)

                                                                      tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                    >  1:567224/1‑93 (MQ=255)

                                                                       cGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGAc                   <  1:4511224/93‑1 (MQ=255)

                                                                        gTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACg                  >  1:602875/1‑93 (MQ=255)

                                                                          cgcgacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACATGGTGGCGTATGAACCCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtg                <  1:125029/91‑1 (MQ=255)

                                                                           tcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtgg               >  2:4639401/1‑93 (MQ=255)

                                                                           tcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtgg               >  2:3577899/1‑93 (MQ=255)

                                                                           tcGACGAGCTGGCGCGGTATTTCCAGCATGGACTGGAATACAAGGAGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtgg               >  2:2418493/1‑93 (MQ=255)

                                                                            cgacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggt              >  2:4683600/1‑93 (MQ=255)

                                                                             gacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGACGACgtgttg             >  2:1006914/1‑93 (MQ=255)

                                                                             gacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggtg             >  1:957605/1‑93 (MQ=255)

                                                                              acgaGCTGGCGCGGGATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGa            >  2:2173722/1‑93 (MQ=255)

                                                                                gaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaa          >  2:1101211/1‑93 (MQ=255)

                                                                                gaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaa          <  1:3337341/93‑1 (MQ=255)

                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa         <  1:1571581/93‑1 (MQ=255)

                                                                                  gCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTCAAGCCTCCATTGAAGACGTGGTGacacc        >  1:1997886/1‑89 (MQ=255)

                                                                                    tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGATACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc      <  2:2811588/93‑1 (MQ=255)

                                                                                     ggCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgct     >  2:3948349/1‑93 (MQ=255)

                                                                                     ggCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgct     <  2:2618411/93‑1 (MQ=255)

                                                                                       cgcgGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctgc   >  2:3860827/1‑93 (MQ=255)

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCT  >  NC_000913/2586468‑2586648

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