BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,525	(C)_6→7	coding (402/630 nt)	yqiE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,519	1	.	C	96.6%	724.6 / 16.2	238	coding (408/630 nt)	yqiE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (5/3); new base C (133/97); total (138/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 290 total aligned reads displayed.

CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACA  >  NC_000913/3175427‑3175606

 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑c                                                                                         <  1:451625/93‑1 (MQ=255)

 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑c                                                                                         >  1:7304437/1‑93 (MQ=255)

  aaTATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCATTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cc                                                                                        <  1:3900076/93‑1 (MQ=255)

   atatCTTCTTTTTCATCAGCCAGCCCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑ccc                                                                                       <  2:5736007/93‑1 (MQ=255)

    tatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccc                                                                                      >  1:5891315/1‑93 (MQ=255)

     atCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑ccccc                                                                                     >  2:528629/1‑93 (MQ=255)

     atCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑ccccc                                                                                     >  2:7344907/1‑93 (MQ=255)

      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccc                                                                                    >  2:2036651/1‑93 (MQ=255)

      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccc                                                                                    >  1:2121954/1‑93 (MQ=255)

      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGG‑G‑‑cccccc                                                                                    >  1:6718535/1‑93 (MQ=255)

        ttCTTTTTCATCAGCCAGACCGTGAATCCCGCTGGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccccg                                                                                  <  1:1556360/93‑3 (MQ=255)

        ttCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑cccccccg                                                                                  >  1:7042018/1‑91 (MQ=255)

             tttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCtt                                                                             >  1:3646459/1‑93 (MQ=255)

               tcatcaCCCAGGCCGTGAATCCCGCTTGCGGTCGTGGCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGc                                                                           <  2:3109070/93‑1 (MQ=255)

                  tcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAg                                                                        >  1:4551044/1‑93 (MQ=255)

                    aGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGa                                                                      >  2:1780060/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCGGGaa                                                                     >  2:1442823/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGaa                                                                     >  2:3205241/1‑93 (MQ=255)

                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGaa                                                                     >  1:4372362/1‑93 (MQ=255)

                      ccAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGaaa                                                                    <  1:4294057/93‑1 (MQ=255)

                       cAGACGGTGAATACCGATTGGGGTCGTGGCATCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAAc                                                                   <  2:4437948/93‑1 (MQ=255)

                       cAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGaaa                                                                    >  1:5011711/1‑93 (MQ=255)

                       cAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAAc                                                                   >  2:1719887/1‑93 (MQ=255)

                        aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACt                                                                  >  2:3492329/1‑93 (MQ=255)

                        aGACCGTGAATACCGCTTGCGGCCGTGGCGTCCACTTCGCTCACTCTAATTAACGAACGCTCACTGGG‑GT‑CCCCCC‑GGGCTTGCCAGGACACt                                                                  >  1:2404271/1‑93 (MQ=255)

                           ccGTGAATCCCGCTTGCGTTCGTGCCGTCCACTTCCCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTaa                                                               <  1:5963308/93‑1 (MQ=255)

                           ccGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTaa                                                               >  1:4659853/1‑93 (MQ=255)

                             gTGATTCCCGCTGGCGGAGGTGGCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAaca                                                             <  2:5140636/93‑1 (MQ=255)

                             gTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAac                                                              >  2:1346595/1‑93 (MQ=255)

                             gTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAaca                                                             >  2:2305459/1‑93 (MQ=255)

                              tGAATCCCGCTTGCGGTCGTGCCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAaca                                                             <  1:1182375/93‑1 (MQ=255)

                               gAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTGGCCAGTAAACTTAACGcc                                                           >  1:1122490/1‑93 (MQ=255)

                                 atcccGCTTGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCgg                                                         <  2:6400213/90‑1 (MQ=255)

                                 aTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCgg                                                         >  1:1691110/1‑93 (MQ=255)

                                   cccccTTGCGGTCGTGTCGCCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGtt                                                       <  1:4801911/89‑1 (MQ=255)

                                   cccGCTTGGGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGtt                                                       <  2:2575056/92‑1 (MQ=255)

                                   aCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGtt                                                       >  1:7734941/1‑93 (MQ=255)

                                    ccGCTTGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGttt                                                      <  1:3549400/93‑1 (MQ=255)

                                    ccGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGttt                                                      <  1:2746735/93‑1 (MQ=255)

                                     cGCTTGCGGTGGGGGCTCCCCCTGCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTg                                                     <  1:3360525/93‑1 (MQ=255)

                                     cGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGttt                                                      >  1:5264413/1‑93 (MQ=255)

                                     cGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTg                                                     >  2:4936881/1‑93 (MQ=255)

                                      gCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                    >  2:7544975/1‑93 (MQ=255)

                                      gCTTGCGGTCGGGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                    <  2:2436323/93‑1 (MQ=255)

                                      ccTTGCGGTGGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                    <  1:1768701/92‑1 (MQ=255)

                                       cTTGCGGTGGTGGCTTCCACTCCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                   <  2:4769874/93‑1 (MQ=255)

                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCTAGGAAACTTAACACCGGTTTGGt                                                   >  1:985842/1‑93 (MQ=255)

                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                   >  1:7109735/1‑93 (MQ=255)

                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                   >  2:4259733/1‑93 (MQ=255)

                                        ttGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTc                                                  <  1:5450773/93‑1 (MQ=255)

                                        ttGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTc                                                  >  1:6624980/1‑93 (MQ=255)

                                         tGCGGTCGTGGCGTCCCCTCCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTcc                                                 <  2:6404100/93‑1 (MQ=255)

                                          gCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                >  2:7688246/1‑93 (MQ=255)

                                          gCGGTCGTGGCGTCCACTCCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                <  2:4892035/93‑1 (MQ=255)

                                             gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTCACACCGGTTTGGTCCGttt                                             >  2:1844905/1‑93 (MQ=255)

                                             gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                             >  2:5892236/1‑93 (MQ=255)

                                             gTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                             >  1:3164347/1‑93 (MQ=255)

                                               cGTGCCCTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGa                                           <  2:1094803/93‑1 (MQ=255)

                                                gTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                          <  1:4466071/93‑1 (MQ=255)

                                                gTGGCGTCCACCTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                          <  1:2811290/93‑1 (MQ=255)

                                                 tgccGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACt                                         <  2:2203831/90‑1 (MQ=255)

                                                 tGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                          >  1:3026827/1‑93 (MQ=255)

                                                  ggCGTCCCCTCCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACt                                         <  1:867454/93‑1 (MQ=255)

                                                   gCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                       >  1:6283633/1‑93 (MQ=255)

                                                     gTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                     >  2:6541601/1‑93 (MQ=255)

                                                     ctCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                     <  2:3607964/92‑1 (MQ=255)

                                                      tCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAg                                    >  2:7851279/1‑93 (MQ=255)

                                                       ccACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                   >  2:7787896/1‑93 (MQ=255)

                                                       ccACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                   >  2:3710697/1‑93 (MQ=255)

                                                       ccACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                   >  2:5417904/1‑93 (MQ=255)

                                                         ccTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTc                                  <  1:2005343/92‑1 (MQ=255)

                                                         ccTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                 <  2:1285745/92‑1 (MQ=255)

                                                         ccTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                 <  2:6213136/92‑1 (MQ=255)

                                                         aCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTc                                  >  1:4135538/1‑93 (MQ=255)

                                                         aCTTCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                 >  1:3116471/1‑93 (MQ=255)

                                                           ttCGCCCACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                               >  2:343410/1‑93 (MQ=255)

                                                            tcccccACCATAATTGACGAACGCTCACTGGTCC‑‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                              <  1:5143836/90‑1 (MQ=255)

                                                             cccccACCATAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGcc                             <  2:5545925/91‑1 (MQ=255)

                                                                  accaTAATTGACGAACGCTCACTGGG‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTc                        >  1:5635739/1‑93 (MQ=255)

                                                                   ccaTAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCa                       >  1:7220375/1‑93 (MQ=255)

                                                                    caTAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCaa                      >  1:2530375/1‑93 (MQ=255)

                                                                    caTAATTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCaa                      >  2:4144291/1‑93 (MQ=255)

                                                                       aaTTGACGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                   >  1:7068253/1‑93 (MQ=255)

                                                                          tggcGAACGCTCACTGGTCC‑‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCtt                <  2:468870/90‑1 (MQ=255)

                                                                          tGACGACCGCTCACTGGTCC‑‑CCCCCC‑GGGCTTCCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCtt                <  2:859548/93‑1 (MQ=255)

                                                                           gACGAACGCTCACTGGG‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc               <  1:580260/93‑1 (MQ=255)

                                                                           gACGAACGCTCACGGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc               <  1:4086210/93‑1 (MQ=255)

                                                                            aCGAACGCTCACTGGT‑GCCCCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc               <  1:3999551/93‑1 (MQ=255)

                                                                             cGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgc             <  2:4881869/93‑1 (MQ=255)

                                                                             cGAACGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgc             <  2:2527605/93‑1 (MQ=255)

                                                                              gAACGCTCACTGGT‑G‑‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           >  2:1982936/1‑93 (MQ=255)

                                                                               accgctcgctGGT‑G‑‑CCCCCCGGGGCTTGCGAGGAAACTTAGCGCGGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           <  1:3316826/91‑1 (MQ=255)

                                                                               aaCGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           <  2:3909287/93‑1 (MQ=255)

                                                                               aaCGCTCACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa           <  1:2900033/93‑1 (MQ=255)

                                                                                     cACTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCa     >  2:3316674/1‑93 (MQ=255)

                                                                                      ccTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCaa    <  2:6679688/92‑1 (MQ=255)

                                                                                      aCTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGCCGGGCac    >  2:3138223/1‑92 (MQ=255)

                                                                                       cTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc   >  2:2971958/1‑93 (MQ=255)

                                                                                       cTGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc   >  1:5554/1‑93 (MQ=255)

                                                                                        tGGT‑GC‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACa  >  2:3079502/1‑93 (MQ=255)

CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGT‑G‑‑CCCCCC‑GGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACA  >  NC_000913/3175427‑3175606

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