BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)_7→8	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.5%	546.1 / 16.6	154	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (5/2); new base C (79/68); total (84/70)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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

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

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

 gCACATCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                  <  1:1703190/100‑1 (MQ=255)

 gCACATCGTAGTAAAAGAGACCGCCAGCGGCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                     <  2:3783265/99‑1 (MQ=255)

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

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                 >  2:614139/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cc                                                                                              >  2:3306235/1‑100 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccc                                                                                             <  2:2175546/100‑1 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccc                                                                                             <  2:1401504/100‑1 (MQ=255)

       cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccccc                                                                                           <  2:2230634/100‑1 (MQ=255)

       cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccccc                                                                                           >  2:4058702/1‑100 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cccccc                                                                                          <  1:1749555/100‑2 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cccccc                                                                                          >  1:3416088/1‑99 (MQ=255)

         tagtaAAGGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cccccc                                                                                          <  2:3358281/99‑2 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccccccgc                                                                                        >  1:584568/1‑97 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑cccccc                                                                                          >  2:628583/1‑97 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccccc‑‑ccc                                                                                      >  1:580474/1‑100 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc‑ccccccgcc                                                                                       >  2:3150454/1‑96 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAg                                                                                    >  2:751644/1‑100 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAgcg                                                                                  <  1:255865/100‑1 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAgcg                                                                                  <  2:3731884/100‑1 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAgcg                                                                                  >  1:19832/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAgcgc                                                                                 >  2:2627322/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAgcgc                                                                                 >  2:3255494/1‑100 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGcc                                                                                >  1:1521669/1‑100 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGcc                                                                                <  1:3468477/100‑1 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCa                                                                               <  2:2662105/100‑1 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCa                                                                               <  1:1086716/100‑1 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAAc                                                                             <  1:2628054/100‑1 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAAcc                                                                            >  2:2676797/1‑99 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAAcc                                                                            <  2:903486/99‑1 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCa                                                                           >  1:2754117/1‑100 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCAt                                                                          >  1:1651525/1‑100 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTGTAACGTCGGTGAATATTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATa                                                                         >  1:103765/1‑100 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATa                                                                         >  2:1945837/1‑100 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAg                                                                        <  1:2429088/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGtt                                                                      <  1:943136/100‑1 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGttt                                                                     >  2:1799846/1‑100 (MQ=255)

                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTa                                                                    >  1:2223171/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATcc                                                                 >  1:1693578/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATcc                                                                 <  2:1998896/100‑1 (MQ=255)

                                  gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCAt                                                               >  1:2796105/1‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCAttt                                                             >  2:3089108/1‑99 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:3163678/100‑1 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:678755/1‑100 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  2:2819621/1‑100 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  1:2299185/1‑100 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:2055525/1‑100 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGAATTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:2745889/1‑100 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  1:191961/1‑100 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:2372405/100‑1 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  2:167740/100‑1 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  2:3832154/1‑100 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCCCCCACCATAGTTTATTCATTTTCGCCCAGGttt                                                 >  2:575196/1‑98 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                              <  2:1483333/100‑1 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  2:1450144/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  2:3167196/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  2:3912867/100‑1 (MQ=255)

                                                      cTAATTGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  2:2516917/99‑1 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  2:2462030/1‑100 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTATAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTTATCGt                                           >  2:476806/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  1:4132047/1‑100 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  1:329724/100‑1 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTTCGCCCGGGTTTCATCGTGaa                                        >  2:1782145/1‑97 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:3832197/100‑1 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  1:4075641/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGGCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:816541/100‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:3083502/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                 <  2:496720/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                 >  2:1201489/1‑100 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  2:2249722/1‑100 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  1:3267538/1‑100 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:3067008/1‑100 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  1:2900874/100‑1 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:3076463/1‑100 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:479903/100‑1 (MQ=255)

                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  2:2299165/100‑1 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  1:3204533/1‑100 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  1:2162530/1‑100 (MQ=255)

                                                                           ntAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  2:4105630/2‑98 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  1:3255582/1‑100 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  2:1650197/1‑100 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  1:2053000/1‑100 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  1:1201498/100‑1 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:3592456/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTT‑‑CCACCCCC‑‑GCCAGCGCCAACCATAGCTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:1412220/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTT‑‑CC‑CCCCCCGGCCAGCGCCAACAATGGTTTATCTATTTGCGCCCCGGTTTAATCGTGAATTGGACAAAGCCACCGATgcgc                  >  2:3583280/1‑100 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  2:1071982/1‑100 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCAATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:3416050/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  1:2764157/1‑100 (MQ=255)

                                                                                     tGAATTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  2:3611337/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgtg         <  1:2819653/100‑3 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:431774/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:509901/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:2762150/100‑1 (MQ=255)

                                                                                           ttCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:517176/1‑100 (MQ=255)

                                                                                           ttCTCGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:3205000/1‑100 (MQ=255)

                                                                                            tctcGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  1:26530/100‑1 (MQ=255)

                                                                                              tcGTTTC‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   >  2:2861214/1‑100 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CC‑CCCCC‑‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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