BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	284,871	(G)_6→7	coding (253/1383 nt)	yagG →	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	284,865	1	.	G	96.3%	649.0 / 15.4	214	coding (247/1383 nt)	yagG	putative permease
Reads supporting (aligned to +/- strand): ref base . (1/7); new base G (86/120); total (87/127)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 269 total aligned reads displayed.

CCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGC  >  NC_000913/284774‑284951

ccGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCt‑g‑t‑‑‑g                                                                                         <  1:2901390/93‑1 (MQ=255)

ccGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCt‑g‑t‑‑‑g                                                                                         >  1:6944711/1‑93 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  2:5615607/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  2:1381627/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  1:6058966/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  2:5358968/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  2:6355993/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      <  2:792951/93‑1 (MQ=255)

  gcgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg                                                                                      >  2:6138474/1‑93 (MQ=255)

   cgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑g                                                                                    <  2:759678/93‑1 (MQ=255)

   cgTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑g                                                                                    <  1:370724/93‑1 (MQ=255)

    gTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑gg                                                                                   <  2:7331889/93‑1 (MQ=255)

    gTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑gg                                                                                   >  2:2953030/1‑93 (MQ=255)

    gTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑gg                                                                                   <  1:3669552/93‑1 (MQ=255)

    gTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑gg                                                                                   <  2:1077922/93‑1 (MQ=255)

     tGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑ggg                                                                                  <  2:2501153/93‑1 (MQ=255)

     tGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑ggg                                                                                  <  1:5622454/93‑1 (MQ=255)

       cTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑ggggc                                                                                <  1:7425406/93‑3 (MQ=255)

       cTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑ggggc                                                                                <  2:5209519/93‑3 (MQ=255)

        tCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑g‑gg‑ggggcc                                                                               >  1:3366275/1‑90 (MQ=255)

             gCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATccc                                                                           <  2:4385987/93‑1 (MQ=255)

               cGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGt                                                                         >  1:1492757/1‑93 (MQ=255)

                gTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGtt                                                                        <  2:4861290/93‑1 (MQ=255)

                ctCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGtt                                                                        <  2:5828242/92‑1 (MQ=255)

                 tCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTc                                                                       >  1:7633909/1‑93 (MQ=255)

                 tCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTc                                                                       >  2:5961733/1‑93 (MQ=255)

                  cACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGG‑G‑GG‑GGGCCATCCCGTTc                                                                       >  2:5883896/1‑93 (MQ=255)

                   aCCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCgg                                                                     >  1:6798478/1‑93 (MQ=255)

                   aCCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCgg                                                                     <  2:4222878/93‑1 (MQ=255)

                   aCCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCgg                                                                     >  1:1112141/1‑93 (MQ=255)

                     cGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCa                                                                   <  2:2579756/93‑1 (MQ=255)

                     cGACCCGCTGATGGCGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCa                                                                   <  1:220977/93‑1 (MQ=255)

                      gACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGG‑G‑GG‑GGGCCATCCCGTTCGGCa                                                                   <  2:2991223/93‑1 (MQ=255)

                      gACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCAt                                                                  <  2:6629340/93‑1 (MQ=255)

                       aCCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATc                                                                 <  1:3483258/93‑1 (MQ=255)

                         ccGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGt                                                               >  1:123149/1‑93 (MQ=255)

                         ccGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTc                                                              <  1:3314082/93‑1 (MQ=255)

                          cGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTc                                                              <  2:1048434/93‑1 (MQ=255)

                           gccgATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCt                                                             <  1:6677137/90‑1 (MQ=255)

                             cgATGGGTCTGCTGGTAGTCCGCACCCGCACGCGGCACGCCCAGTTCCGCCCGTTCCTGCT‑GTT‑‑‑G‑GG‑GGGCCTTCCCGTTCGGCATCTTCTGc                                                           <  2:4131093/92‑1 (MQ=255)

                               aTGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGt                                                         <  1:1012621/93‑1 (MQ=255)

                                  gggCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCt                                                      >  1:2944405/1‑93 (MQ=255)

                                    gctgctGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGa                                                    <  1:2586802/93‑1 (MQ=255)

                                     ctgctgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGAc                                                   >  1:1613902/1‑93 (MQ=255)

                                       gctgGTAGACCGCTCCCGCGCGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGGGG‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGAc                                                   <  2:3566198/93‑1 (MQ=255)

                                         tgGTAGACCGCCCCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTc                                               <  2:1190599/93‑1 (MQ=255)

                                         tgGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTc                                               <  2:574163/93‑1 (MQ=255)

                                          ggTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCt                                              <  1:2997339/93‑1 (MQ=255)

                                            tAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTac                                            <  1:1753735/93‑1 (MQ=255)

                                             aGACCGCCCCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTaca                                           <  1:5448741/93‑1 (MQ=255)

                                             aGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTaca                                           <  1:2661482/93‑1 (MQ=255)

                                             aGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCCGACCTTCTaca                                           <  2:4564281/93‑1 (MQ=255)

                                              gACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGG‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTaca                                           >  1:7724810/1‑93 (MQ=255)

                                              gACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTacac                                          <  2:3052651/93‑1 (MQ=255)

                                                ccgcgcccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGc                                        <  2:3768209/93‑1 (MQ=255)

                                                 cgcacccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGG‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGc                                        <  2:3453316/93‑1 (MQ=255)

                                                  ccccccgcgcGCGGCGCGGCCAGTTCCGCCCGTTCCTGCTGG‑G‑‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCg                                      <  1:1977929/90‑1 (MQ=255)

                                                    acccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGa                                    >  1:4416492/1‑93 (MQ=255)

                                                     cccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGAc                                   >  2:6458271/1‑93 (MQ=255)

                                                     cccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGAc                                   >  1:4519404/1‑93 (MQ=255)

                                                     cccgcacGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGAc                                   <  1:3086760/93‑1 (MQ=255)

                                                          acGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtctc                              <  2:2412805/93‑1 (MQ=255)

                                                          acGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtctc                              <  1:2492793/93‑1 (MQ=255)

                                                          acGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTtctc                              <  2:3276146/93‑1 (MQ=255)

                                                          acGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGAACTTCTACACGCCGGACTtctc                              <  1:1393562/93‑1 (MQ=255)

                                                           cgcgGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTcc                             <  2:3786826/93‑1 (MQ=255)

                                                             cgGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTATACGCCGGACTGCTCCGc                           >  2:7010442/1‑93 (MQ=255)

                                                               gCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGcac                         <  2:921941/93‑1 (MQ=255)

                                                                cACGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGcaca                        <  2:1178049/93‑1 (MQ=255)

                                                                 aCGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAg                       <  2:1099370/93‑1 (MQ=255)

                                                                 aCGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAg                       <  1:1348177/93‑1 (MQ=255)

                                                                  cGGCCAGTTCCGCCCGTTCCTGCT‑G‑TGGGG‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGcaca                        >  1:3871350/1‑93 (MQ=255)

                                                                  cGGCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAgg                      <  2:2544643/93‑1 (MQ=255)

                                                                   ggCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAggg                     >  1:1081286/1‑93 (MQ=255)

                                                                   ggCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAggg                     <  2:5082289/93‑1 (MQ=255)

                                                                    gCCAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGc                    <  2:3203597/93‑1 (MQ=255)

                                                                     ccAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCCCAGGGCa                   >  1:2803574/1‑93 (MQ=255)

                                                                     ccAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACCGGGCa                   >  1:2803577/1‑93 (MQ=255)

                                                                      cAGTTCCGCCCGTTGCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCACGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCaa                  <  1:868813/93‑1 (MQ=255)

                                                                      cAGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCaa                  <  2:188344/93‑1 (MQ=255)

                                                                      cAGTTCCGCCCGTGCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCaa                  <  1:960018/93‑1 (MQ=255)

                                                                       aGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAg                 >  1:4777760/1‑93 (MQ=255)

                                                                       aGTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAg                 >  1:4318750/1‑93 (MQ=255)

                                                                        gTTCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGa                <  1:7253948/93‑1 (MQ=255)

                                                                         ttCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGCCTTCTCCGCACAGGGCAAGat               >  1:1719490/1‑93 (MQ=255)

                                                                         ttCCGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGat               >  1:636336/1‑93 (MQ=255)

                                                                            cGCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatcat            >  1:1549143/1‑93 (MQ=255)

                                                                             gCCCGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatcatc           <  2:2381493/93‑1 (MQ=255)

                                                                             gCCCGTTACTGCT‑G‑T‑‑‑GAGG‑CGGACATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatcatc           <  2:3884710/93‑1 (MQ=255)

                                                                              cccGTTCCTGCT‑G‑TGG‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGatcatc           <  2:2191638/93‑1 (MQ=255)

                                                                              cccGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCt          >  2:7721648/1‑93 (MQ=255)

                                                                               ccGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTa         >  1:512272/1‑93 (MQ=255)

                                                                               ccGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTa         <  2:6709729/93‑1 (MQ=255)

                                                                                cGTTCCTGCTGG‑G‑‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc        <  2:4485318/93‑1 (MQ=255)

                                                                                cGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc        <  1:7833926/93‑1 (MQ=255)

                                                                                cGTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTAc        <  2:1629140/93‑1 (MQ=255)

                                                                                 gTTCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACg       >  1:5203308/1‑93 (MQ=255)

                                                                                  ttCCTGCT‑G‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGc      >  2:5517833/1‑93 (MQ=255)

                                                                                     ctgct‑g‑TG‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTg   <  2:2321357/93‑1 (MQ=255)

                                                                                      tgct‑t‑T‑‑‑G‑GGTGGTCCATCCCGTTTTGCTTCGTCTGCGTGTTTACCTTCTACACGCCGGACTTCTCTGTACAGGGCTAGATCATCCACTCCttc  >  2:7248702/1‑91 (MQ=255)

CCGCGTGCTCGACGCCGTCACCGACCCGCTGATGGGGCTGCTGGTAGACCGCACCCGCACGCGGCACGGCCAGTTCCGCCCGTTCCTGCT‑G‑T‑‑‑G‑GG‑GGGCCATCCCGTTCGGCATCGTCTGCGTGCTGACCTTCTACACGCCGGACTTCTCCGCACAGGGCAAGATCATCTACGCCTGC  >  NC_000913/284774‑284951

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