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.5%	560.9 / 29.4	168	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (5/6); new base C (73/84); total (78/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

gCACAGCGTAGGAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGAGGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                  >  1:199061/1‑100 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                 <  2:2751531/100‑1 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                 <  1:355140/100‑1 (MQ=255)

  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccc                                                                                                >  2:1715005/1‑100 (MQ=255)

  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccc                                                                                                <  1:973332/100‑1 (MQ=255)

    atcGTAGTAAAAGATACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTATAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c                                                                                             >  2:364317/3‑100 (MQ=255)

         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCTTTT‑‑gtcc‑c‑cccc                                                                                        >  2:940698/1‑97 (MQ=255)

          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c‑cccgcc                                                                                      >  1:3411050/1‑96 (MQ=255)

          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c‑cccgcc                                                                                      <  1:862891/100‑5 (MQ=255)

          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c‑cccgcc                                                                                      >  1:1950253/1‑96 (MQ=255)

           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c‑cccgcca                                                                                     >  1:231372/1‑95 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAg                                                                                    >  2:2232258/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAg                                                                                    >  2:1845366/1‑100 (MQ=255)

             aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAgc                                                                                   >  1:4343524/1‑100 (MQ=255)

                gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGcc                                                                                >  2:1348043/1‑100 (MQ=255)

                   aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAAc                                                                             <  1:3107052/100‑1 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAAc                                                                             <  1:601506/100‑1 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCAt                                                                          <  2:541631/100‑1 (MQ=255)

                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATa                                                                         <  2:1859452/100‑1 (MQ=255)

                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAg                                                                        <  1:4252566/100‑1 (MQ=255)

                         aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAg                                                                        <  1:3152636/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTa                                                                    <  1:4360991/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTa                                                                    <  2:2179020/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTACTAATGGCCGATGATTGATTTCTAACGTCCGAGAATTTTCTAGTTT‑‑CCCC‑C‑CCAGCCAGCGCCAACCATAGTTTa                                                                    <  1:364316/100‑1 (MQ=255)

                             tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGttt                                                                     >  2:162449/1‑98 (MQ=255)

                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATc                                                                  >  2:1956815/1‑100 (MQ=255)

                               tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTAt                                                                   >  1:4119304/1‑98 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATcc                                                                 >  1:260872/1‑100 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGAAGGTTGATTTCTAACGTCGGTGAATTTTTTCGTTTC‑CCCC‑C‑CC‑GCCAGCCCCAACCATAGTTTATCCa                                                                >  1:3413864/1‑100 (MQ=255)

                                 gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCa                                                                <  2:441856/100‑1 (MQ=255)

                                  ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCAt                                                               >  2:4205320/1‑100 (MQ=255)

                                  ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCAtt                                                              >  1:325931/1‑100 (MQ=255)

                                  ttCTCGACGCCCATTCTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCAtt                                                              <  2:4343556/100‑1 (MQ=255)

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

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

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

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCAttt                                                             >  1:3688501/1‑99 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:1821494/100‑1 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:251008/99‑1 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:417281/99‑1 (MQ=255)

                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:4246330/100‑1 (MQ=255)

                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCAAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  1:2364721/1‑100 (MQ=255)

                                       ggcGCCCAGTGTGGCTAATGGCCGAGGGTTGATTTCTAACGTCGGGGAATTTTCGCGGTT‑‑CCCC‑CGCC‑GCCAGCGCCAACCATATTTTATCCATTTg                                                            <  1:923977/95‑1 (MQ=255)

                                       gccGCCCATTGTTACTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  2:2993442/3‑100 (MQ=255)

                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑CCCAGGGCCAACCATAGTTTATCCttt                                                              >  2:2001614/1‑96 (MQ=255)

                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  2:3903249/1‑100 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGTCAACCATAGTTTATCCATTTTCGCCCg                                                      >  1:343798/1‑100 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      <  1:1072603/100‑1 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTGC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      <  1:2881647/100‑1 (MQ=255)

                                           cccTTTGTGGCTAATGGCGGAGGGTTCATTCCGAACCTGGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:645950/96‑1 (MQ=255)

                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:1374431/1‑100 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑TCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  2:3116514/100‑1 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  2:3942074/100‑1 (MQ=255)

                                              aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:29019/100‑1 (MQ=255)

                                              aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  >  1:3317379/1‑100 (MQ=255)

                                                 gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               >  1:3085113/1‑100 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                              >  2:371035/1‑100 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:244261/1‑100 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑CCCAGCGCCAACCATAGTTTATCCCTTTTCGCCCGGGTTTCCTc                                             >  2:3060824/1‑100 (MQ=255)

                                                    gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:553015/1‑100 (MQ=255)

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

                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTTAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:2005311/99‑1 (MQ=255)

                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  1:1219979/100‑1 (MQ=255)

                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  1:1275154/100‑1 (MQ=255)

                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  1:1596899/100‑1 (MQ=255)

                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGTGCCCGGGTTTCATCGTGa                                         >  2:3865149/1‑100 (MQ=255)

                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATGTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:1190671/100‑1 (MQ=255)

                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  2:954003/1‑100 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:2496592/1‑100 (MQ=255)

                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  2:2865321/1‑100 (MQ=255)

                                                             cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:1867704/100‑1 (MQ=255)

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

                                                               aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                 >  2:258730/1‑100 (MQ=255)

                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCGGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:2794360/1‑100 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATCGTTTATCCACTTGCGACCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          >  1:2979274/1‑100 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCTC‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  1:3946862/100‑1 (MQ=255)

                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:2772087/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:4025360/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:1374426/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  2:2916080/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  2:556417/100‑1 (MQ=255)

                                                                         tCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCAACAGa                       >  2:1115540/1‑100 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  2:221364/100‑1 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:2672366/1‑100 (MQ=255)

                                                                            aaCGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCTTCGTGAAGTTGCGAAAGCCACAGATgc                    >  2:59181/1‑100 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCTTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:2106101/100‑1 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:1152732/100‑1 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:4193817/100‑1 (MQ=255)

                                                                              cGTCGGTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  <  2:1243125/100‑1 (MQ=255)

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

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

                                                                                tCGGTGAATTTTCTCGTTTCCCCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:1827506/100‑1 (MQ=255)

                                                                                  ggTGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGaagatcggaa              >  2:2057752/1‑91 (MQ=255)

                                                                                  ggTGAATTTTCTCGTTCC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              <  1:3226821/100‑1 (MQ=255)

                                                                                    tGAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  1:1115810/100‑1 (MQ=255)

                                                                                     gAATTTTCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtn           <  2:1225959/100‑2 (MQ=255)

                                                                                          ttCTCGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:2731481/1‑100 (MQ=255)

                                                                                           tctcGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:3429457/1‑100 (MQ=255)

                                                                                            ctcGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTaa    >  1:3847810/1‑99 (MQ=255)

                                                                                              cGTTTC‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGGTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:851606/100‑1 (MQ=255)

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCC‑C‑CC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  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: `‑`