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	100.0%	458.9 / NA	117	L168F (TTG→TTT)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base T (48/69); total (48/69)

Only 100 of 120 total aligned reads displayed.

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                                                                >  1:1148257/1‑93 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                                                                <  1:2423527/93‑1 (MQ=255)

tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                                                                <  2:3052296/93‑1 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                                                                              <  1:6351818/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                                                                              >  1:5279004/1‑91 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                                                                              <  2:2720365/93‑3 (MQ=255)

  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                                                                              <  2:4945205/93‑3 (MQ=255)

   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttgg                                                                             <  1:1390840/93‑4 (MQ=255)

   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttgg                                                                             >  2:2368692/1‑90 (MQ=255)

    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑AAGAGACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttggg                                                                            >  2:3799169/1‑89 (MQ=255)

    gAAGTGATGGCGGGCGGCGTCTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttggt                                                                            <  2:2138944/93‑5 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                      <  2:3643840/93‑1 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                      <  1:4539782/93‑1 (MQ=255)

          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                      >  1:4895227/1‑93 (MQ=255)

           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                                                                     <  2:2652227/93‑1 (MQ=255)

            ggcgcgcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACT‑GGCCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                                                                    <  1:819012/93‑1 (MQ=255)

             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                                                                   <  2:4096718/93‑1 (MQ=255)

              cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                                                                  <  1:3726700/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                                                                 <  1:3501601/93‑1 (MQ=255)

               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                                                                 >  2:201803/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                                                                >  2:2412012/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                                                                <  1:4668012/93‑1 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                                                                <  2:2472633/93‑1 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                                                                >  2:4169501/1‑93 (MQ=255)

                ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                                                                >  2:1145412/1‑93 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                                                               <  1:3946254/93‑1 (MQ=255)

                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                                                               <  2:733206/93‑1 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                                                              <  1:3733356/93‑1 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                                                              <  1:5029312/93‑1 (MQ=255)

                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                                                              <  1:3061014/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             >  1:838731/1‑93 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             <  1:1387788/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             >  2:506903/1‑93 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             <  2:6002290/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             <  2:6349720/93‑1 (MQ=255)

                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                             <  1:2279083/93‑1 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg                                                            >  1:1938896/1‑93 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg                                                            >  1:550341/1‑93 (MQ=255)

                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg                                                            <  1:1704662/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  2:5795906/1‑93 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          <  2:3718784/93‑1 (MQ=255)

                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                          >  2:3054417/1‑93 (MQ=255)

                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                      <  1:2486652/91‑1 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                     >  2:2970238/2‑93 (MQ=255)

                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                     >  2:269946/2‑93 (MQ=255)

                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                    >  2:1611781/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   >  1:5171996/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   >  1:501772/1‑93 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   <  2:4329317/93‑1 (MQ=255)

                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                   >  1:1121834/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  2:2466575/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  2:6392472/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  <  2:871705/93‑1 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:3051662/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:5530256/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:616170/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:4138537/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:4308666/1‑93 (MQ=255)

                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                  >  1:4406304/1‑93 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                                 <  1:3417397/93‑1 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                                 <  2:3128863/93‑1 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                                 >  1:317445/1‑93 (MQ=255)

                               ggCCATGCGCTGGAAGCGGTGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                                 <  2:1872397/93‑1 (MQ=255)

                                                          aa‑‑‑‑aCT‑GGTCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACTGGATGACCGACGTCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATg                      <  2:5137555/92‑1 (MQ=255)

                                                                actggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg                    >  1:2061326/5‑93 (MQ=255)

                                                                 ctggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   >  2:3258484/4‑93 (MQ=255)

                                                                 ctggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                   <  2:1841940/90‑1 (MQ=255)

                                                                  tggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCa                  <  1:4044542/91‑1 (MQ=255)

                                                                   ggg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                 >  2:2165022/2‑93 (MQ=255)

                                                                    gg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                <  1:9065/93‑1 (MQ=255)

                                                                    gg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                <  2:4212740/93‑1 (MQ=255)

                                                                    gg‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                <  2:2856839/93‑1 (MQ=255)

                                                                     g‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg               <  2:3894347/93‑1 (MQ=255)

                                                                       cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc              <  1:1314607/93‑2 (MQ=255)

                                                                       cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc              <  2:6266855/93‑2 (MQ=255)

                                                                        tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             <  2:4395037/93‑3 (MQ=255)

                                                                        tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             <  1:4886730/93‑3 (MQ=255)

                                                                        tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             >  1:2762028/1‑91 (MQ=255)

                                                                        tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             >  1:1621121/1‑91 (MQ=255)

                                                                        tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg             >  1:2160175/1‑91 (MQ=255)

                                                                            gcgtGAAGGATTGTATTGCGTTGGGTGACGGAATGAACGACGCCGAAATGCTGTCAATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         <  1:1656708/90‑1 (MQ=255)

                                                                            gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  1:5509161/3‑93 (MQ=255)

                                                                            gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc         >  1:2591349/3‑93 (MQ=255)

                                                                             ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca        >  1:1682782/2‑93 (MQ=255)

                                                                             ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacca        <  2:5886704/92‑1 (MQ=255)

                                                                              cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcaccac       >  2:1370507/1‑92 (MQ=255)

                                                                              cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg       <  1:2215785/93‑1 (MQ=255)

                                                                               tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc      <  2:443092/93‑1 (MQ=255)

                                                                               tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGc      >  1:2156490/1‑93 (MQ=255)

                                                                                gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg     <  2:5212621/93‑1 (MQ=255)

                                                                                gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCg     >  2:2799496/1‑93 (MQ=255)

                                                                                  aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  1:4746228/93‑1 (MQ=255)

                                                                                  aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:6133189/93‑1 (MQ=255)

                                                                                  aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:2312465/93‑1 (MQ=255)

                                                                                  aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:4358748/93‑1 (MQ=255)

                                                                                  aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:1387586/93‑1 (MQ=255)

                                                                                   ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  1:1624860/1‑93 (MQ=255)

                                                                                   ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  <  1:5708665/93‑1 (MQ=255)

                                                                                   ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  <  2:6165831/93‑1 (MQ=255)

                                                                                   ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  <  2:1178316/93‑1 (MQ=255)

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503

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