BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)_7→8	coding (732/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,259	1	.	C	93.8%	374.3 / 16.5	113	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (4/3); new base C (42/64); total (46/67)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 132 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGtt‑t                                                                                                    >  1:2792733/1‑100 (MQ=255)

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGtt‑t                                                                                                    >  1:1505701/1‑100 (MQ=255)

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGtt‑t                                                                                                    >  1:55338/1‑100 (MQ=255)

 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c                                                                                                 <  2:3030134/100‑1 (MQ=255)

 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c                                                                                                 >  1:1770709/1‑100 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c                                                                                               >  1:23588/1‑100 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c                                                                                               >  1:556837/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑cc                                                                                            >  1:2386019/1‑100 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccc                                                                                           <  2:2298661/100‑1 (MQ=255)

       cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccccc                                                                                         <  2:2384028/100‑1 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccccccg                                                                                       <  2:375750/100‑3 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccccccgcc                                                                                     >  1:2274064/1‑96 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccccccgcca                                                                                    >  2:2752051/1‑95 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAg                                                                                    >  1:32414/1‑100 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAg                                                                                    >  1:2841357/1‑100 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAgcg                                                                                  >  1:2889921/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAgcgc                                                                                 >  1:729750/1‑100 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑c‑ccccccgaccgcgcc                                                                               >  2:2832075/1‑90 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCa                                                                               <  1:3076859/100‑1 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCaa                                                                              <  2:644554/100‑1 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAAcc                                                                            >  1:594468/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAAcc                                                                            <  1:2103148/100‑1 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATTGCCGATGGTGGATTTCTAACGTCGGTGAATGTTCTCGTG‑T‑‑CTC‑CCCCCGCCAGTTCCAACTAt                                                                          >  1:139382/1‑100 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATa                                                                         >  1:3000608/1‑100 (MQ=255)

                          aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAg                                                                        <  2:2889892/99‑1 (MQ=255)

                          aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGt                                                                       <  1:2796484/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGtt                                                                      <  1:1845769/100‑1 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGttt                                                                     <  2:1010156/100‑1 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTAt                                                                   <  2:597514/99‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 >  2:580365/1‑99 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 >  1:2720286/1‑99 (MQ=255)

                                  gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                               >  1:2653351/1‑100 (MQ=255)

                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              >  1:1649585/1‑100 (MQ=255)

                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              <  1:1859347/100‑1 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             <  1:1256320/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  2:1331182/1‑99 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            >  2:1743890/1‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:1321421/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1852068/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2686564/100‑1 (MQ=255)

                                         aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  2:2836776/100‑1 (MQ=255)

                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:2511412/1‑100 (MQ=255)

                                           gCCCATTGTGGCTAATGGCCGAGGGTTGATTTCTAAGGTGGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  2:1008213/96‑1 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:225161/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:2996942/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     <  1:299911/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:1353185/100‑1 (MQ=255)

                                              cATTGTTGCTAATGGACGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   <  1:3037638/100‑1 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               <  1:1331193/100‑1 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                              <  1:1121890/100‑1 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                              >  1:2302049/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            >  2:137358/1‑100 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  1:386084/100‑1 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  1:2382529/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:1659723/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  2:3000576/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑CGCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCTCCCGGGTTTCATCGTGa                                         >  2:1721540/1‑100 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  1:2635919/1‑100 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:471296/100‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:2554268/100‑1 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCGCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGGGCCCGGGTTTCATCGTGAAGTg                                     <  2:2665533/100‑1 (MQ=255)

                                                             ccGATGATTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  2:2049168/100‑1 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  2:2985558/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:2124134/100‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  >  1:2000265/1‑100 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:1045635/100‑1 (MQ=255)

                                                                 tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                >  2:1806289/1‑100 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:2845068/1‑100 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  1:434864/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑TCCC‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:1770690/100‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCATCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:1979729/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:2383719/100‑1 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  2:1873898/1‑100 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:382673/1‑100 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  1:1539656/100‑1 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:467622/100‑1 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  <  1:1659741/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:2432235/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:497574/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 >  1:1511993/1‑100 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:2112376/100‑1 (MQ=255)

                                                                                   ggTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  2:1956326/1‑100 (MQ=255)

                                                                                    gTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:1202317/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCAAAAGCCACAGATGCGCCtgag             <  2:1056422/100‑1 (MQ=255)

                                                                                      gAATTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:345342/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2450594/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:124881/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:238624/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2450367/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2788682/100‑1 (MQ=255)

                                                                                         tgttCTCGTTCC‑‑C‑C‑CCCCCGCCAGCGCCACCCATCGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:255150/97‑1 (MQ=255)

                                                                                          tttCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:1262785/100‑1 (MQ=255)

                                                                                          tttCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:779840/100‑1 (MQ=255)

                                                                                          tttCTCGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:1578812/100‑1 (MQ=255)

                                                                                            tctcGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:1511978/100‑1 (MQ=255)

                                                                                            tctcGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:55338/100‑1 (MQ=255)

                                                                                            tctcGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATATGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:139382/100‑1 (MQ=255)

                                                                                             ctcGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGCGAGTAc    >  1:2111440/1‑100 (MQ=255)

                                                                                             ctcGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  1:203629/1‑100 (MQ=255)

                                                                                               cGTT‑TC‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  1:2413113/100‑1 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑C‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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