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%	1063.8 / NA	214	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (125/89); total (125/89)

Only 100 of 231 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:3035309/1‑50 (MQ=255)

tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:190431/1‑50 (MQ=255)

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

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

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

  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:2718698/50‑3 (MQ=255)

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

     gCGTGAAACAAATCGCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         <  2:2316324/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATGTCTGC‑GCATTccc                                       <  1:3424747/50‑1 (MQ=255)

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

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

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

         gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     >  1:4118888/1‑50 (MQ=255)

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

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

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

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

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

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

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

             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:2561181/1‑50 (MQ=255)

             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:555827/1‑50 (MQ=255)

             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  2:2021162/1‑50 (MQ=255)

             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:2103918/1‑50 (MQ=255)

             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  2:761642/1‑50 (MQ=255)

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

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

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

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

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

                aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  2:2700159/50‑1 (MQ=255)

                aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  2:1009823/50‑1 (MQ=255)

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

                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             >  2:2917971/1‑50 (MQ=255)

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

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

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

                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            <  2:3740938/49‑1 (MQ=255)

                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           >  2:719944/1‑50 (MQ=255)

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

                    cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          <  1:1249405/50‑1 (MQ=255)

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

                     tATTAAAGAAAGTAATCTATATTTTTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:4512814/1‑49 (MQ=38)

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

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

                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       <  2:4107407/50‑1 (MQ=255)

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

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

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

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

                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:170575/50‑1 (MQ=255)

                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:3477957/1‑50 (MQ=255)

                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  2:4022066/50‑1 (MQ=255)

                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:2340217/50‑1 (MQ=255)

                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  2:38820/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  2:658389/1‑48 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                       atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGCTTTCCAACGAGTATaac       >  2:3945939‑M1/1‑47 (MQ=255)

                                       atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:1800290‑M1/1‑47 (MQ=255)

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

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

                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  2:464988‑M1/50‑5 (MQ=255)

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

                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  2:3560511‑M1/49‑4 (MQ=255)

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

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

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

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

                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  2:1385227‑M1/50‑7 (MQ=255)

                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  2:3445738‑M1/49‑6 (MQ=255)

                                          tttCTGC‑GCATTCCCAGCCTCCTGTGCTGATTTCCAACGAGTATaacttt    >  2:1174440‑M1/1‑44 (MQ=255)

                                            tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:2238163‑M1/1‑42 (MQ=255)

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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