BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,777,543	G→A	E31K (GAA→AAA)	ypjF →	CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,777,543	0	G	A	100.0%	975.0 / NA	235	E31K (GAA→AAA)	ypjF	CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (124/111); total (124/111)

Only 100 of 239 total aligned reads displayed.

GAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGT  >  NC_000913/2777445‑2777640

gagCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGaaa                                                                                                 >  2:3589642/1‑101 (MQ=255)

  gcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAaaca                                                                                               >  1:3589643/1‑101 (MQ=255)

   caagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAaacaa                                                                                              <  1:3589269/101‑1 (MQ=255)

    aagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAaacaac                                                                                             >  2:3589646/1‑101 (MQ=255)

     agcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAacaaca                                                                                            <  1:3589067/101‑1 (MQ=255)

     agcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAacaaca                                                                                            <  2:3589362/101‑1 (MQ=255)

     agcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAacaaca                                                                                            >  1:3589650/1‑101 (MQ=255)

     agcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAacaaca                                                                                            <  1:3589259/101‑1 (MQ=255)

      gcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAacac                                                                                           <  1:3589279/101‑1 (MQ=255)

       caTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACt                                                                                          >  2:3589651/1‑101 (MQ=255)

        aTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTa                                                                                         >  2:3589653/1‑101 (MQ=255)

         tGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTAt                                                                                        >  1:3589656/1‑101 (MQ=255)

           aaCACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATgg                                                                                      >  1:3589658/1‑101 (MQ=255)

             cacTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGcc                                                                                    >  2:3589659/1‑101 (MQ=255)

              acTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCt                                                                                   <  1:3589240/101‑1 (MQ=255)

                tctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGa                                                                                 >  1:3589663/1‑101 (MQ=255)

                 ctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGac                                                                                <  2:3589260/101‑1 (MQ=255)

                  tACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGaca                                                                               >  1:3589665/1‑101 (MQ=255)

                    ccTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAacac                                                                                           <  1:3589386/87‑1 (MQ=255)

                     cTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGcc                                                                                    <  1:3589226/93‑1 (MQ=255)

                     cTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCt                                                                                   >  2:3589673/1‑94 (MQ=255)

                     cTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTg                                                                            >  2:3589670/1‑101 (MQ=255)

                     cTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTg                                                                            >  1:3589671/1‑101 (MQ=255)

                       gCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGaa                                                                          <  1:3589404/101‑1 (MQ=255)

                        cTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAAc                                                                         <  1:3589108/101‑1 (MQ=255)

                        cTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAAc                                                                         >  2:3589678/1‑101 (MQ=255)

                        cTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAAc                                                                         <  1:3589342/101‑1 (MQ=255)

                        cTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAAc                                                                         >  1:3589676/1‑101 (MQ=255)

                         tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACg                                                                        >  2:3589681/1‑101 (MQ=255)

                           cAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGac                                                                      >  1:3589683/1‑101 (MQ=255)

                            aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGaca                                                                     >  2:3589684/1‑101 (MQ=255)

                            aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGaca                                                                     >  2:3589688/1‑101 (MQ=255)

                             aTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGacac                                                                    >  2:3589690/1‑101 (MQ=255)

                                tCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGcc                                                                 <  1:3589454/101‑1 (MQ=255)

                                tCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGcc                                                                 >  2:3589693/1‑101 (MQ=255)

                                tCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGcc                                                                 >  1:3589694/1‑101 (MQ=255)

                                  gCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCt                                                                                   >  2:3589697/1‑81 (MQ=255)

                                   cAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGtt                                                              >  2:3589699/1‑101 (MQ=255)

                                    aGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTc                                                             <  1:3589292/101‑1 (MQ=255)

                                      gcggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAg                                                           >  1:3589701/1‑101 (MQ=255)

                                      gcggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAg                                                           <  2:3589394/101‑1 (MQ=255)

                                       cggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGt                                                          >  2:3589705/1‑101 (MQ=255)

                                        ggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtg                                                         <  2:3589282/101‑1 (MQ=255)

                                        ggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtg                                                         >  2:3589708/1‑101 (MQ=255)

                                          cgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                       >  1:3589709/1‑101 (MQ=255)

                                               cccTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtg                                                         <  2:3589275/94‑1 (MQ=255)

                                                cctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACt                                                 <  2:3589177/101‑1 (MQ=255)

                                                cctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACt                                                 >  1:3589712/1‑101 (MQ=255)

                                                 ctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtga                                                        >  2:3589713/1‑93 (MQ=255)

                                                 cggcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACg                                                                   <  2:3589069/80‑1 (MQ=255)

                                                   gcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                              <  1:3589122/101‑1 (MQ=255)

                                                     ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGa                                                                                 <  1:3589328/64‑1 (MQ=255)

                                                     ctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGttat                                            >  2:3589717/1‑101 (MQ=255)

                                                      tgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGttatt                                           <  1:3589247/101‑1 (MQ=255)

                                                        tCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTaa                                         <  2:3589241/101‑1 (MQ=255)

                                                        tCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTaa                                         <  2:3589541/101‑1 (MQ=255)

                                                         cGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAg                                        >  2:3589719/1‑101 (MQ=255)

                                                          gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAgg                                       >  2:3589723/1‑101 (MQ=255)

                                                          gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAgg                                       >  1:3589722/1‑101 (MQ=255)

                                                           ccAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGa                                      <  2:3589337/101‑1 (MQ=255)

                                                            cAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGaa                                     >  1:3589727/1‑101 (MQ=255)

                                                             aGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAAc                                    <  1:3589301/101‑1 (MQ=255)

                                                             aGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAAc                                    <  2:3589405/101‑1 (MQ=255)

                                                             aGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAAc                                    >  2:3589728/1‑101 (MQ=255)

                                                              gTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACa                                   <  2:3589173/101‑1 (MQ=255)

                                                              gTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACa                                   <  2:3589268/101‑1 (MQ=255)

                                                              gTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACa                                   <  2:3589451/101‑1 (MQ=255)

                                                               tGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACat                                  >  2:3589730/1‑101 (MQ=255)

                                                               tGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACat                                  <  1:3589544/101‑1 (MQ=255)

                                                                ggCTGTCTGGCAAATGCAACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACata                                 >  2:3589734/1‑101 (MQ=255)

                                                                  ctttCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATc                               >  1:3589736/4‑101 (MQ=255)

                                                                  cTGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATc                               <  2:3589370/101‑1 (MQ=255)

                                                                   tGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCg                              >  2:3589737/1‑101 (MQ=255)

                                                                   tGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCg                              <  2:3589412/101‑1 (MQ=255)

                                                                   tGTCTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCg                              <  1:3589363/101‑1 (MQ=255)

                                                                      cTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATg                           <  2:3589443/101‑1 (MQ=255)

                                                                      cTGGCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATg                           <  1:3589323/101‑1 (MQ=255)

                                                                         gCAAATGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                        >  1:3589740/1‑101 (MQ=255)

                                                                            aaTGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTa                     <  2:3589315/101‑1 (MQ=255)

                                                                            aaTGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTa                     >  2:3589743/1‑101 (MQ=255)

                                                                             aTGCTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTAt                    <  2:3589252/101‑1 (MQ=255)

                                                                                cTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCAc                 <  1:3589395/101‑1 (MQ=255)

                                                                                cTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCAc                 <  1:3589205/101‑1 (MQ=255)

                                                                                cTACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCAc                 >  2:3589745/1‑101 (MQ=255)

                                                                                 tACTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCAct                >  1:3589747/1‑101 (MQ=255)

                                                                                  aCTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCActc               <  1:3589376/101‑1 (MQ=255)

                                                                                  aCTGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCActc               >  2:3589750/1‑101 (MQ=255)

                                                                                    tGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg             <  1:3589596/101‑1 (MQ=255)

                                                                                    tGACACGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg             <  2:3589332/101‑1 (MQ=255)

                                                                                       cacGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGcc          <  2:3589203/101‑1 (MQ=255)

                                                                                       cacGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGcc          >  2:3589753/1‑101 (MQ=255)

                                                                                       cacGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGcc          >  1:3589754/1‑101 (MQ=255)

                                                                                         cGCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGa        >  1:3589758/1‑101 (MQ=255)

                                                                                          gCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGAt       <  2:3589222/101‑1 (MQ=255)

                                                                                          gCCTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGAt       >  2:3589759/1‑101 (MQ=255)

                                                                                           ccTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg      <  1:3589423/101‑1 (MQ=255)

                                                                                           ccTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg      >  1:3589761/1‑101 (MQ=255)

                                                                                           ccTGCTGAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg      >  1:3589763/1‑101 (MQ=255)

                                                                                            ctgctgAAACAACACTATGGCCTGACACTTAACGACACGCCGTTCATTGATGAAACTGTTATTaa                                         <  1:3589192/65‑1 (MQ=255)

                                                                                             tgctgAAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCa    >  1:3589768/1‑101 (MQ=255)

GAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGT  >  NC_000913/2777445‑2777640

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