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	96.0%	334.3 / 7.3	99	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (41/54); total (45/54)
Fisher's exact test for biased strand distribution p-value = 3.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 112 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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

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

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

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

 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                  >  2:1473401/1‑100 (MQ=255)

 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                  >  1:2032215/1‑100 (MQ=255)

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

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

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑gcc‑c‑c‑ccc                                                                                        >  1:2891326/1‑99 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccc‑c‑c‑cccg                                                                                       <  1:2001813/100‑3 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGGTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCC‑C‑C‑CCGc                                                                                       >  1:2678393/1‑100 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCCCGCTT‑‑ccc‑c‑c‑cccc                                                                                       >  1:2898240/1‑96 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccc‑c‑c‑cccgcc                                                                                     <  2:1002413/100‑5 (MQ=255)

            taAAAGAGACCGCCAGCGAGTGGTTCTCGACGCCCATTGTTGATAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGGTT‑‑ccc‑c‑c‑ccc                                                                                        >  1:557926/1‑95 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgc                                                                                   >  2:2477777/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgc                                                                                   >  2:2763061/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgc                                                                                   <  1:1919604/100‑1 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgcg                                                                                  <  1:2804627/100‑1 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgc                                                                                   >  1:2364163/1‑99 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAgcgc                                                                                 >  1:2046383/1‑100 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCa                                                                               >  1:975545/1‑99 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAAc                                                                             <  1:2681362/100‑1 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAAc                                                                             <  2:900470/100‑1 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAAc                                                                             <  2:2894775/100‑1 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAAcc                                                                            >  2:2187743/1‑100 (MQ=255)

                      caccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCa                                                                           >  2:2214471/3‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCa                                                                           <  1:2187937/100‑1 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCAt                                                                          >  1:1342641/1‑100 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAg                                                                        >  2:1260413/1‑100 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAg                                                                        <  1:819574/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGtt                                                                      <  2:2774681/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGtt                                                                      <  2:550201/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCAATTGTTGCTAATCGCCGCTGGTTGATTTCTAACGACGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTAt                                                                   >  1:2006304/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:2187730/100‑1 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATc                                                                  <  2:1538557/100‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCCGCGCCAACAATCGTTTAt                                                                   >  1:2487732/1‑98 (MQ=255)

                                tGGTTCTCGACGCACATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTGGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATcc                                                                 >  1:1247628/1‑100 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCAttt                                                             <  2:1573931/100‑1 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  1:426936/1‑100 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCAC‑C‑CTCCCCGCGCCAAACATAGTTTAt                                                                   >  1:1510343/1‑94 (MQ=255)

                                    nctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  2:2623567/2‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            >  1:1411465/1‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:2259769/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1391591/100‑1 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:2851319/1‑100 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:84073/100‑1 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  1:1938480/1‑100 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:2209968/1‑100 (MQ=255)

                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  1:1849243/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:1489722/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:1254705/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:887690/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     <  2:456435/100‑1 (MQ=255)

                                            cccATTGTGGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGTATGTTCTCGTTT‑‑CCC‑C‑CGCCGCCAGCGCCAACCATAGTTTATCCATTTTCGCCCgg                                                     >  1:448283/1‑100 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  2:1373321/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  2:2105939/1‑100 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:1492665/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:818220/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  2:1835355/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:816493/100‑1 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTCCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCCGCGCCAACCATAGTTTATCCCTTTGCGCCCGGGTTTc                                                >  1:2865547/1‑95 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  2:1179527/1‑100 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  1:294672/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  2:1988066/100‑1 (MQ=255)

                                                               gATGGTTGTTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:2689787/100‑1 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:985412/100‑1 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  2:2987291/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:2334221/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:871278/100‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:1005732/1‑100 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           <  1:1180851/100‑1 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:1374523/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          >  1:2192467/1‑100 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:388209/100‑1 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  2:2417167/100‑1 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  2:646931/100‑1 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  2:638114/100‑1 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCAATTGCGCCCGGTTTTCATCGTGAAGTGGCTAAAGCCACAGacg                     <  2:2063190/100‑3 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:2976017/100‑1 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:707181/100‑1 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   >  1:1303225/1‑100 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:1651892/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:2364176/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:2902890/100‑1 (MQ=255)

                                                                                   ggTGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGACAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTCTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:1841955/97‑1 (MQ=255)

                                                                                     tGAATTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  2:1620307/1‑100 (MQ=255)

                                                                                       aaTTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgcgtg           >  1:608152/1‑99 (MQ=255)

                                                                                        gttttCTCGTTT‑‑CCC‑CGC‑CCGCCCGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtaagtgag         >  1:2137415/2‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:1473395/100‑1 (MQ=255)

                                                                                         ttttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  2:1152400/100‑1 (MQ=255)

                                                                                         ttttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  1:675070/100‑1 (MQ=255)

                                                                                         ttttCTCGGTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  2:3040210/100‑1 (MQ=255)

                                                                                          tttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       >  2:1200935/1‑100 (MQ=255)

                                                                                          tttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:895404/100‑1 (MQ=255)

                                                                                           ttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:1020748/99‑1 (MQ=255)

                                                                                           ttCTCGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      <  1:1505790/100‑1 (MQ=255)

                                                                                             ctcGTTTC‑CCC‑C‑C‑CCGCCAGCGCCATCCATAGTTTGTCCATTTGCGCCCGTGTTTCATCGTGCAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  1:2920777/100‑1 (MQ=255)

                                                                                             ctcGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  1:310504/100‑1 (MQ=255)

                                                                                              tcGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:426940/100‑1 (MQ=255)

                                                                                               cGTTTC‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACgn  <  1:2623552/100‑2 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCC‑C‑C‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑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: `‑`