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	95.8%	331.8 / 5.6	96	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (41/51); total (42/54)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                 <  2:1926059/100‑1 (MQ=255)

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

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATTGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cccc‑cctg                                                                                       <  1:1613713/100‑3 (MQ=255)

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cccc‑cccg                                                                                       <  1:2094356/100‑3 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAg                                                                                    >  2:611631/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cccc‑cccgcc                                                                                     >  1:1491075/1‑94 (MQ=255)

              aaGAGACCGCCAGCGTCTGGTTCTCGACGTCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAgcg                                                                                  <  2:2193143/100‑1 (MQ=255)

              aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAgcg                                                                                  <  1:1742857/100‑1 (MQ=255)

               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAgcgc                                                                                 <  1:998768/100‑1 (MQ=255)

                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCCGCGCCaa                                                                              >  1:199817/1‑100 (MQ=255)

                  gaCCGCCAGCGTATGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCaa                                                                              >  1:1565792/1‑100 (MQ=255)

                   aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAAc                                                                             <  2:892704/100‑1 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCTCGTTT‑‑c‑cccc‑ccccccagcggc                                                                               >  2:1024832/1‑89 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAAcc                                                                            >  1:253541/1‑100 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAAcc                                                                            >  1:709087/1‑100 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCAt                                                                          >  2:1153899/1‑100 (MQ=255)

                           cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGttt                                                                     >  2:1935407/1‑100 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCTCCAGCGCCAACCATAGtt                                                                      >  2:855097/1‑98 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTa                                                                    >  1:2174641/1‑100 (MQ=255)

                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:2165293/100‑1 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCCC‑CCGCCAGCGCCAACCATAGTTTATcc                                                                 <  1:2305589/100‑1 (MQ=255)

                                 gTTCTCAACGCCCATTGTTTCTAATGGCCGATGGTTGATTTCCACCGTCGCTGAATTTTCTCGTGC‑‑CACCCC‑CCGCAAGCGCCCACCATAGCTTATc                                                                  <  2:1153211/97‑1 (MQ=255)

                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1843268/100‑1 (MQ=255)

                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1362762/100‑1 (MQ=255)

                                      cGACGCTCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:511300/100‑1 (MQ=255)

                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:1773955/1‑100 (MQ=255)

                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  1:1627755/1‑100 (MQ=255)

                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  2:2064421/1‑99 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:334870/1‑100 (MQ=255)

                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:74621/1‑100 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACAATAGTTTATCCATTTGCGCCCgg                                                     >  2:411496/1‑99 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:1010524/1‑100 (MQ=255)

                                                     cTAATCGCCCATTGTTGATCTCTCACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  2:2099994/98‑1 (MQ=255)

                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  1:1691488/1‑100 (MQ=255)

                                                       aaTGGCCGATGGTTGATTTCAAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  2:1455416/1‑100 (MQ=255)

                                                         tttccGTTGGTTGATTTCTAACTTCGTTGAATCTTCTCGTTTC‑C‑CCCC‑CCTCCAGCTCCAACCTTAGTTTATCCTTTTTCGCCCGGGTTTCATCGTg                                          <  2:609170/94‑1 (MQ=255)

                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:329566/100‑1 (MQ=255)

                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  1:475460/1‑100 (MQ=255)

                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  1:475707/1‑100 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  2:1546802/99‑1 (MQ=255)

                                                            ccGATGGTTGATTTCTAACGTCGGTGTATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCACCCGGGTTTCATCGTGAAgggg                                    <  2:135937/100‑1 (MQ=255)

                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  1:1596696/1‑100 (MQ=255)

                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:854616/100‑1 (MQ=255)

                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  2:1596684/100‑1 (MQ=255)

                                                             cgatagctcaTTTCTACCGTCCGTAAATTTTCTCTTTTC‑C‑CCCC‑CCGCCAGCGCCAGCCATAGTTTATCCGTTTGCGCCCGGGTTTCATCGCGAAGTAGc                                   <  2:1662983/91‑1 (MQ=255)

                                                             cGATGGTTGATTTCCAACGTCGGTGAAATTTCTCGTTTCCC‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:247929/100‑1 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:2356391/100‑1 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  >  2:728118/1‑100 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:509088/100‑1 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGGGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:269914/100‑1 (MQ=255)

                                                               aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:2313673/99‑1 (MQ=255)

                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCgca                                >  1:2314339/1‑98 (MQ=255)

                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                >  1:529417/1‑100 (MQ=255)

                                                                 ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:2185891/100‑1 (MQ=255)

                                                                 ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTAAAGTGGCGaaa                               >  1:43526/1‑100 (MQ=255)

                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:9102/1‑100 (MQ=255)

                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:1535424/100‑1 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:74622/100‑1 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:576830/100‑1 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCCTGAAGTGGCCaaa                               >  2:2193224/1‑98 (MQ=255)

                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  1:2197474/1‑100 (MQ=255)

                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTCCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  1:2267486/100‑1 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:406521/1‑100 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGTCGAAAGCcac                          <  1:113748/100‑1 (MQ=255)

                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:543828/100‑1 (MQ=255)

                                                                        ttCTAAGGTAGGTGAATTTTCTCGTTG‑‑C‑CCCCGCCGCCAGGGCCAACCATAGTTTATCCATGTGCGCGCGGGTTTCAGCGTGAAGTGGCGAAAGCCACAg                        <  2:826218/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTCTCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  2:393199/1‑100 (MQ=255)

                                                                        gatctaCGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  1:2137772/95‑1 (MQ=255)

                                                                          cTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  1:1948669/100‑1 (MQ=255)

                                                                          cTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  2:1413911/1‑100 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCACTTGCGCCCGGCTTTCATCGTGAAGTGGCTAAAGCCACAGacg                     <  1:1611157/100‑3 (MQ=255)

                                                                            aaCGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:2096836/100‑1 (MQ=255)

                                                                            aaCGTCGGTGAATTTTCTCGTTCC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:1637303/100‑1 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGagatc                   >  2:2137755/1‑96 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:497854/1‑99 (MQ=255)

                                                                             aCGTCGGTCAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:181575/100‑1 (MQ=255)

                                                                              cGTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGGGGCGAAAGCCACAGATgcg                   >  2:10020/1‑99 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:713497/100‑1 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:244944/100‑1 (MQ=255)

                                                                                 cGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:1153905/100‑1 (MQ=255)

                                                                                 cGGTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:1913068/100‑1 (MQ=255)

                                                                                  ggTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  1:469431/1‑100 (MQ=255)

                                                                                   gTGACTTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGTGTTTCTTCGTGAAGTGGCGAAAGCCACAGATGCGCCtgaa             <  2:2073056/100‑2 (MQ=255)

                                                                                   gTGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:928862/100‑1 (MQ=255)

                                                                                    tGAATTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  2:2150802/100‑1 (MQ=255)

                                                                                       aTTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:866828/100‑1 (MQ=255)

                                                                                       aTTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:1688630/100‑1 (MQ=255)

                                                                                       aTTTTCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:897046/100‑1 (MQ=255)

                                                                                        ttttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:47037/1‑100 (MQ=255)

                                                                                         tttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:995919/100‑1 (MQ=255)

                                                                                         tttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:475458/100‑1 (MQ=255)

                                                                                         tttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:475705/100‑1 (MQ=255)

                                                                                          ttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:1041707/1‑100 (MQ=255)

                                                                                          ttCTCGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:133973/1‑100 (MQ=255)

                                                                                           tctcGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  1:938821/100‑1 (MQ=255)

                                                                                             tcGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:533482/100‑1 (MQ=255)

                                                                                              cGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  1:2316262/100‑1 (MQ=255)

                                                                                              cGTTTC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:625226/1‑100 (MQ=255)

                                                                                              cGTTTC‑C‑CCCC‑CCGCCAGCGACAACCATAGTTTATCCATTTGCGCCCGGGTTTCCTCGTGAAGTGGCGAAAGCCAAAGATGCGCCTGAGTGAGTGTAc    >  2:633528/1‑98 (MQ=255)

                                                                                              cGTTT‑‑C‑CCCCTCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGTGATAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:1565781/100‑1 (MQ=255)

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑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: `‑`