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	93.8%	414.7 / 20.8	129	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (6/2); new base C (69/50); total (77/52)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑ccc                                                                                             <  2:53030/100‑1 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc                                                                                            >  1:1076734/1‑100 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc                                                                                            <  1:2008626/100‑1 (MQ=255)

      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑c                                                                                          >  2:2768845/1‑100 (MQ=255)

      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑c                                                                                          >  2:2393612/1‑100 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑ccc                                                                                        >  2:3678867/1‑99 (MQ=255)

        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑ccc                                                                                        >  1:2579852/1‑99 (MQ=255)

         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑cccg                                                                                       <  2:1572043/100‑3 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑cccgc                                                                                      >  2:2867093/1‑97 (MQ=255)

           gtaGAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑cccgcc                                                                                     >  1:1707834/1‑96 (MQ=255)

           gtaAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑cccgcc                                                                                     <  2:3326912/100‑5 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c‑cccc‑cccg                                                                                       >  1:3104256/1‑96 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAgc                                                                                   <  1:1194846/100‑1 (MQ=255)

               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAgcg                                                                                  <  2:1267662/100‑1 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAgcgc                                                                                 <  2:1040287/100‑1 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGGTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGcc                                                                                >  1:774889/1‑100 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGcc                                                                                >  1:2708699/1‑100 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCaa                                                                              >  2:712862/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAAcc                                                                            >  1:2370237/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCa                                                                           >  2:441153/1‑100 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCAt                                                                          >  2:2794603/1‑100 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCAt                                                                          >  2:263279/1‑100 (MQ=255)

                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATa                                                                         >  2:1626132/1‑100 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAg                                                                        <  1:2726200/100‑1 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAg                                                                        <  2:2837810/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCTTT‑TC‑C‑CCCC‑CCGCCAGCGCCaan                                                                             <  1:244981/93‑2 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGtt                                                                      <  2:1085337/100‑1 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGttt                                                                     >  2:2019405/1‑100 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGttt                                                                     >  2:135644/1‑100 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCACCGCCAACCATa                                                                         >  2:3750897/1‑96 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:1514916/100‑1 (MQ=255)

                                 ggTTCTCGACGCGCATTGTTGCTAATGGCCGATGGTGGATTTCTAACGTAGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTAt                                                                   <  1:1291570/97‑1 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCa                                                                >  1:1795296/1‑100 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCa                                                                >  1:3531953/1‑100 (MQ=255)

                                 cgTTCTCGACATCCATTGTTGGAGATGGCCGATGGTTGCTTTCTAACGTCGGTGGATTTTCTCGTT‑T‑‑C‑CCCCTCCGTCAGCGCCAACCATAGTTTAttta                                                                <  1:2489016/99‑4 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  2:2587667/1‑100 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TG‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            >  2:2666270/1‑100 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCCGCGCCAACCATAGTTTGTCCCTTTgcg                                                          >  2:3636377/1‑100 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:1301747/100‑1 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:682267/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:2263471/1‑100 (MQ=255)

                                              naTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGAT‑T‑‑C‑CCCCTCCGCCAACGCCAACCACAGTTTATCCATTTGCGCCCg                                                      >  1:817658/2‑97 (MQ=255)

                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   <  2:1542219/100‑1 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAACGCCAACCATAGTTTATCCATTTCCGCCCGGGtt                                                  >  1:1384689/1‑100 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAACTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:2667368/100‑1 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑CACCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  2:3008816/1‑99 (MQ=255)

                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:2955126/1‑100 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  2:2139901/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  1:3340026/100‑1 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:1417661/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGCGAATTTGCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  1:2020635/1‑100 (MQ=255)

                                                            gctgATGCTTGATTTCTAACGACGGTGAACTTTTTCGTT‑TA‑C‑ACCC‑CCACCAGCGCCAACCAAAGTTTTTCCATTTGCGCCCGTGTTTCATGGTGAAGTg                                     >  1:900970/4‑100 (MQ=255)

                                                            gCCGGTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTGC‑‑C‑CCCC‑CCACCGGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:3605328/96‑1 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  2:2037903/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  1:2060133/1‑100 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:100138/100‑1 (MQ=255)

                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:714279/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCAAAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  >  1:3791270/1‑99 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  1:3782/1‑100 (MQ=255)

                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  2:3252594/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:2975465/1‑100 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:3792928/1‑100 (MQ=255)

                                                                    ttcaTTTCTAACGTCGGTTAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:3352175/95‑1 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:1849170/1‑100 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  1:859997/100‑1 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:3680003/1‑100 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:1642758/1‑100 (MQ=255)

                                                                         ttCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  2:2372211/1‑100 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGCGGCGAAAGCCACAGa                       >  2:3747264/1‑100 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  2:1729679/1‑100 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  2:185552/100‑1 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  1:1496151/1‑100 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGAGCCAACTATTGTTTATCCATTTGCGCCCGGGGTTCATCGTGAAGTGGCGCAAGCCACAGATgcg                   >  1:2604784/1‑100 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  >  2:3369108/1‑100 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTT‑TCCC‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  <  1:2768828/100‑1 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:2407420/100‑1 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  2:1310879/1‑100 (MQ=255)

                                                                                   ggTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  1:220247/1‑100 (MQ=255)

                                                                                    gTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:2064962/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:322377/1‑100 (MQ=255)

                                                                                    gTGAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:3273927/100‑1 (MQ=255)

                                                                                      gcaTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACATATGCGCCtgagtg           <  2:2069737/98‑1 (MQ=255)

                                                                                      gAATTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:3200932/100‑1 (MQ=255)

                                                                                       aaTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACGATAGTTTATCTATTTGTGCCCGGGTTTCAACGGGAAGTGGCGAAAGGCACAGCTGCGCCtg               >  1:3726924/1‑95 (MQ=255)

                                                                                       aaTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          >  1:1375543/1‑100 (MQ=255)

                                                                                       aaTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          >  1:74033/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:3110770/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:1206661/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2126769/100‑1 (MQ=255)

                                                                                        aTTTTCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2126959/1‑100 (MQ=255)

                                                                                         ttttCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  2:3531980/100‑1 (MQ=255)

                                                                                         ttttCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  1:3548949/100‑1 (MQ=255)

                                                                                          tttCTCGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:2642636/100‑1 (MQ=255)

                                                                                            tctcGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:2274048/1‑100 (MQ=255)

                                                                                             ctcGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  2:2859105/100‑1 (MQ=255)

                                                                                              tcGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:3138251/100‑1 (MQ=255)

                                                                                               cGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGATTGTACGc  >  1:2537188/1‑100 (MQ=255)

                                                                                               cGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGCGTAAGc  >  2:618142/1‑100 (MQ=255)

                                                                                               cGTT‑TC‑C‑CCCC‑CCGCCAGCGCCAACCATAGTTGATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:3608176/1‑100 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑C‑CCCC‑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: `‑`