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

Only 100 of 146 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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

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

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

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

  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:448388/1‑48 (MQ=255)

  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:2494755/1‑48 (MQ=255)

  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:3188842/1‑48 (MQ=255)

  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:1336222/1‑48 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:871814/1‑50 (MQ=255)

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

           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:3224563/1‑50 (MQ=255)

           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:125553/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                         aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     >  1:5116406/1‑50 (MQ=255)

                         aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     >  1:2595181/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:3702996/50‑1 (MQ=255)

                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:2568713/50‑1 (MQ=255)

                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:3598600/50‑1 (MQ=255)

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

                             aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 <  1:3247286/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  1:4240314‑M1/49‑4 (MQ=255)

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

                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  1:574192‑M1/50‑7 (MQ=255)

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

                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:3007151‑M1/49‑6 (MQ=255)

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

                                            tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:4408072‑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` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`