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	99.6%	795.5 / ‑6.1	252	S863P (TCC→CCC)	acrB	acridine efflux pump
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (128/123); total (128/124)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 260 total aligned reads displayed.

CGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCG  >  NC_000913/480951‑481130

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

cGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCa                                                                                         >  1:2729774/1‑93 (MQ=255)

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

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

  tACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATc                                                                                       >  1:5177422/1‑93 (MQ=255)

  tACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATc                                                                                       >  1:7986665/1‑93 (MQ=255)

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

   aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATcc                                                                                      >  2:1817202/1‑93 (MQ=255)

   aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATcc                                                                                      >  2:5319033/1‑93 (MQ=255)

   aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATcc                                                                                      >  2:599228/1‑93 (MQ=255)

   aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCACGTACAGGGAAGGTGCCTGGATGCCTGAGAGACGTTGCTGATAGGGCATcc                                                                                      >  2:6985198/1‑93 (MQ=255)

    cAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATccc                                                                                     >  1:1777090/1‑93 (MQ=255)

    cAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATccc                                                                                     <  2:5668130/93‑1 (MQ=255)

    cAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATccc                                                                                     >  1:434548/1‑93 (MQ=255)

           gccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAg                                                                              >  1:6624290/1‑93 (MQ=255)

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

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

            ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGt                                                                             <  1:4735690/93‑1 (MQ=255)

            ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGt                                                                             <  1:609665/93‑1 (MQ=255)

            ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGt                                                                             <  2:9567514/93‑1 (MQ=255)

             cAGACACAGGAACACGACAATCAACTAAATCGCGTACAGTGAAGGTGCCTGGTTTCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTc                                                                            <  1:3283403/93‑1 (MQ=255)

             cAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTc                                                                            >  1:1609294/1‑93 (MQ=255)

              aGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCa                                                                           <  2:982069/93‑1 (MQ=255)

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

               gACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCAt                                                                          >  2:8525907/1‑93 (MQ=255)

               gACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCAt                                                                          <  1:9231254/93‑1 (MQ=255)

                acacAGGAACACGACAATCAACGAAGTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATa                                                                         <  1:5928272/93‑1 (MQ=255)

                  acaGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGc                                                                       <  1:4394449/93‑1 (MQ=255)

                  acaGGAACAAGACAATCAACGAAATCGCGTACAATGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGAGCATCCCCGTCCAGTCATAGc                                                                       <  2:8965682/93‑1 (MQ=255)

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

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

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

                    aGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGCTAGGGCATCCCCGTCCAGTCATAGCCa                                                                     >  1:2439282/1‑93 (MQ=255)

                     ggAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTAATAGCCaa                                                                    >  2:8796864/1‑93 (MQ=255)

                      gAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAac                                                                   >  1:2058062/1‑93 (MQ=255)

                       aaCACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAaca                                                                  >  1:618910/1‑93 (MQ=255)

                        acacGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAacac                                                                 <  1:1023432/93‑1 (MQ=255)

                         aacGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACAcc                                                                <  2:8804618/92‑1 (MQ=255)

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

                            gACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCgga                                                             >  1:7230304/1‑92 (MQ=255)

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

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

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

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

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

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

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

                                  cAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAg                                                       <  2:8751250/93‑1 (MQ=255)

                                  cAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAg                                                       <  1:885552/93‑1 (MQ=255)

                                    aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGtt                                                     <  2:7317011/93‑1 (MQ=255)

                                    aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGtt                                                     >  1:4106442/1‑93 (MQ=255)

                                       aaaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGc                                                  >  2:3716612/1‑93 (MQ=255)

                                       aaaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGc                                                  <  2:3665612/93‑1 (MQ=255)

                                         aTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACTCCGGTAGGCAGTTTGCTc                                                <  1:9537812/93‑1 (MQ=255)

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

                                         aTCGCGGACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTc                                                <  2:1525985/93‑1 (MQ=255)

                                           cgcgTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGc                                              <  1:2876133/93‑1 (MQ=255)

                                            gcgTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGcc                                             <  1:3803395/93‑1 (MQ=255)

                                             cgTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCa                                            <  2:7694139/93‑1 (MQ=255)

                                               tACAGTGAAGGTGCCTGGTTGCCTGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgt                                          <  2:7222251/93‑1 (MQ=255)

                                               tACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgt                                          >  1:8034891/1‑93 (MQ=255)

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

                                                aCAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgtt                                         >  2:9191472/1‑93 (MQ=255)

                                                 cAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttg                                        <  2:399370/93‑1 (MQ=255)

                                                  aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgt                                       <  1:10023029/93‑1 (MQ=255)

                                                  aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgt                                       <  1:4669354/93‑1 (MQ=255)

                                                   gTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgtt                                      <  1:3414602/93‑1 (MQ=255)

                                                   gTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgtt                                      <  1:3275679/93‑1 (MQ=255)

                                                     gAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTcc                                    >  2:8305570/1‑93 (MQ=255)

                                                     gAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCAAAGCCAACACCGGTGGGCAGTTTGCTCGCCAGTTGTTcc                                    <  2:1734502/93‑1 (MQ=255)

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

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

                                                        ggTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATc                                 >  2:4988268/1‑93 (MQ=255)

                                                         gTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCa                                >  2:4805946/1‑93 (MQ=255)

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

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

                                                          tGCCTGGTTGCCGGAGAGACGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAg                               >  2:2989011/1‑93 (MQ=255)

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

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

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

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

                                                               ggTTGCCGGAGAGACGGTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCa                          >  1:420720/1‑93 (MQ=255)

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

                                                                         gagaCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCAcc                >  1:1463423/1‑93 (MQ=255)

                                                                         gagaCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCAcc                >  2:449381/1‑93 (MQ=255)

                                                                          agaCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCg               <  1:6899120/93‑1 (MQ=255)

                                                                          agaCGTCCCTGATAGGGCATCCCCGGCCAGTCACAGCCAACACCGGTAGGCAGTTTGCTCGCCAGGCGTTCCATCAGCGCCATTGCGTCACCg               >  1:8942019/1‑93 (MQ=255)

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

                                                                            aCGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGt             <  2:8472838/93‑1 (MQ=255)

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

                                                                             cGTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTa            <  1:5414737/93‑1 (MQ=255)

                                                                              gTTCCTGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTAc           <  2:3942690/93‑1 (MQ=255)

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

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

                                                                                   tGATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTa      >  2:8059183/1‑93 (MQ=255)

                                                                                    gATAGGGCATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAc     >  2:5786362/1‑93 (MQ=255)

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

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

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

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

CGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCG  >  NC_000913/480951‑481130

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