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	99.3%	439.1 / ‑4.8	140	Q315R (CAG→CGG)	acrD	AcrD
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (65/74); total (65/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

Only 100 of 144 total aligned reads displayed.

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586468‑2586647

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

  gccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       >  2:1003313/1‑93 (MQ=255)

  gacTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       >  2:2917507/3‑93 (MQ=255)

      ccGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTcc                                                                                   <  1:1133520/93‑1 (MQ=255)

       cGGTCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTccc                                                                                  <  1:4018280/93‑1 (MQ=255)

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

           gtGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCAt                                                                              <  1:2290741/92‑1 (MQ=255)

            tGGGGGTAAAACTGGCCTCCGGCGATAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGGCTCGACGAGCTGGCGCGGTATTTCCCGCATg                                                                             <  2:3918859/93‑1 (MQ=255)

                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                         <  1:2675920/93‑1 (MQ=255)

                 gTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTg                                                                        >  1:430873/1‑93 (MQ=255)

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

                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGaa                                                                     >  2:4404614/1‑93 (MQ=255)

                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGaa                                                                     <  2:1854940/93‑1 (MQ=255)

                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGaa                                                                     <  2:2405882/93‑1 (MQ=255)

                     aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    >  2:2598376/1‑93 (MQ=255)

                     aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    <  1:511949/93‑1 (MQ=255)

                     aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    >  1:952372/1‑93 (MQ=255)

                     aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    <  2:3651010/93‑1 (MQ=255)

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

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

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

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

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

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

                         ggCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACaa                                                                <  2:2320516/93‑1 (MQ=255)

                         ggCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACaa                                                                <  1:3695637/93‑1 (MQ=255)

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

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

                            cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                             >  2:2051/1‑93 (MQ=255)

                              ccGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTgg                                                           >  2:2090679/1‑93 (MQ=255)

                               cGGCGCTAACGAAAGGGCGACAGCGGAGCTTGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGc                                                          <  1:329738/93‑1 (MQ=255)

                                 gcgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGt                                                        >  2:1265222/1‑93 (MQ=255)

                                  cgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTa                                                       <  2:811156/93‑1 (MQ=255)

                                  cgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTa                                                       <  2:3921216/93‑1 (MQ=255)

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

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

                                    cTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATg                                                     >  1:2905547/1‑93 (MQ=255)

                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    >  2:4038676/1‑93 (MQ=255)

                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    >  1:4534624/1‑93 (MQ=255)

                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    >  1:4385349/1‑93 (MQ=255)

                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaa                                                   <  1:3160031/93‑1 (MQ=255)

                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGAACGGAATACAAGTTGGCGTagaga                                                   >  1:3286921/1‑89 (MQ=255)

                                          aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                               <  2:2663410/93‑1 (MQ=255)

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

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

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

                                              ggCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTc                                           <  2:3527941/93‑1 (MQ=255)

                                               gCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCg                                          <  1:3769216/93‑1 (MQ=255)

                                                  aCAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGttt                                       <  1:2202874/93‑1 (MQ=255)

                                                  aCAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGttt                                       <  1:1560233/93‑1 (MQ=255)

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

                                                       ggAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaa                                  <  2:1867158/93‑1 (MQ=255)

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

                                                            tGGTGCTCAATCGTCTTGACGAGCTGGCGGGGTATTTCCCGCATGGACTGGAATACGAGGTAGCGTATGGAACCACATCGTTTGTTAAAGCCt                             >  2:3182442/1‑93 (MQ=255)

                                                            tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                             <  2:86681/93‑1 (MQ=255)

                                                            tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                             <  2:4335187/93‑1 (MQ=255)

                                                             ggTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTc                            <  2:1972378/93‑1 (MQ=255)

                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGGATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                           >  1:82531/1‑93 (MQ=255)

                                                                gCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGACTCCAt                         >  2:3421505/1‑93 (MQ=255)

                                                                 cttaaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                        <  2:3809867/90‑1 (MQ=255)

                                                                 cTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                        <  2:4392827/93‑1 (MQ=255)

                                                                 cTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                        <  2:4299716/93‑1 (MQ=255)

                                                                  tCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTg                       <  2:4306647/93‑1 (MQ=255)

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

                                                                   cAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGa                      >  2:3070818/1‑93 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  2:3020982/93‑1 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  1:3153681/93‑1 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  1:4252861/93‑1 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  1:4609143/93‑1 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTGCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  1:1814563/93‑1 (MQ=255)

                                                                     aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                    >  1:1891339/1‑93 (MQ=255)

                                                                      tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACACGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                   <  1:3040032/93‑1 (MQ=255)

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

                                                                      acGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                   <  1:3664275/92‑1 (MQ=255)

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

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

                                                                        gTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACg                 <  2:1465783/93‑1 (MQ=255)

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

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

                                                                         tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                <  1:315782/93‑1 (MQ=255)

                                                                         tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                >  1:3134692/1‑93 (MQ=255)

                                                                         tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                >  1:2718825/1‑93 (MQ=255)

                                                                          ctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtg               <  2:4347864/93‑1 (MQ=255)

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

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

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

                                                                           tcGACGAGCTGGCGCGGTATTCCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtgg              <  2:2481410/93‑1 (MQ=255)

                                                                            cgacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggt             <  1:3938109/93‑1 (MQ=255)

                                                                            cgacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggt             <  1:3143994/93‑1 (MQ=255)

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

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

                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa        >  2:388647/1‑93 (MQ=255)

                                                                                  gCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAAc       >  1:4438918/1‑93 (MQ=255)

                                                                                  gCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAAc       >  2:4400385/1‑93 (MQ=255)

                                                                                   cTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGCTTGTTAAAGCCTCCATTGAAGACGTGCTGACAACg      >  1:2105559/1‑93 (MQ=255)

                                                                                    tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     >  2:4135778/1‑93 (MQ=255)

                                                                                    tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     >  1:2926592/1‑93 (MQ=255)

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

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

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

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586468‑2586647

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