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.4%	544.4 / 29.2	167	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (4/7); new base C (71/84); total (76/91)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                      <  1:3536123/100‑1 (MQ=255)

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                      >  1:1673696/1‑100 (MQ=255)

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

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

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCAATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                   <  1:757635/99‑1 (MQ=255)

   acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑c                                                                                               <  1:2516283/100‑1 (MQ=255)

   acaGCCTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑c                                                                                               <  1:3487891/100‑1 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑c                                                                                            <  1:2378075/100‑1 (MQ=255)

      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑cc                                                                                           <  1:2850467/100‑1 (MQ=255)

      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑cc                                                                                           >  2:601676/1‑100 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑cccc                                                                                         >  1:393317/1‑99 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTCGATTTCTAAAGTCGGTGGATTTTCTCGTTT‑‑c‑c‑gc‑cccc                                                                                         >  2:3595829/1‑99 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑ccccg                                                                                        <  1:3705979/100‑3 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑ccccgcc                                                                                      <  1:1240219/100‑5 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑ccccgcca                                                                                     >  1:1951365/1‑95 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc‑ccccgcc                                                                                      >  2:3700944/1‑95 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAgc                                                                                   >  1:2216734/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAgc                                                                                   <  2:2823015/100‑1 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGCCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAgcg                                                                                  <  1:952118/99‑1 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAgcgc                                                                                 <  1:366687/100‑1 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAgcgc                                                                                 >  1:1274799/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑c‑cc                                                                                              >  2:2598271/1‑88 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGccca                                                                              >  2:1716636/1‑96 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAAcc                                                                            <  2:250340/100‑1 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAAcc                                                                            >  2:2734010/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑C‑CC‑CCC‑GCCAGCGCCAACCAt                                                                          <  2:1889387/100‑1 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑C‑CC‑CCC‑GCCAGCGCCAACCAt                                                                          <  2:13294/100‑1 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCCACCATAg                                                                        >  2:231980/1‑100 (MQ=255)

                          aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGt                                                                       <  2:810139/100‑1 (MQ=255)

                            cGTCTCGTTCTCTACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGtt                                                                      <  1:497104/99‑1 (MQ=255)

                             gTCTGGTTCGCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTa                                                                    >  2:2620300/1‑100 (MQ=255)

                              tCTGGTTCTCGCCGCCCATTGTTGCTAATTGCCGATAGTTGATCCCTACCGTCCGTGCATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTAt                                                                   <  2:1821861/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTAt                                                                   >  2:124266/1‑100 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTAt                                                                   <  2:889913/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTAt                                                                   <  2:1287308/100‑1 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATc                                                                  <  2:294230/100‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATcc                                                                 >  1:2740603/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATcc                                                                 <  1:664241/100‑1 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATcc                                                                 >  1:2834000/1‑100 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCa                                                                >  1:156285/1‑100 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCa                                                                >  1:1920803/1‑100 (MQ=255)

                                  gTTGTCGACGCTCATTGTTGCTAATCGAGGATGGTTGATTTCCAATGTCAGTGAATTTTCTCGTTT‑‑C‑CTCC‑TCT‑GCCAGCGGCAACAATAGTTTTTcc                                                                 <  2:3353919/98‑1 (MQ=255)

                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCAt                                                               >  1:1387330/1‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2794850/100‑1 (MQ=255)

                                         aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCCGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:415974/96‑1 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:46488/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:1737514/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:1429670/1‑100 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:335833/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCTGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:1837866/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:3049763/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:1487784/100‑1 (MQ=255)

                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  1:567239/1‑100 (MQ=255)

                                                     gCTAATGGCCGATTGTTCATCCCTAACGTCGGTAAATTTTCCCGTTA‑‑C‑C‑CCACCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  2:766722/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  1:88014/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  2:2818666/100‑1 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCACCCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  2:1012247/1‑99 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  1:3523398/1‑100 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  2:163788/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:2548434/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:2714865/100‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTGTCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  2:928954/1‑100 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  1:3757053/100‑1 (MQ=255)

                                                           ggCGGATGTTTGATTTCTACCGCCGGGGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  1:3412520/100‑1 (MQ=255)

                                                           ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  1:2959552/100‑1 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:1951354/100‑1 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGGGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  1:3634209/1‑100 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:1810472/99‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:1925057/100‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:3210718/100‑1 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  2:2957782/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:1462614/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGGGAAGTGGCGAAAg                              >  1:1319248/1‑100 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:3198735/100‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  1:332724/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          >  2:112625/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  1:2734030/100‑1 (MQ=255)

                                                                          tCTAATGTCGGTGAATTTTCTCGTTC‑‑C‑C‑CC‑CCCTTTCTGCGCTAACCATAGTTTATCCATTTGCGCCCGTGTTTCATCGTGGAGTGGCGAAAGCCACAGa                       <  2:322035/100‑1 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTCCC‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCGTCGTGAAGTGGCGAAAGCCACAg                        >  2:1670663/1‑100 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  1:103319/1‑100 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:357299/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:2254324/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTCC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:3320622/100‑1 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:800041/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:2984583/100‑1 (MQ=255)

                                                                                     tGAATTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  1:2697801/1‑100 (MQ=255)

                                                                                       aaTTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          <  1:2031789/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGGCtgatg            >  2:1322320/1‑95 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:3111752/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:265393/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:305237/100‑1 (MQ=255)

                                                                                         ttttCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGGGGCGAAAGCCACAGATGCGCCtgag             <  2:3551236/95‑1 (MQ=255)

                                                                                          tttCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:1141912/100‑1 (MQ=255)

                                                                                          tttCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGCGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       >  2:705146/1‑100 (MQ=255)

                                                                                           ttCTCGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgaattagggt      >  2:2050394/1‑92 (MQ=255)

                                                                                           ttCTCGTTT‑‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  2:1776395/1‑100 (MQ=255)

                                                                                             ctcGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTTGCGAAAGCCACAGCTGCGCCTGAGTGAGTGTa     >  2:1102025/1‑99 (MQ=255)

                                                                                             ctcGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCAAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:2860274/1‑100 (MQ=255)

                                                                                             ctcGTTT‑‑CGC‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  2:1420508/100‑1 (MQ=255)

                                                                                              tcGTTTC‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:328117/100‑1 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑C‑CC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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