BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,421	+T	coding (504/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,421	1	.	T	99.2%	437.0 / ‑3.7	119	coding (504/516 nt)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (60/58); total (61/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATT‑G‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                                >  2:1132533/1‑93 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                                >  2:6713870/1‑93 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAGGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                                <  1:2981338/93‑1 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                               >  1:990386/1‑92 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                               <  1:1567329/93‑2 (MQ=255)

 cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                               <  1:1673509/93‑2 (MQ=255)

   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                             >  2:6385580/1‑90 (MQ=255)

   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                             >  2:5544515/1‑90 (MQ=255)

   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                             >  2:3952663/1‑90 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                            <  2:6498533/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                            <  1:5372490/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                            <  2:4370028/93‑5 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggg                                                                            >  2:1014114/1‑89 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                                      <  1:3213125/93‑1 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                                      <  1:6058621/93‑1 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                                      <  2:3791411/93‑1 (MQ=255)

            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGAt                                                                    <  2:5203299/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATg                                                                   >  1:5841435/1‑93 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATg                                                                   >  1:4191617/1‑93 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCTAAGA‑‑‑AACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                                  <  2:912291/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                                  <  1:31479/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                                 >  2:806217/1‑93 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                                 <  2:3841035/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                                 <  1:47385/93‑1 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacg                                                               >  2:2993890/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacg                                                               >  2:274444/1‑93 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacga                                                              <  2:3073690/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                             >  2:5179197/1‑93 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                             <  1:5688350/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                             <  2:2367591/93‑1 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgacg                                                            >  1:2108765/1‑93 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgacg                                                            <  2:1665924/93‑1 (MQ=255)

                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGc                                                           <  1:4782314/93‑1 (MQ=255)

                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGc                                                           >  2:777217/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGG‑ATTGCGTTTGGTGTCGGGATGAACGACGcc                                                          >  2:6280598/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  2:6297100/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          <  2:3562144/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  2:2100640/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  2:4869242/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  1:3137376/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          <  1:4708802/93‑1 (MQ=255)

                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGa                                                        <  1:6379746/89‑1 (MQ=255)

                        ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGa                                                        <  1:5667309/89‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                      <  2:5979161/91‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                      <  1:2934883/91‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                      <  2:1507336/91‑1 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                      <  2:3475704/91‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                     <  2:5637763/92‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                     <  2:3991366/92‑1 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                    >  1:721326/1‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                    <  2:6010772/93‑1 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                    >  1:2205579/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   <  1:6289377/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   <  2:5981706/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   <  1:3895288/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  2:2767946/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  2:2939531/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  <  2:3264646/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  2:6691263/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGAAGCCGAAATGcg                                                  >  1:1083384/1‑92 (MQ=255)

                                                               actgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg                    >  1:1339363/5‑93 (MQ=255)

                                                               actgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg                    <  2:2395438/89‑1 (MQ=255)

                                                                ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGCTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   >  2:4704086/4‑93 (MQ=255)

                                                                ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   <  2:6247674/90‑1 (MQ=255)

                                                                ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   <  2:3666627/90‑1 (MQ=255)

                                                                ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   <  1:1812278/90‑1 (MQ=255)

                                                                 tgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCa                  <  1:1193761/91‑1 (MQ=255)

                                                                   ggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                <  1:4147958/93‑1 (MQ=255)

                                                                   ggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                <  2:1739758/93‑1 (MQ=255)

                                                                    gCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg               >  1:3847057/1‑93 (MQ=255)

                                                                    gCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg               >  1:1354706/1‑93 (MQ=255)

                                                                     cTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc              >  2:3390041/1‑92 (MQ=255)

                                                                     cTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc              >  1:2566403/1‑92 (MQ=255)

                                                                      tACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             >  1:269785/1‑91 (MQ=255)

                                                                      tACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             <  2:5981713/93‑3 (MQ=255)

                                                                          gccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  2:4888485/3‑93 (MQ=255)

                                                                          gccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  2:2154263/3‑93 (MQ=255)

                                                                          gccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  2:1566474/3‑93 (MQ=255)

                                                                          gccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  1:6279783/3‑93 (MQ=255)

                                                                          gccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  2:1796222/3‑93 (MQ=255)

                                                                           ccTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca        >  1:5192517/2‑93 (MQ=255)

                                                                            cTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg       >  2:5673178/1‑93 (MQ=255)

                                                                             tGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc      >  1:6343208/1‑93 (MQ=255)

                                                                             tGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc      >  1:5044605/1‑93 (MQ=255)

                                                                              gAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg     <  1:1865152/93‑1 (MQ=255)

                                                                              gAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg     >  2:3227353/1‑93 (MQ=255)

                                                                               aaGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGt    >  1:1171613/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCCTTATCGGCAGTGCGCACCAGCGTc   >  2:2894029/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:5446439/93‑1 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:4562313/93‑1 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:3983998/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:3590897/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:3311561/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:2222735/1‑93 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  1:950398/93‑1 (MQ=255)

                                                                                aGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  1:852417/1‑93 (MQ=255)

                                                                                 ggATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  1:2578659/1‑93 (MQ=255)

                                                                                 ggATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  1:3088337/1‑93 (MQ=255)

                                                                                 ggATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  1:3428999/1‑93 (MQ=255)

                                                                                 ggATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  <  2:3544748/93‑1 (MQ=255)

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATT‑G‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503

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