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	94.3%	277.6 / 12.3	88	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (5/0); new base C (38/45); total (43/45)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACG  >  NC_000913/3379160‑3379352

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

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

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTATTGCTAATGGCAGATGGATGATTTCTAACATCGGTGAATTTTCTCGttt                                                                                                       >  1:270175/1‑100 (MQ=255)

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

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

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

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

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑c                                                                                             <  2:1166890/100‑1 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑ccg                                                                                      <  2:188802/100‑3 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTATCGTCGGTGAATTTTCTCTTTT‑‑‑c‑‑cc‑cc‑c‑cctg                                                                                     <  1:1785425/100‑4 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccg                                                                                     >  2:302446/1‑97 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccg                                                                                     <  2:611310/100‑4 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑ccgcc                                                                                    >  1:789785/1‑96 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccgcc                                                                                   >  2:2621432/1‑95 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTAAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccgc                                                                                    >  2:778912/1‑95 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTAGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑ccgcca                                                                                   >  1:2522745/1‑94 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAgc                                                                                  <  1:1736563/100‑1 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccgccag                                                                                 >  2:1989919/1‑93 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑CCACC‑C‑CGCCAgc                                                                                  >  1:1867613/1‑99 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAgcgc                                                                                >  1:1904589/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAgcgc                                                                                <  1:2571119/100‑1 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑cc‑cc‑c‑cccgccagcg                                                                               <  2:1431469/100‑10 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGcc                                                                               <  2:1818119/100‑1 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCTCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCa                                                                              <  1:469747/100‑1 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAAc                                                                            >  2:564917/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAAcc                                                                           <  1:1680327/100‑1 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAAcc                                                                           >  1:2588576/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAAcc                                                                           <  1:1073531/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGt                                                                      >  1:1413141/1‑100 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGtt                                                                     <  2:2650716/100‑1 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGtt                                                                     <  1:1928032/100‑1 (MQ=255)

                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTa                                                                   <  1:88251/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTAt                                                                  >  1:505885/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCCTTGTTGCTTATGGCCGATCGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCACAGtt                                                                     <  2:2636583/97‑1 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTAt                                                                  >  2:1232288/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTAt                                                                  >  2:108700/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATc                                                                 <  1:2355484/100‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAATTTAt                                                                  >  2:1140122/1‑99 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATcc                                                                >  1:2243356/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATcc                                                                <  2:1713475/100‑1 (MQ=255)

                                  gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCa                                                               <  1:954277/100‑1 (MQ=255)

                                   ttCTCGCCGAGCATTGTTGCTAATGGCCCATGATTGGTTTCTAACGTCTGTGGATTTTCTCGTTT‑‑‑C‑‑CC‑CC‑CATGACAGCGCCACCAATACCTTTTCCAt                                                              >  2:1846045/1‑99 (MQ=255)

                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CACCAGCGCCAACCATAGTTTATCCAtt                                                             >  1:2520173/1‑100 (MQ=255)

                                     ctcGAAGCGCATTGTTGCTAATGGCCGATTGTTGATTTCTAACGTCGGTAAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCAttt                                                            <  2:1155726/99‑1 (MQ=255)

                                      ttgatgCCCTTTTTTTCTTATTGCCGATGGTTGATTTCTTACGTCGGTGAATTTTCTCGGTTCCCC‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCAtt                                                             <  1:2192074/94‑1 (MQ=255)

                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTAGTTT‑‑‑CCTCC‑CC‑C‑CGCCAGCAAAAAAAATAGTTTATACAttt                                                            >  1:1296897/1‑99 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGCCCGATGGTTGATTTCTAACCTCGGTGAATGTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCCACCACAGTTTATCCATCTg                                                           >  2:281085/1‑98 (MQ=255)

                                           gCCCATTGTTGCTAATGGCCGATGTTTGATTTCTAACGTCGGTTAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                     <  1:704360/100‑1 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                      >  2:2106320/1‑99 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                     >  2:742827/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                    >  1:676035/1‑100 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                   <  1:1305135/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATACATTTGCGCCCggg                                                   >  1:652409/1‑100 (MQ=255)

                                              cATTGTTGTTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                   <  1:1311218/99‑1 (MQ=255)

                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                   <  2:1033126/100‑1 (MQ=255)

                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑CC‑CCAC‑CGCCAGCACCAACCATAGTTTATCCATTTTCGCCCggg                                                   >  1:2212811/1‑99 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTATCGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                   <  1:115653/99‑1 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTTAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                 <  1:844191/100‑1 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGTGCCAACCATAGTTTATCTATTTGCGCCCgtg                                                   <  1:1785796/99‑3 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                 <  2:1372454/100‑1 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCCTAGTTTATCCATTTGCGCCCGGGTTTc                                               >  2:15906/1‑100 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                               >  2:1017805/1‑100 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                               >  2:394755/1‑100 (MQ=255)

                                                     gccaaTGGCCGATGGTTGATTTCTAATGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑GC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                           >  1:1833183/4‑99 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                            <  1:1938433/100‑1 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                          <  1:2052932/100‑1 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                         >  1:113724/1‑100 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                        <  2:1310057/100‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                       <  1:2332196/100‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                      <  2:238434/100‑1 (MQ=255)

                                                           ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCAAAATTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                       >  1:1555708/1‑99 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                   >  1:2170192/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                   <  1:1823317/100‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCGTTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                 <  1:1922289/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                 <  2:648037/100‑1 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CCCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                               >  1:1631212/1‑100 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                              >  2:395848/1‑100 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                             <  2:505881/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                             <  1:395850/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTATCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAAACCTAGTTTTTCCACTTGCGCCCGGGTTTCATCCTTAAGCGGCGaaa                              >  2:801122/1‑99 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                           >  1:1383643/1‑100 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                          >  1:2653408/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                          >  1:2228160/1‑100 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                         <  2:1573035/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCA‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACAATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                         >  2:842526/1‑100 (MQ=255)

                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTCCCC‑‑CC‑CC‑C‑CTCCAACGCCAACCCTAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                         >  2:1776092/1‑100 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTTTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                      <  1:606357/100‑1 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                      >  2:1508741/1‑100 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCCGCGCCAACCATAGTTTATCCATTTGCGCCCAGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                     >  1:414441/1‑99 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                    >  2:1898082/1‑100 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                  <  2:2243329/100‑1 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCTAGCGTCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                  <  1:384940/100‑1 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGAGCCCGGGTTTCATCGTGAAGTGGCGAAAACCACAGATgcg                  >  1:403621/1‑100 (MQ=255)

                                                                                   ggTGAATTTTCTCGTTTCC‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAAAATGAAGTGGCGAAAGCCACAGATGCGCCt               >  1:57182/1‑99 (MQ=255)

                                                                                    gTGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag            <  1:2264861/100‑1 (MQ=255)

                                                                                     tGAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt           <  1:2290419/100‑1 (MQ=255)

                                                                                      gAATTTTCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg          <  1:2097143/100‑1 (MQ=255)

                                                                                           ttCTCGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt     >  1:1240356/1‑100 (MQ=255)

                                                                                            tctcGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa    <  2:789777/100‑1 (MQ=255)

                                                                                               cGTTTC‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCTTTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg  <  2:2169228/99‑1 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑CC‑CC‑C‑CGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACG  >  NC_000913/3379160‑3379352

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