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

Only 100 of 149 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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

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

tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:3444358/50‑1 (MQ=255)

 aaaaGCGTGAAACAACTCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:6387726/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

               aaTCACTATTACAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               <  1:635670/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                         <  1:4837890/48‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                             aaagTAATCTATATTTCTGC‑GCATTTCCAGCCTCCTGTGTTGATTTCCaa                 >  1:2931995/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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