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.6%	808.0 / ‑6.1	245	Q315R (CAG→CGG)	acrD	AcrD
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (98/146); minor base T (0/1); total (98/147)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 250 total aligned reads displayed.

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586468‑2586647

cagccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGt                                                                                         >  2:4024311/1‑93 (MQ=255)

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

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

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

  gccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       <  1:935491/93‑1 (MQ=255)

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

   ccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAtt                                                                                      >  1:149971/1‑93 (MQ=255)

   ccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAAGGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAtt                                                                                      <  2:5023249/93‑1 (MQ=255)

     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                    >  2:2918964/1‑93 (MQ=255)

     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                    <  2:3865273/93‑1 (MQ=255)

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

            tGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCCTg                                                                             >  1:5779956/1‑93 (MQ=255)

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

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

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

                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                         >  2:5747452/1‑93 (MQ=255)

                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                         >  1:3823435/1‑93 (MQ=255)

                 gTAAAGCTCACCGCCGTCGCTAACGAAATGCCGACAGCGGAGCTGGTGCTCAATCGCCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTg                                                                        <  2:5632818/93‑1 (MQ=255)

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

                  tAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTgg                                                                       <  2:5387972/93‑1 (MQ=255)

                   aaaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGa                                                                      <  2:2126115/93‑1 (MQ=255)

                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGTGGTATTTCCCGCATGGACTGGaa                                                                     >  2:1057095/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                          gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                               <  1:5584078/93‑1 (MQ=255)

                             tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                            >  1:777935/1‑93 (MQ=255)

                              ccGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTgg                                                           <  2:135511/93‑1 (MQ=255)

                                ggCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCg                                                         <  2:5968309/93‑1 (MQ=255)

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

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

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

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

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

                                    cTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTTCCGCATGGACTGGAATACAAGGTGGCGTATg                                                     >  1:4948607/1‑93 (MQ=255)

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

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

                                       aCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaaa                                                  >  2:6600833/1‑93 (MQ=255)

                                       aCGAAAGGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaaa                                                  <  1:210087/93‑1 (MQ=255)

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

                                        cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAac                                                 <  2:42598/93‑1 (MQ=255)

                                         gAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacc                                                >  2:4328769/1‑93 (MQ=255)

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

                                            aTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccacc                                             >  1:6202640/1‑93 (MQ=255)

                                             tGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCt                                            <  1:4684502/93‑1 (MQ=255)

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

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

                                             aggCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCt                                            <  1:4976037/92‑1 (MQ=255)

                                                cGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGCGGCGTATGAAACCACCTCGt                                         <  2:4634291/93‑1 (MQ=255)

                                                  aCAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCTCCTCGttt                                       <  1:6505191/93‑1 (MQ=255)

                                                   aaGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTg                                      <  1:4784236/92‑1 (MQ=255)

                                                    aGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGt                                     <  2:3605998/93‑1 (MQ=255)

                                                      cGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTa                                   <  1:5050207/93‑1 (MQ=255)

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

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

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

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

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

                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                           <  1:3903039/93‑1 (MQ=255)

                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                           <  1:4260319/93‑1 (MQ=255)

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

                                                               tGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCa                          >  2:6015722/1‑93 (MQ=255)

                                                                gCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAt                         >  1:2675796/1‑93 (MQ=255)

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

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

                                                                 cTCAATCGTCTCGACGAGCTGGCGCGGCATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                        >  2:1753302/1‑93 (MQ=255)

                                                                    aaTCGTCTCGACGATCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGGATGAAACCACCTCGTTTGTTAAAGACTCCATTGaa                     >  2:2708673/1‑93 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     >  2:6441647/1‑93 (MQ=255)

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     >  2:95299/1‑93 (MQ=255)

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

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

                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     >  1:5713968/1‑93 (MQ=255)

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

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

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

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

                                                                       cGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGAc                  >  2:2368738/1‑93 (MQ=255)

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

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

                                                                          ctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtg               >  1:4627870/1‑93 (MQ=255)

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

                                                                               cgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaa          >  2:2431389/1‑93 (MQ=255)

                                                                               cgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaa          <  2:1888538/93‑1 (MQ=255)

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

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

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

                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAATGAAGACGTGGTGaaaa        >  2:1003356/1‑93 (MQ=255)

                                                                                  gCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAAc       <  1:5856861/93‑1 (MQ=255)

                                                                                  gCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAAc       <  2:2009950/93‑1 (MQ=255)

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

                                                                                       cgcgGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctgc  <  2:3691211/93‑1 (MQ=255)

                                                                                       cgcgGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctgc  <  2:4964837/93‑1 (MQ=255)

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586468‑2586647

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