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	95.4%	451.3 / 14.4	131	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (64/60); total (67/64)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379162‑3379353

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

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

  caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc                                                                                              <  2:631548/100‑1 (MQ=255)

   aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑c                                                                                            <  2:295473/100‑1 (MQ=255)

    gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑cc                                                                                           <  2:2703526/100‑1 (MQ=255)

      gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑cccc                                                                                         <  1:2724207/100‑2 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑ccccgcc                                                                                      <  2:3132272/100‑5 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑ccccgcc                                                                                      <  1:778601/100‑5 (MQ=255)

         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑ccccgcc                                                                                      <  1:1097539/100‑5 (MQ=255)

          taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccc‑ccccgcca                                                                                     >  2:2147250/1‑95 (MQ=255)

           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAg                                                                                    <  1:1334475/100‑1 (MQ=255)

            aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAgc                                                                                   >  1:1124943/1‑100 (MQ=255)

              agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAgcgc                                                                                 >  2:928105/1‑100 (MQ=255)

              agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAgcgc                                                                                 >  2:895119/1‑100 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAAcc                                                                            >  1:953501/1‑100 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAAcc                                                                            <  1:2307241/100‑1 (MQ=255)

                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑ACCAGCGCCAAcc                                                                            >  1:65804/1‑100 (MQ=255)

                    cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑CCC‑GCCAGCGCCAAcc                                                                            >  2:2014137/1‑100 (MQ=255)

                     gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCa                                                                           <  1:502714/99‑1 (MQ=255)

                        agtgTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCC‑CCC‑GCCAGCGCCAACCATAg                                                                        <  1:254828/97‑1 (MQ=255)

                         gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGtt                                                                      >  1:1754075/1‑100 (MQ=255)

                             cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATc                                                                  <  2:2069960/100‑1 (MQ=255)

                             cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCTCCC‑GCCAGCGCCAACCATAGTTTATc                                                                  <  1:1971064/100‑1 (MQ=255)

                              tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATcc                                                                 <  2:1389662/100‑1 (MQ=255)

                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCa                                                                >  2:439298/1‑100 (MQ=255)

                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCa                                                                >  2:2731857/1‑100 (MQ=255)

                                gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTGTCTCGTTT‑‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCAtt                                                              >  2:2281576/1‑100 (MQ=255)

                                 ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCAtt                                                              <  2:1000241/100‑1 (MQ=255)

                                  tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCAttt                                                             <  2:2552217/100‑1 (MQ=255)

                                  tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCAttt                                                             >  1:1450931/1‑100 (MQ=255)

                                   ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:3179096/100‑1 (MQ=255)

                                   ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1805572/100‑1 (MQ=255)

                                    tccgcGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2883621/95‑1 (MQ=255)

                                     cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  2:1237676/1‑100 (MQ=255)

                                        cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:398836/1‑100 (MQ=255)

                                        cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:1691441/1‑100 (MQ=255)

                                          cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     <  1:175748/100‑1 (MQ=255)

                                           ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  1:2099203/1‑100 (MQ=255)

                                           ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  1:2503140/1‑100 (MQ=255)

                                           ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:2812181/100‑1 (MQ=255)

                                           ccATTGTTGCCAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  2:1877456/100‑1 (MQ=255)

                                            taTTGTTGCTAATGGTCGATGGTTTATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:264899/98‑1 (MQ=255)

                                              ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  2:490422/100‑1 (MQ=255)

                                               tgttgCTAATGGCCGATTGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  1:899993/100‑1 (MQ=255)

                                               tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  1:2632861/1‑100 (MQ=255)

                                                gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               >  1:1625205/1‑100 (MQ=255)

                                                gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               <  1:1098727/100‑1 (MQ=255)

                                                 ttgCTAATGGCCGATGGTTGATTTATAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCACGGGttt                                                 >  2:233945/1‑97 (MQ=255)

                                                  tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:3152163/1‑100 (MQ=255)

                                                  tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:2206542/1‑100 (MQ=255)

                                                  tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:3115885/1‑100 (MQ=255)

                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  1:1661843/1‑100 (MQ=255)

                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  2:1523460/1‑100 (MQ=255)

                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  2:2099188/100‑1 (MQ=255)

                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  2:608063/100‑1 (MQ=255)

                                                     tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  1:736879/100‑1 (MQ=255)

                                                      aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:492559/100‑1 (MQ=255)

                                                        tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:65802/100‑1 (MQ=255)

                                                         ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  2:1291376/1‑100 (MQ=255)

                                                           ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  2:2476560/100‑1 (MQ=255)

                                                            cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:1754063/100‑1 (MQ=255)

                                                              aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                 >  1:2810461/1‑100 (MQ=255)

                                                               tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                <  2:2568776/100‑1 (MQ=255)

                                                                ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  2:1092651/1‑100 (MQ=255)

                                                                 gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:38105/1‑100 (MQ=255)

                                                                 gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:341194/100‑1 (MQ=255)

                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTAGCGAAAGc                             >  2:3233569/1‑100 (MQ=255)

                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:1295288/100‑1 (MQ=255)

                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  1:1502729/1‑100 (MQ=255)

                                                                    gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:2676067/1‑100 (MQ=255)

                                                                    gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:3061404/1‑100 (MQ=255)

                                                                      tttCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:1106185/100‑1 (MQ=255)

                                                                        tCTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  1:414273/1‑100 (MQ=255)

                                                                         cTAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  1:785451/1‑100 (MQ=255)

                                                                          tAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  1:2148278/1‑100 (MQ=255)

                                                                          tAACGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:2730313/1‑100 (MQ=255)

                                                                            aCGTCGGTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:903715/100‑1 (MQ=255)

                                                                                 ggTGAATTTTTTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:901456/99‑1 (MQ=255)

                                                                                 ggTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              <  1:439300/100‑1 (MQ=255)

                                                                                 ggTAAATTTTCTCGTTT‑‑CCCC‑CCCTGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:3227313/99‑1 (MQ=255)

                                                                                  gTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:754199/1‑100 (MQ=255)

                                                                                  gTGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:2264004/100‑1 (MQ=255)

                                                                                   tGAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  1:2854101/1‑100 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:1135672/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:1506424/100‑1 (MQ=255)

                                                                                    gAATTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATAGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:1596692/100‑1 (MQ=255)

                                                                                      aTTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:1156642/100‑1 (MQ=255)

                                                                                      aTTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:2503122/100‑1 (MQ=255)

                                                                                      aTTTTCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2286614/1‑100 (MQ=255)

                                                                                       ttttCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:1648726/1‑100 (MQ=255)

                                                                                        tttCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:1098072/100‑1 (MQ=255)

                                                                                         ttCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCGTTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:2041286/1‑100 (MQ=255)

                                                                                         ttCTCGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:2314841/1‑100 (MQ=255)

                                                                                          tctcGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:1205302/1‑100 (MQ=255)

                                                                                           ctcGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:1355849/1‑100 (MQ=255)

                                                                                           ctcGTTTC‑CCCC‑CCC‑GCCAGCGCCAAACATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:351520/1‑100 (MQ=255)

                                                                                           ctcGTTT‑‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:2384165/100‑1 (MQ=255)

                                                                                            tcGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:2384125/100‑1 (MQ=255)

                                                                                            tcGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:1912749/100‑1 (MQ=255)

                                                                                             cGTTTC‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGAGAGAGTActc  >  1:58701/1‑98 (MQ=255)

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCC‑CCC‑GCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379162‑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: `‑`