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	.	81.6%	352.9 / 107.4	114	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand): ref base G (16/5); new base . (45/48); total (61/53)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACC  >  NC_000913/3720548‑3720643

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

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

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

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

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

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

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

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

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

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCa                                        <  2:4090813/49‑1 (MQ=255)

         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAg                                       >  1:3320556/1‑50 (MQ=255)

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

          aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGc                                      >  2:201622/1‑50 (MQ=255)

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

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

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

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

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

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

                tCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCtg                                >  2:4990616/1‑50 (MQ=255)

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

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

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

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

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

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

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

                 cACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCtgt                               >  2:2528276/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

                      ttAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGa                          >  1:4841230/1‑50 (MQ=255)

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

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

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

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

                        aaagaaagTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGAtt                        >  2:1540402/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

                           gaaagTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTc                      >  2:1024399/1‑49 (MQ=255)

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

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

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

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

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

                             aagTAATCTATAATTCTAC‑GCATTCCCGGCCTCCTGTGTTGATTTCCAAc                  <  2:3569396/50‑1 (MQ=255)

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

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

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

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

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

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

                                  aTCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg               <  2:211433/49‑1 (MQ=255)

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

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

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

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

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

                                    cTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt            >  2:2795519/1‑50 (MQ=255)

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

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

                                     tataTTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATg           >  2:1391686/1‑50 (MQ=255)

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

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

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

                                       tatTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt            <  2:5206807/47‑1 (MQ=255)

                                       tatTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt            <  1:3279237/47‑1 (MQ=255)

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

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

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

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

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

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

                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa          <  2:3196555‑M1/48‑3 (MQ=255)

                                        atTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTAc        >  2:250427/1‑50 (MQ=255)

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

                                         tttCTGC‑GCATTCCCAGCGGCCTGTGTTGATTTCCAACGAGTATaactt       <  2:2393615‑M1/49‑6 (MQ=255)

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

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

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

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

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

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

                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg     <  2:1936780‑M1/49‑8 (MQ=255)

                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg     <  2:3559026‑M1/49‑8 (MQ=255)

                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg     <  2:1682552‑M1/49‑8 (MQ=255)

                                               cGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCAcc  >  1:1141364/1‑49 (MQ=255)

                                               cGGCATTCCCAGCCTACTGTGTTGATTTCCAACGAGTATGTACTGCAcc  >  2:2642024/1‑49 (MQ=255)

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACC  >  NC_000913/3720548‑3720643

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