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.8%	401.1 / 15.4	116	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (4/2); new base C (44/66); total (48/68)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 139 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCA  >  NC_000913/3379160‑3379354

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

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

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

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

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑cc                                                                                              >  1:2438611/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑cc                                                                                              >  1:1169640/1‑100 (MQ=255)

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

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑cc                                                                                              >  1:443532/1‑100 (MQ=255)

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

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑cc                                                                                              <  2:5988/99‑1 (MQ=255)

      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑ccc‑c                                                                                           <  1:1549398/100‑1 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑ccc‑cccg                                                                                        >  1:1237184/1‑98 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑ccc‑ccccgc                                                                                      >  2:1585551/1‑96 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAg                                                                                     >  2:1067855/1‑100 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAg                                                                                     >  1:2607889/1‑100 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑c‑‑c‑ccc‑ccccgccagc                                                                                  <  2:2478405/100‑9 (MQ=255)

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

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTAGCAAATGGCCGATGGTTGATTTCTAACGTAGATGAATTTTCTCGTTT‑‑‑a‑‑c‑ccc‑ccccgc                                                                                      >  1:1447577/1‑91 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCaa                                                                               >  1:2895768/1‑100 (MQ=255)

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

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAAc                                                                              >  2:1733560/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAAc                                                                              >  2:2536680/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAAcc                                                                             <  1:2505298/100‑1 (MQ=255)

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

                      cGCAAGCGTCTGGTTCTAGACGCTCATTGTTGCTAATGGCCGATGGTTAATTTCTGACGTCGGTTAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAAcc                                                                             <  1:577037/100‑1 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATa                                                                          >  1:2614825/1‑100 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑C‑CCCACCGCCAGCGCCAACCATa                                                                          >  1:3193979/1‑100 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCATTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCAAAACATAg                                                                         >  2:1341855/1‑98 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGttt                                                                      <  2:1388496/100‑1 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGCTGGTTGATTTCTAACGCCGCTGCATTTTCTCGTTT‑‑‑C‑‑CGCGC‑CCGCCAGCGCAAAACATAAATTATc                                                                   >  1:1138932/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGATGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATc                                                                   >  1:1847411/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATcc                                                                  <  2:2911906/100‑1 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATcc                                                                  >  1:1731379/1‑100 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCAtt                                                               >  2:3004006/1‑100 (MQ=255)

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

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

                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             >  2:2782170/1‑100 (MQ=255)

                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             <  2:1742812/100‑1 (MQ=255)

                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCAGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             <  2:157541/100‑1 (MQ=255)

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

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

                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                       <  1:530723/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                      <  2:264807/100‑1 (MQ=255)

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

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                   <  2:2918240/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                 <  2:2808763/100‑1 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                                <  2:498203/100‑1 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                               >  1:2771509/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAATTTATCCATTTGCGCCCGGGTTTCacc                                              >  1:257060/1‑98 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                          >  2:2377073/1‑100 (MQ=255)

                                                        aaTGGCCGTTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                           <  1:38618/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCC‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTTCCCCCGGGTTTCCTCGt                                            >  1:1733905/1‑100 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                             <  1:2536661/98‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                        >  1:1300366/1‑100 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCATCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                    <  1:2292456/99‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                   <  1:1968074/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                  <  2:2950469/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCCTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                                <  1:331658/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                               <  2:1208488/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                               <  1:240371/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                              <  2:883814/100‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                             <  2:927670/100‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                             <  1:653097/100‑1 (MQ=255)

                                                                       aTTTCTAACGTCGGTTAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                             <  1:911216/99‑1 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                            >  2:2208706/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                           <  2:1859256/100‑1 (MQ=255)

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

                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                         <  2:149286/100‑1 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                         <  1:1733544/100‑1 (MQ=255)

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

                                                                           cTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                        <  1:1113026/100‑1 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                        >  2:683218/1‑100 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                       >  2:1447382/1‑100 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                        >  2:1401530/1‑98 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                      <  2:2440376/100‑1 (MQ=255)

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

                                                                                   ggTGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                  <  2:113554/98‑1 (MQ=255)

                                                                                   ggTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCACCAGCGCCAACCATTGTTTATCCATTTGCCCCCGGGTTTCTTCGTGAAGTTGCAAAAGCCACAGATGCGCCtga               >  2:288461/1‑100 (MQ=255)

                                                                                    gTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  2:1300378/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  2:2771531/100‑1 (MQ=255)

                                                                                     tGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  1:1379575/100‑1 (MQ=255)

                                                                                     tGAATTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  1:931627/100‑1 (MQ=255)

                                                                                      gAATTTTCTCTTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGTCAACCATAGTTTCTCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            <  1:1341843/100‑1 (MQ=255)

                                                                                       aaTTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            <  2:1731395/100‑1 (MQ=255)

                                                                                       aaTTTTCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            >  2:74466/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag          >  2:2956599/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag          >  2:3210978/1‑100 (MQ=255)

                                                                                         ttttCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag          <  1:559825/100‑1 (MQ=255)

                                                                                         ttttCTCGTTT‑‑‑CCGC‑CCC‑CCGCCAGCGCCAACCATAGTGTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag          >  1:2132790/1‑100 (MQ=255)

                                                                                          ttttTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            <  1:1455660/96‑1 (MQ=255)

                                                                                          tttCTCGTTTCC‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt         >  1:1535158/1‑100 (MQ=255)

                                                                                          tttCTCGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtggctgag          >  2:2385918/1‑98 (MQ=255)

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

                                                                                            tctcGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa      >  2:3248287/1‑100 (MQ=255)

                                                                                            tctcGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa      <  2:1001620/100‑1 (MQ=255)

                                                                                            tctcGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCCTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa      >  1:1432429/1‑100 (MQ=255)

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

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

                                                                                               cGTTTC‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc   <  2:455140/100‑1 (MQ=255)

                                                                                                ttct‑‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATGCATTTGCGCCCGGGTTTCAACGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa  <  1:2782063/96‑1 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑‑C‑‑C‑CCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCA  >  NC_000913/3379160‑3379354

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: `‑`