BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	481,041	A→G	S863P (TCC→CCC)	acrB ←	acridine efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	481,041	0	A	G	97.8%	427.0 / ‑0.6	139	S863P (TCC→CCC)	acrB	acridine efflux pump
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (69/67); total (70/69)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 143 total aligned reads displayed.

CGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCG  >  NC_000913/480951‑481130

cGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCa                                                                                         <  2:2807846/93‑1 (MQ=255)

cGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCAAGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGAGAGTAcc                                                                                         >  2:1594154/1‑92 (MQ=255)

 gTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCAt                                                                                        >  2:4557674/1‑93 (MQ=255)

 gTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCAt                                                                                        >  1:2070060/1‑93 (MQ=255)

 gTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCAt                                                                                        >  1:2070059/1‑93 (MQ=255)

  tACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGTTTGCCGCAGAGACGCTCCTCATAGGGCATc                                                                                       >  1:688801/1‑93 (MQ=37)

  tACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATc                                                                                       <  2:4314728/93‑1 (MQ=255)

  tACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATc                                                                                       <  1:4517926/93‑1 (MQ=255)

   aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATcc                                                                                      <  1:1122135/93‑1 (MQ=255)

      gcgcCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCg                                                                                   >  1:1220786/1‑93 (MQ=255)

       cgccgcCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGt                                                                                  >  2:4542673/1‑93 (MQ=255)

         ccgccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTcc                                                                                <  2:1207781/93‑1 (MQ=255)

         ccgccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTcc                                                                                <  2:778727/93‑1 (MQ=255)

         ccgccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTcc                                                                                <  1:1214640/93‑1 (MQ=255)

          cgccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCa                                                                               >  1:2124068/1‑93 (MQ=255)

           gccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAg                                                                              <  1:4012833/93‑1 (MQ=255)

            ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGt                                                                             >  1:2331782/1‑93 (MQ=255)

             cAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGTAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTc                                                                            >  2:3940317/1‑93 (MQ=255)

              aGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCa                                                                           >  1:1470955/1‑93 (MQ=255)

              aGACACAGGAACACGACAATCAACGAAATCGCGGACAGTGAAGGCGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCa                                                                           <  1:4577588/93‑1 (MQ=255)

                acacAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATa                                                                         <  1:3802162/93‑1 (MQ=255)

                   caGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGcg                                                                      >  2:3639576/1‑92 (MQ=255)

                   caGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGcc                                                                      <  2:3497177/93‑1 (MQ=255)

                    aGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCa                                                                     >  2:1957780/1‑93 (MQ=255)

                     ggAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCaa                                                                    >  1:1619541/1‑93 (MQ=255)

                      gAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAac                                                                   <  1:1057240/93‑1 (MQ=255)

                          acGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCg                                                               >  1:4313739/1‑93 (MQ=255)

                           cGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCgg                                                              <  2:811006/93‑1 (MQ=255)

                           cGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCgg                                                              >  2:3012274/1‑93 (MQ=255)

                            gACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGt                                                             <  1:3661116/93‑1 (MQ=255)

                            gACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGt                                                             <  2:234979/93‑1 (MQ=255)

                            gACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGt                                                             >  2:731255/1‑93 (MQ=255)

                             aCAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTa                                                            <  2:597899/93‑1 (MQ=255)

                             aCAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTa                                                            <  2:2999454/93‑1 (MQ=255)

                              cAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAg                                                           <  1:4695978/93‑1 (MQ=255)

                              cAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAg                                                           >  1:1492055/1‑93 (MQ=255)

                               aaTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAgg                                                          >  1:3310409/1‑93 (MQ=255)

                               aaTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAgg                                                          >  1:2904/1‑93 (MQ=255)

                                aTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGc                                                         >  1:3585601/1‑93 (MQ=255)

                                aTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGc                                                         >  2:1579173/1‑93 (MQ=255)

                                aTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGc                                                         <  2:3312653/93‑1 (MQ=255)

                                 tCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCa                                                        >  1:612437/1‑93 (MQ=255)

                                 tCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCa                                                        >  1:3390741/1‑93 (MQ=255)

                                 tCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCa                                                        <  1:2623985/93‑1 (MQ=255)

                                  cAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGAAGGCAg                                                       >  1:3743110/1‑93 (MQ=255)

                                   atggAAAGCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGt                                                      <  1:2695566/90‑1 (MQ=255)

                                   aaCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGt                                                      >  1:1095793/1‑93 (MQ=255)

                                      gAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTg                                                   >  2:2163768/1‑93 (MQ=255)

                                        aaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCt                                                 <  2:1725503/93‑1 (MQ=255)

                                         aTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTc                                                >  2:1268866/1‑93 (MQ=255)

                                             cgTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCa                                            >  2:2751609/1‑93 (MQ=255)

                                              gTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAg                                           >  2:3200174/1‑93 (MQ=255)

                                                aCAGTGAGGGTGCCTGGTTTCCGGAGATACGTTCCTTATATGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgtt                                         <  1:2316671/93‑1 (MQ=255)

                                                aCAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgtt                                         <  2:1142267/93‑1 (MQ=255)

                                                   gTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgtt                                      >  2:349677/1‑93 (MQ=255)

                                                     gAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACTCCGGTAGGCAGTTTGCTCGCCAGTTGTTcc                                    <  2:3882455/93‑1 (MQ=255)

                                                      aaGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCa                                   <  2:3694600/93‑1 (MQ=255)

                                                      aaGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCGCGCCAGTTGTTCCa                                   <  2:3694599/93‑1 (MQ=255)

                                                         gTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCa                                <  2:3766569/93‑1 (MQ=255)

                                                         gTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCa                                <  2:1936196/93‑1 (MQ=255)

                                                          tGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAg                               <  1:4384195/93‑1 (MQ=255)

                                                           gCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGc                              <  2:1319836/93‑1 (MQ=255)

                                                           gCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGc                              <  2:2158811/93‑1 (MQ=255)

                                                             cTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                            >  2:2844692/1‑93 (MQ=255)

                                                             cTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                            <  2:911293/93‑1 (MQ=255)

                                                             cTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                            >  2:1727387/1‑93 (MQ=255)

                                                             cTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                            >  1:3368291/1‑93 (MQ=255)

                                                              tGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTcc                           >  2:3229031/1‑93 (MQ=255)

                                                              tGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTcc                           >  1:3600168/1‑93 (MQ=255)

                                                              tGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTca                           >  2:3329149/1‑92 (MQ=255)

                                                               ggTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCa                          <  2:3664494/93‑1 (MQ=255)

                                                               ggTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCa                          >  1:3949071/1‑93 (MQ=255)

                                                               ggTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCa                          <  2:465331/93‑1 (MQ=255)

                                                                gTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAt                         >  2:1163964/1‑93 (MQ=255)

                                                                 ttGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAtt                        <  1:3469769/93‑1 (MQ=255)

                                                                 ttGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAtt                        <  2:1301019/93‑1 (MQ=255)

                                                                 ttGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAtt                        <  1:1359406/93‑1 (MQ=255)

                                                                 gtGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAtt                        <  1:2608054/92‑1 (MQ=255)

                                                                    ccGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCt                     <  2:1690635/93‑1 (MQ=255)

                                                                    ccGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCt                     <  1:1075542/93‑1 (MQ=255)

                                                                     cGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCtt                    >  1:4579104/1‑93 (MQ=255)

                                                                      ggAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTc                   >  1:300163/1‑93 (MQ=255)

                                                                      ggAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTc                   >  2:2347971/1‑93 (MQ=255)

                                                                      ggAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTc                   <  1:3682082/93‑1 (MQ=255)

                                                                       gagagaCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCa                  >  2:596418/1‑93 (MQ=255)

                                                                           gaCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCgg              >  1:4653310/1‑93 (MQ=255)

                                                                            aCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGt             >  1:3804531/1‑93 (MQ=255)

                                                                               ttCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACt          >  1:1033695/1‑93 (MQ=255)

                                                                                tCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtt         >  2:124456/1‑93 (MQ=255)

                                                                                 ccTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACttt        <  2:4302290/93‑1 (MQ=255)

                                                                                  cTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtttt       <  1:4275554/93‑1 (MQ=255)

                                                                                  cTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtttt       >  1:2585178/1‑93 (MQ=255)

                                                                                  cTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtttt       >  2:1481017/1‑93 (MQ=255)

                                                                                  cTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtttt       <  2:2028256/93‑1 (MQ=255)

                                                                                   tGATAGGGCATACCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTa      >  2:3356752/1‑93 (MQ=255)

                                                                                     aTAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAcc    <  2:1677284/93‑1 (MQ=255)

                                                                                      tAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAccc   >  1:3928575/1‑93 (MQ=255)

                                                                                      tAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAccc   <  1:1620334/93‑1 (MQ=255)

                                                                                       aGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCg  >  1:4464096/1‑93 (MQ=255)

                                                                                       aGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCg  >  1:2191218/1‑93 (MQ=255)

CGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCG  >  NC_000913/480951‑481130

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: `‑`