BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,234,288	A→G	T150A (ACC→GCC)	b4027 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,234,288	0	A	G	99.5%	744.0 / ‑4.8	219	T150A (ACC→GCC)	b4027	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (104/114); total (104/115)
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.77e-01

Only 100 of 223 total aligned reads displayed.

GATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGC  >  NC_000913/4234198‑4234377

taTGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGcc                                                                                         <  2:1066758/92‑1 (MQ=255)

gATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGcc                                                                                         >  1:1289119/1‑93 (MQ=255)

gATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGcc                                                                                         <  2:5074659/93‑1 (MQ=255)

 gtGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCg                                                                                        <  1:2562269/92‑1 (MQ=255)

 aTGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCg                                                                                        >  2:3317351/1‑93 (MQ=255)

  tGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGa                                                                                       <  1:1597563/93‑1 (MQ=255)

   ggCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGAt                                                                                      <  1:6337902/93‑1 (MQ=255)

   ggCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGAt                                                                                      <  2:3084238/93‑1 (MQ=255)

     cgcAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATAc                                                                                    >  2:5268408/1‑93 (MQ=255)

      gcAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATAcc                                                                                   >  2:6224913/1‑93 (MQ=255)

      gcAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATAcc                                                                                   >  2:6320817/1‑93 (MQ=255)

      gcAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATAcc                                                                                   <  2:2877211/93‑1 (MQ=255)

       cAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCg                                                                                  <  1:4315062/93‑1 (MQ=255)

       cAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCg                                                                                  <  2:5624182/93‑1 (MQ=255)

        aaGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGt                                                                                 >  1:619856/1‑93 (MQ=255)

        aaGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGt                                                                                 <  1:2415067/93‑1 (MQ=255)

         aGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTc                                                                                <  1:213658/93‑1 (MQ=255)

         aGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTc                                                                                >  1:1434764/1‑93 (MQ=255)

         aGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTc                                                                                <  2:4719337/93‑1 (MQ=255)

          gCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCa                                                                               >  2:4215135/1‑93 (MQ=255)

            tGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCaaa                                                                             <  1:6429162/93‑1 (MQ=255)

             ggACGCGCACTATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGGCGATACCGTCAAAg                                                                            >  1:4418069/1‑93 (MQ=255)

             ggACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAg                                                                            >  1:5592742/1‑93 (MQ=255)

              gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAgc                                                                           >  2:2143772/1‑92 (MQ=255)

              gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGt                                                                           <  1:1495411/93‑1 (MQ=255)

              gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGt                                                                           <  1:5817284/93‑1 (MQ=255)

              gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGt                                                                           >  1:4541887/1‑93 (MQ=255)

              gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGt                                                                           >  1:1225923/1‑93 (MQ=255)

                     aCGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGc                                                                    <  1:6457094/93‑1 (MQ=255)

                      cGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCg                                                                   <  1:1175369/93‑1 (MQ=255)

                       gATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGa                                                                  >  1:1844841/1‑93 (MQ=255)

                       gATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGa                                                                  >  2:1793012/1‑93 (MQ=255)

                        aTGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTGAAAGTCGGCAGCGAt                                                                 >  1:2743116/1‑93 (MQ=255)

                        aTGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGAt                                                                 >  2:3760329/1‑93 (MQ=255)

                        aTGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGAt                                                                 <  2:207143/93‑1 (MQ=255)

                         tGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATg                                                                >  2:5404550/1‑93 (MQ=255)

                         tGGGCTGGACCGAGTACCAGCAGGTACGCGACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATg                                                                <  2:4237295/93‑1 (MQ=255)

                              tGGACCGAGTACCAGCAGGTACGCTACGCCGCCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAAcc                                                           <  1:1766736/93‑1 (MQ=255)

                              tGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAAcc                                                           >  2:5960601/1‑93 (MQ=255)

                               ggACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAAccc                                                          >  2:5543501/1‑93 (MQ=255)

                               ggACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAAccc                                                          <  2:4780209/93‑1 (MQ=255)

                               ggACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAAccc                                                          <  1:1302499/93‑1 (MQ=255)

                                gACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAAcccc                                                         <  2:5455191/93‑1 (MQ=255)

                                  ccGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCgg                                                       <  1:4825336/93‑1 (MQ=255)

                                  ccGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCgg                                                       >  2:2741232/1‑93 (MQ=255)

                                  ccGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCgg                                                       >  1:4650901/1‑93 (MQ=255)

                                  ccGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCgg                                                       >  2:5040274/1‑93 (MQ=255)

                                   cGAGTACCAGCAGGTACGCTACGCCACCGCACNCTNAGNNTTCAAATGGGATGGNGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGt                                                      <  1:5084862/93‑1 (MQ=255)

                                    gAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGtt                                                     >  2:2535359/1‑93 (MQ=255)

                                    gAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGtt                                                     >  1:1064904/1‑93 (MQ=255)

                                     aGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTc                                                    >  2:2700651/1‑93 (MQ=255)

                                      gTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTcg                                                   >  2:4030333/1‑93 (MQ=255)

                                      gTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTcg                                                   >  1:311460/1‑93 (MQ=255)

                                       tACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTcgc                                                  <  1:4600522/93‑1 (MQ=255)

                                          cagcagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTg                                               >  1:6247439/1‑93 (MQ=255)

                                           agcagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGc                                              <  2:3912695/93‑1 (MQ=255)

                                           agcagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGc                                              <  2:4492797/93‑1 (MQ=255)

                                            gcagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCt                                             >  1:5155839/1‑93 (MQ=255)

                                             cagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTg                                            <  2:127200/93‑1 (MQ=255)

                                              agGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTgg                                           <  1:5418453/93‑1 (MQ=255)

                                              agGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTgg                                           >  1:1938452/1‑93 (MQ=255)

                                              agGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGCATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAACCCCCGGTTCGCGTGCTgg                                           <  2:3969680/93‑1 (MQ=255)

                                               ggTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGa                                          >  2:16392/1‑93 (MQ=255)

                                                   cgctacgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACgaa                                      <  1:1448334/93‑1 (MQ=255)

                                                     ctacgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACgaaga                                    <  2:4641729/93‑1 (MQ=255)

                                                      tacgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACgaagaa                                   >  2:167573/1‑93 (MQ=255)

                                                       acgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGaagaag                                  >  2:4352558/1‑93 (MQ=255)

                                                       acgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGaagaag                                  >  1:4764866/1‑93 (MQ=255)

                                                        cgcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGt                                 >  1:2053247/1‑93 (MQ=255)

                                                         gcCACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGtc                                <  2:5832490/93‑1 (MQ=255)

                                                           cACCGCACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGtctc                              <  1:720072/93‑1 (MQ=255)

                                                                cACGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCg                         >  2:2373944/1‑93 (MQ=255)

                                                                cACGCCCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCg                         <  2:4975431/93‑1 (MQ=255)

                                                                 aCGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCgc                        >  2:4182091/1‑92 (MQ=255)

                                                                  cGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAc                       >  2:3514213/1‑93 (MQ=255)

                                                                  cGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAc                       <  1:2622679/93‑1 (MQ=255)

                                                                  cGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAc                       <  1:1700350/93‑1 (MQ=255)

                                                                  cGCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAc                       <  1:2299473/93‑1 (MQ=255)

                                                                   gCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAcc                      >  1:3432797/1‑93 (MQ=255)

                                                                   gCTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGAcc                      <  1:5122082/93‑1 (MQ=255)

                                                                    cTCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACccacca                     >  1:3328652/1‑93 (MQ=255)

                                                                     tCAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAATCCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAg                    <  2:2118126/93‑1 (MQ=255)

                                                                      cAGTCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAgg                   <  2:1905022/93‑1 (MQ=255)

                                                                        gTCTTCAAATGGGATGGCGCCGATCCCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGgcg                 <  1:3046608/93‑1 (MQ=255)

                                                                         tCTTCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGgcgc                >  2:3163388/1‑93 (MQ=255)

                                                                           ttCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGgcgcgc              <  1:2500368/93‑1 (MQ=255)

                                                                           ttCAAATGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGgcgcgc              <  2:1348530/93‑1 (MQ=255)

                                                                               aaTGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGc          >  1:9983/1‑93 (MQ=255)

                                                                               aaTGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGc          >  2:1496130/1‑93 (MQ=255)

                                                                               aaTGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGc          >  2:1034920/1‑93 (MQ=255)

                                                                                aTGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCa         >  1:3740663/1‑93 (MQ=255)

                                                                                 tGGGATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCAt        >  2:1411652/1‑93 (MQ=255)

                                                                                   ggATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATaa      >  1:5470814/1‑93 (MQ=255)

                                                                                    gATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGGGCGCTGGCATAAc     >  1:1734142/1‑93 (MQ=255)

                                                                                    gATGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAAc     <  1:4152721/93‑1 (MQ=255)

                                                                                     aTGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAAcc    <  2:3332115/93‑1 (MQ=255)

                                                                                     aTGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAAcc    >  2:2209097/1‑93 (MQ=255)

                                                                                      tGGCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCg   >  1:2865383/1‑93 (MQ=255)

                                                                                       ggCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGc  <  1:3592878/93‑1 (MQ=255)

                                                                                       ggCGCCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGc  >  2:645911/1‑93 (MQ=255)

GATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGC  >  NC_000913/4234198‑4234377

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