BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,720,597	Δ1 bp	coding (4/153 nt)	hokA ←	small toxic polypeptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,596	0	G	.	100.0%	517.1 / NA	104	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (59/45); total (59/45)

Only 100 of 106 total aligned reads displayed.

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640

aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                              >  1:2613275/1‑49 (MQ=255)

 aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                             <  1:44618/49‑2 (MQ=255)

    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        >  1:2269789/1‑50 (MQ=255)

    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        <  1:1267891/50‑1 (MQ=255)

    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        <  1:5600322/50‑1 (MQ=255)

    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        <  1:3982506/50‑1 (MQ=255)

    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        <  1:7727205/50‑1 (MQ=255)

     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       <  1:8130722/50‑1 (MQ=255)

     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:780368/1‑50 (MQ=255)

      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:225826/50‑1 (MQ=255)

      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:1478505/1‑49 (MQ=39)

      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:3358066/1‑49 (MQ=39)

       tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:7079736/1‑48 (MQ=255)

       tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                     <  1:5104685/50‑1 (MQ=255)

       tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                     <  1:3195104/50‑1 (MQ=255)

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                    <  1:5234234/49‑1 (MQ=255)

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   <  1:6531555/50‑1 (MQ=255)

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  1:2701623/1‑50 (MQ=255)

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  1:949794/1‑50 (MQ=255)

          aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                  <  1:3227522/50‑1 (MQ=255)

           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 >  1:6142208/1‑50 (MQ=255)

           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 <  1:1777515/50‑1 (MQ=255)

           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 <  1:4275377/50‑1 (MQ=255)

           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 <  1:4352115/50‑1 (MQ=255)

            cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                <  1:5820819/50‑1 (MQ=255)

            cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                <  1:5746471/50‑1 (MQ=255)

             aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                               >  1:5376475/1‑50 (MQ=255)

              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              >  1:6992606/1‑50 (MQ=255)

              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:4414220/50‑1 (MQ=255)

              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:3584108/50‑1 (MQ=255)

              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:8283653/50‑1 (MQ=255)

               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                             >  1:2925302/1‑50 (MQ=255)

                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:6868771/50‑1 (MQ=255)

                 cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           >  1:5430844/1‑50 (MQ=255)

                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtgg                          >  1:568942/1‑49 (MQ=255)

                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          <  1:6568811/50‑1 (MQ=255)

                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:5184795/1‑50 (MQ=255)

                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:76259/1‑50 (MQ=255)

                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:916100/1‑50 (MQ=255)

                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:2212099/1‑50 (MQ=255)

                    tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:5120005/1‑49 (MQ=255)

                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                       >  1:2155330/1‑50 (MQ=255)

                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                       >  1:1340503/1‑50 (MQ=255)

                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                       >  1:1840571/1‑50 (MQ=255)

                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:6195006/1‑50 (MQ=255)

                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:2033097/1‑50 (MQ=255)

                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:2489955/1‑50 (MQ=255)

                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:2756575/1‑50 (MQ=255)

                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:1783774/1‑50 (MQ=255)

                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:1070963/1‑50 (MQ=255)

                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:6384707/1‑50 (MQ=255)

                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     <  1:5717473/49‑1 (MQ=255)

                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    <  1:6670817/49‑1 (MQ=255)

                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:8000401/1‑49 (MQ=255)

                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  >  1:6475016/1‑49 (MQ=255)

                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                >  1:3715516/1‑50 (MQ=255)

                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:4697309/50‑1 (MQ=255)

                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:4544738/50‑1 (MQ=255)

                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:3208770/50‑1 (MQ=255)

                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:7406404/1‑50 (MQ=255)

                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:4805253/1‑49 (MQ=255)

                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:1496485/1‑49 (MQ=255)

                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:1960476/1‑49 (MQ=255)

                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg              <  1:5772335/50‑1 (MQ=255)

                               gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa             <  1:1168211/50‑1 (MQ=255)

                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:5769088/50‑1 (MQ=255)

                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:875542/50‑1 (MQ=255)

                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:1461617/1‑50 (MQ=255)

                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:8067839/1‑50 (MQ=255)

                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:7514240/1‑50 (MQ=255)

                                 aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:3653935/1‑49 (MQ=255)

                                 aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:2928114/1‑49 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:7444955/49‑1 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:8097445/49‑1 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:598913/49‑1 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  1:939890/50‑1 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  1:918003/50‑1 (MQ=255)

                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:362466/1‑50 (MQ=255)

                                   tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:5055848/1‑49 (MQ=255)

                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:4485490/1‑49 (MQ=255)

                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:8263686/1‑49 (MQ=255)

                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  1:3347162‑M1/1‑48 (MQ=255)

                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  1:1917461‑M1/1‑48 (MQ=255)

                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  1:3772080‑M1/1‑48 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  1:4021362/46‑1 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  1:8033979/46‑1 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     <  1:1761994‑M1/50‑5 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     >  1:2398184‑M1/1‑46 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     >  1:7287363‑M1/1‑46 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       <  1:3336421‑M1/48‑3 (MQ=255)

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa        <  1:167421‑M1/47‑2 (MQ=255)

                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:3835855‑M1/1‑45 (MQ=255)

                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:1558500‑M1/1‑45 (MQ=255)

                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:7517250‑M1/1‑45 (MQ=255)

                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:3961924‑M1/1‑44 (MQ=255)

                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:718581‑M1/1‑44 (MQ=255)

                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:5183306‑M1/1‑43 (MQ=255)

                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:7785817‑M1/1‑43 (MQ=255)

                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:8175203‑M1/49‑8 (MQ=255)

                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:1310560‑M1/49‑8 (MQ=255)

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640

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