BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)_7→6	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,260	0	C	.	98.1%	490.4 / 1.4	104	coding (738/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base C (2/0); new base . (60/42); total (62/42)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTc                                                                                               <  2:1067143/100‑1 (MQ=255)

gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTc                                                                                               <  2:1974359/100‑1 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTcc                                                                                              >  1:901888/1‑100 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTcc                                                                                              >  2:1064216/1‑100 (MQ=255)

     gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTcccccc                                                                                          >  1:570743/1‑100 (MQ=255)

      cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTccccccg                                                                                         <  1:2076410/100‑2 (MQ=255)

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTccccccgcc                                                                                       <  1:688706/100‑4 (MQ=255)

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTccccccgcc                                                                                       <  1:2085356/100‑4 (MQ=255)

         agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTccccccgcca                                                                                      <  1:686179/100‑5 (MQ=255)

          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAg                                                                                    >  2:1107247/1‑100 (MQ=255)

          gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGACGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCGCCGCCAg                                                                                    >  2:2257372/1‑100 (MQ=255)

           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgc                                                                                   <  1:2385545/100‑1 (MQ=255)

           taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgc                                                                                   >  2:1874266/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgcg                                                                                  >  2:577976/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgcg                                                                                  >  1:1817941/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgcg                                                                                  >  1:2355379/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAgcg                                                                                  >  2:434228/1‑100 (MQ=255)

              aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGcc                                                                                <  1:2029586/100‑1 (MQ=255)

              aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGcc                                                                                >  1:1656689/1‑100 (MQ=255)

               tccgatctccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCa                                                                               <  2:1795462/92‑1 (MQ=255)

               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCa                                                                               >  2:1644794/1‑100 (MQ=255)

                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAAcc                                                                            <  2:901889/100‑1 (MQ=255)

                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAAcc                                                                            <  1:1963954/100‑1 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCAt                                                                          >  2:692148/1‑100 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAg                                                                        <  1:577975/100‑1 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAg                                                                        >  1:1813723/1‑100 (MQ=255)

                       ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCaagatcgga                                                                       >  1:1795461/1‑93 (MQ=255)

                       ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGt                                                                       <  1:36022/100‑1 (MQ=255)

                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGtt                                                                      <  1:2099818/100‑1 (MQ=255)

                         aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGttt                                                                     >  2:1034690/1‑100 (MQ=255)

                           cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTAt                                                                   >  2:1110159/1‑100 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:2172816/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATc                                                                  <  2:1354813/100‑1 (MQ=255)

                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                <  1:75667/100‑1 (MQ=255)

                               tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                               >  2:1733525/1‑100 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              >  1:1940029/1‑100 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              >  1:971115/1‑100 (MQ=255)

                                 gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  1:2268484/1‑100 (MQ=255)

                                 gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             >  2:436525/1‑100 (MQ=255)

                                 gTTCGCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              <  2:2172225/99‑1 (MQ=255)

                                  ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1645106/100‑1 (MQ=255)

                                  ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1644793/100‑1 (MQ=255)

                                  ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1488164/100‑1 (MQ=255)

                                   tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  2:2165779/1‑100 (MQ=255)

                                     ncGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTCTCCATTTgcgc                                                         >  1:1017741/2‑100 (MQ=255)

                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  2:1940028/100‑1 (MQ=255)

                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        >  2:347462/1‑100 (MQ=255)

                                        aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:2108974/1‑100 (MQ=255)

                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:12454/1‑100 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  2:1506005/1‑100 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  1:829602/1‑100 (MQ=255)

                                           nccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   >  1:1912623/2‑100 (MQ=255)

                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   >  2:2184211/1‑100 (MQ=255)

                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   >  1:2085173/1‑100 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:919650/100‑1 (MQ=255)

                                              aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  2:706218/1‑100 (MQ=255)

                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgggagatcg                                              >  1:2360738/1‑94 (MQ=255)

                                                 gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:545598/1‑100 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  2:95813/1‑100 (MQ=255)

                                                    gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:2174671/1‑100 (MQ=255)

                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCGGGGTTTCATCGTGaa                                        >  1:2363435/1‑100 (MQ=255)

                                                        aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  2:922886/1‑100 (MQ=255)

                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:1352868/1‑100 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  2:636376/1‑100 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                >  1:50084/1‑100 (MQ=255)

                                                               cgatctGATTTCTAATGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:2363264/95‑1 (MQ=255)

                                                               aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:1506004/100‑1 (MQ=255)

                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:1279072/1‑100 (MQ=255)

                                                                 ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  2:2085171/100‑1 (MQ=255)

                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:1506469/1‑100 (MQ=255)

                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:815989/1‑100 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:40050/1‑100 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           <  2:829603/100‑1 (MQ=255)

                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  2:673001/100‑1 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaaagatcg                         >  2:2363261/1‑94 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:530165/100‑1 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:302556/100‑1 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:2184213/100‑1 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  1:786081/1‑100 (MQ=255)

                                                                       tttCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  2:944268/1‑100 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  2:511324/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  1:1333723/100‑1 (MQ=255)

                                                                         tCTAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCGACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:642569/1‑100 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   >  2:261759/1‑100 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:651640/100‑1 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgcn                 <  2:1912622/100‑2 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:786082/100‑1 (MQ=255)

                                                                                 cGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:1081811/100‑1 (MQ=255)

                                                                                 cGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  1:181560/1‑100 (MQ=255)

                                                                                 cGGTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  2:258528/100‑1 (MQ=255)

                                                                                  ggTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  2:1297920/1‑100 (MQ=255)

                                                                                   gTGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:2174673/100‑1 (MQ=255)

                                                                                    tGAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          >  2:1905386/1‑100 (MQ=255)

                                                                                     gAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:1944915/100‑1 (MQ=255)

                                                                                     gAATTTTCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:1241514/1‑100 (MQ=255)

                                                                                        ttttCTCGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:689334/1‑100 (MQ=255)

                                                                                           tctcGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:2257374/100‑1 (MQ=255)

                                                                                            ctcGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  1:338856/100‑1 (MQ=255)

                                                                                            ctcGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:1307712/100‑1 (MQ=255)

                                                                                            ctcGTTT‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:1567471/1‑100 (MQ=255)

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

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