BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,266:1	(C)_7→8	100%	coding (732/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,259	1	.	C	93.0%	371.4 / 19.8	114	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (2/6); new base C (60/46); total (62/52)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379162‑3379353

cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c                                                                                                  <  2:2859755/100‑1 (MQ=255)

cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c                                                                                                  <  2:2859731/100‑1 (MQ=255)

cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c                                                                                                  >  2:257806/1‑100 (MQ=255)

 acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCTCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cc                                                                                                 <  2:1200935/100‑1 (MQ=255)

 acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cc                                                                                                 >  2:799155/1‑100 (MQ=255)

 acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cc                                                                                                 <  2:2908776/100‑1 (MQ=255)

  caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc                                                                                                <  1:2808362/100‑1 (MQ=255)

  caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc                                                                                                <  1:886237/100‑1 (MQ=255)

    gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTC‑‑c‑ccc‑‑c‑c                                                                                           >  2:1110964/1‑100 (MQ=255)

      gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc‑‑c‑c‑cc                                                                                        >  2:2107968/1‑99 (MQ=255)

        agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc‑‑c‑c‑ccgc                                                                                      >  1:1011956/1‑97 (MQ=255)

         gtaGAAGAGACCGCCAGCGTCTGGTTCTCGACGCCTATTGTTGCTAATGGCCGAGGGTTGATTTCTAACGTCGGTGAATTTTCTCGttttgc                                                                                                    >  1:1749100/1‑89 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc‑‑c‑c‑ccgcc                                                                                     >  1:155952/1‑96 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc‑‑c‑c‑ccgcc                                                                                     >  2:1931156/1‑96 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccc‑‑c‑c‑ccgcc                                                                                     <  2:643174/100‑5 (MQ=255)

           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAg                                                                                    >  1:1245332/1‑100 (MQ=255)

           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAg                                                                                    >  1:2809180/1‑100 (MQ=255)

           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAg                                                                                    >  1:2315320/1‑100 (MQ=255)

            aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAgc                                                                                   >  2:589518/1‑100 (MQ=255)

             aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAgcg                                                                                  >  1:1788216/1‑100 (MQ=255)

             aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAg                                                                                    >  1:414054/1‑98 (MQ=255)

                 gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGCTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCa                                                                               >  1:1057177/1‑99 (MQ=255)

                  aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCC‑‑C‑C‑CGCCAGCGccca                                                                              >  2:31162/1‑98 (MQ=255)

                  aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAAc                                                                             <  1:2230405/100‑1 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCC‑‑C‑C‑CGCCAGCGCCAAc                                                                             >  2:1732976/1‑100 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAAcc                                                                            >  2:2110326/1‑100 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAAcc                                                                            >  1:2512459/1‑100 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAAcc                                                                            <  2:468765/100‑1 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CTCCC‑‑C‑C‑CGCCAGCGCCAAcc                                                                            <  1:2499002/100‑1 (MQ=255)

                      ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCC‑‑C‑C‑CGCCAGCGCCAACCAt                                                                          >  2:1540729/1‑100 (MQ=255)

                      ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATa                                                                         >  2:768894/1‑100 (MQ=255)

                      ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATa                                                                         <  2:2835284/100‑1 (MQ=255)

                       cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAg                                                                        <  1:1680127/100‑1 (MQ=255)

                       cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAg                                                                        >  1:309568/1‑100 (MQ=255)

                        aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGt                                                                       >  2:2174069/1‑100 (MQ=255)

                             cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATc                                                                  <  2:431143/100‑1 (MQ=255)

                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCa                                                                >  2:2141266/1‑100 (MQ=255)

                                   ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTg                                                            >  2:1763402/1‑100 (MQ=255)

                                      gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  1:1213809/1‑100 (MQ=255)

                                        cGCCCATTCTTGTTAAAGGCCGATGGTTGATTCCTCACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCTC‑C‑CGCCAGCACCAACCATAGTTTATCCATTTgcg                                                          <  1:2741515/98‑1 (MQ=255)

                                         gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      <  2:79865/100‑1 (MQ=255)

                                           cgaTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:2468643/98‑1 (MQ=255)

                                           ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:2132986/100‑1 (MQ=255)

                                            cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   <  2:2011488/100‑1 (MQ=255)

                                            cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCACCCGGGt                                                   >  1:1562843/1‑100 (MQ=255)

                                            cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGAGCCAACCGTAGTTTATCCATTTGCGCCCggg                                                    >  1:181482/1‑99 (MQ=255)

                                              ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  2:2466466/100‑1 (MQ=255)

                                               tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  2:631600/1‑100 (MQ=255)

                                                gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCAAAGTTTATCGATTTGCGCCCGGGGTTc                                                >  1:1366749/1‑99 (MQ=255)

                                                  tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:2550721/1‑100 (MQ=255)

                                                  cctgcacctccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  1:2437576/91‑1 (MQ=255)

                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:2645962/100‑1 (MQ=255)

                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  1:1258395/1‑100 (MQ=255)

                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:426576/100‑1 (MQ=255)

                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGGTTCATc                                             >  1:692138/1‑99 (MQ=255)

                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  1:17753/1‑100 (MQ=255)

                                                     tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:699452/1‑100 (MQ=255)

                                                     tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  2:768879/100‑1 (MQ=255)

                                                     tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:78258/1‑100 (MQ=255)

                                                       aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                        <  1:1468424/100‑1 (MQ=255)

                                                         ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  1:69744/1‑100 (MQ=255)

                                                         ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  2:2093644/100‑1 (MQ=255)

                                                            cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATAGTGAAGTGGc                                   >  1:477837/1‑100 (MQ=255)

                                                                 gctGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  1:477259/3‑99 (MQ=255)

                                                                 gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:54996/1‑100 (MQ=255)

                                                                 gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:50259/100‑1 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGCATTTTCTCGTTTCCC‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:1375037/100‑1 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCCACCATAGTTTATCCATTTGCGCCCGGGTTTCATCTTGAAGTGGCGAAAGc                             >  2:877997/1‑100 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  2:1940753/100‑1 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  1:855692/1‑100 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  2:229876/1‑100 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  2:875015/1‑100 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  1:1864016/1‑100 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTG‑‑C‑CCC‑‑C‑CGCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:2046058/99‑1 (MQ=255)

                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCCTAGTTTATCCATTTGTGCCCGGGTTTCATCGTGAAGTGGCGCAAGcc                            >  1:2913250/1‑100 (MQ=255)

                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:2518650/1‑100 (MQ=255)

                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:1529408/100‑1 (MQ=255)

                                                                     aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGCAGTGGCGAAAGCcac                          <  1:705667/100‑1 (MQ=255)

                                                                     aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  1:99765/100‑1 (MQ=255)

                                                                      tttCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCCTTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         >  2:1445936/1‑100 (MQ=255)

                                                                        tCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCC‑‑CGC‑CGCCAGCGCCACCCATAGGTTATAGATTTGCGCCCGGATTTCATCCTGAAGTTGCGACAGCCACAg                        >  1:1065063/1‑99 (MQ=255)

                                                                          tAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:88587/100‑1 (MQ=255)

                                                                           aaCGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATCGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  1:259583/1‑98 (MQ=255)

                                                                           aaCGTCGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:888452/1‑100 (MQ=255)

                                                                               tCGGTGAATTTTCTCGTTTCCC‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:182213/100‑1 (MQ=255)

                                                                                cGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGccgg               >  2:1500667/1‑98 (MQ=255)

                                                                                cGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:2763015/100‑1 (MQ=255)

                                                                                cGGTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:575070/100‑1 (MQ=255)

                                                                                 ggTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  1:2479489/1‑99 (MQ=255)

                                                                                  gTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCACCCATAGTTTATCCATTTCCGCCCGGGTTTCACCGTCAATTGGACAAAGCCACAGAtacccc                 >  1:2785689/1‑91 (MQ=255)

                                                                                  gTGAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:1252643/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:1104868/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:2794140/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:2550706/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTAGCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:175661/100‑1 (MQ=255)

                                                                                      aTTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:675127/100‑1 (MQ=255)

                                                                                      aTTTTCTCGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGAGGCGAAAGCCACCGATGTGACtgagtta          >  1:2828785/1‑97 (MQ=255)

                                                                                           ctcGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:1619177/1‑100 (MQ=255)

                                                                                            tcGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:2793347/100‑1 (MQ=255)

                                                                                             cGTTTC‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:1562849/100‑1 (MQ=255)

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCC‑‑C‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379162‑3379353

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