BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	547,694	A→G	E38E (GAA→GAG)	b4507 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	547,694	0	A	G	99.2%	328.5 / ‑5.8	127	E38E (GAA→GAG)	b4507	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (71/55); minor base T (0/1); total (71/56)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 129 total aligned reads displayed.

ACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGGCGTGCGTG  >  NC_000913/547622‑547769

 cAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                                         <  2:383792/76‑1 (MQ=255)

 cAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCa                                                                         <  2:2926886/76‑1 (MQ=255)

  aaaTTCGTCGCGCTCGTCCACACTGGCTGGATTTGCAACCGGCTTCTTCCCTTTTCAGCCAACTAAACGCGGGCaa                                                                        >  2:4866358/1‑76 (MQ=255)

  aaaTTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGTACTAAACGAGGGCaa                                                                        >  2:1468870/1‑76 (MQ=255)

  aaaTTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCaa                                                                        <  2:4403351/76‑1 (MQ=255)

     ttCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAac                                                                     <  1:2955859/76‑1 (MQ=255)

      tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCATTTATCAGCGAACTAAACGAGGGCAAAacc                                                                    >  2:2739488/1‑75 (MQ=255)

      tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                                    >  1:1100198/1‑76 (MQ=255)

      tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                                    <  1:3571848/76‑1 (MQ=255)

      tcgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAaca                                                                    >  1:467635/1‑76 (MQ=255)

       cgtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAacac                                                                   <  2:499033/76‑1 (MQ=255)

        gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGACCTAAACGAGGGCAAAACACt                                                                  <  1:1099195/76‑1 (MQ=255)

        gtcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACt                                                                  <  2:384127/76‑1 (MQ=255)

         tcgCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTg                                                                 >  1:4367061/1‑76 (MQ=255)

           gcgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCt                                                               <  1:2829522/76‑1 (MQ=255)

            cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCCGCGAACTAAACGAGGGCAAAACACTGCtc                                                              >  2:2245543/1‑75 (MQ=255)

            cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                              <  2:4171789/76‑1 (MQ=255)

            cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCtt                                                              <  1:2971870/76‑1 (MQ=255)

            cgcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACCGcct                                                              >  2:658811/1‑74 (MQ=255)

             gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTc                                                             >  1:1908309/1‑76 (MQ=255)

             gcTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACAGCTTc                                                             >  2:4067775/1‑76 (MQ=255)

             ccTCGCCCACACTGGCACGATGTGCAAGCGGCTTCTTCACATATCCGCGAACTACGCTAGGGCAAACCACTGCTTc                                                             <  1:2421280/75‑1 (MQ=255)

              cTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCa                                                            <  2:1471364/76‑1 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTccc                                                           >  2:3955022/1‑74 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                           >  1:2791651/1‑76 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                           <  1:2618978/76‑1 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                           >  1:905648/1‑76 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                           >  1:424114/1‑76 (MQ=255)

               tCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCAc                                                           <  1:4480457/76‑1 (MQ=255)

                cGTCCACACTGGCGGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACg                                                          <  2:348383/76‑1 (MQ=255)

                 ttCCACACTGGCTGGATGTGCAACCGTCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                                         <  1:1463287/75‑1 (MQ=255)

                 gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCCCTTATCAGCCAACTAAACGAGGGCAAAACACTGCTTCACGc                                                         >  2:4106751/1‑76 (MQ=255)

                 gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                                         >  2:1951064/1‑76 (MQ=255)

                 gTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGc                                                         >  2:3614261/1‑76 (MQ=255)

                  tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                                        >  1:3178384/1‑76 (MQ=255)

                  tCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGcc                                                        <  2:4235860/76‑1 (MQ=255)

                   ccACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCg                                                       >  2:659167/1‑76 (MQ=255)

                    cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                                      <  2:3786272/76‑1 (MQ=255)

                    cacaCTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCgg                                                      <  2:3758259/76‑1 (MQ=255)

                     acacTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCggg                                                     >  1:2052860/1‑76 (MQ=255)

                       acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                                   >  2:2981515/1‑76 (MQ=255)

                       acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                                   <  2:2045377/76‑1 (MQ=255)

                       acTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                                   <  1:3453502/76‑1 (MQ=255)

                       acTGGCGGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgcc                                                   <  2:2970896/76‑1 (MQ=255)

                          ggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                                <  2:1916457/76‑1 (MQ=255)

                          ggctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGgccgcc                                                >  2:1003491/1‑76 (MQ=255)

                            ctggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaa                                              >  2:2526292/1‑76 (MQ=255)

                             tggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAAt                                             <  1:3021931/76‑1 (MQ=255)

                              ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaacg                                            >  2:2028284/1‑74 (MQ=255)

                              ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCaacg                                            >  2:482856/1‑74 (MQ=255)

                              ggATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATg                                            <  1:2317637/76‑1 (MQ=255)

                               gATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgc                                           >  2:3369148/1‑76 (MQ=255)

                               gATGTGCAACCGGCTTCTTCACTGATCAGCGAACTAAACGATGGCAAACCACTGCTTCACGCCGGGCCGCCAATgc                                           <  1:4536643/76‑1 (MQ=255)

                                 tgtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcgc                                         >  1:2060238/1‑76 (MQ=255)

                                 ggtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATgcgc                                         <  2:4552957/75‑1 (MQ=255)

                                  gtgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCt                                        >  2:3253170/1‑76 (MQ=255)

                                   tgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                                       >  2:4348410/1‑76 (MQ=255)

                                   tgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTg                                       <  1:2551414/76‑1 (MQ=255)

                                   tgCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCCGGCCGCCAATGCGCTg                                       >  2:4490871/1‑76 (MQ=255)

                                    gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                                      <  2:146281/76‑1 (MQ=255)

                                    gCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTgg                                      <  1:523755/76‑1 (MQ=255)

                                     cAACCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGc                                     <  1:4184554/76‑1 (MQ=255)

                                     cAACCGGCTTCTTCACCGATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGGCGGGCCGCCAATGCGCTGGc                                     <  2:2746851/76‑1 (MQ=255)

                                      aaCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCa                                    >  1:1761745/1‑76 (MQ=255)

                                      aaCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGACCGCCAATGCGCTGGCa                                    >  2:2487327/1‑76 (MQ=255)

                                       aCCGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCCATGCGCTGGCAg                                   >  1:2139706/1‑76 (MQ=255)

                                        ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCCgg                                  >  2:697401/1‑76 (MQ=255)

                                        ccGGCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAgg                                  <  1:1778451/76‑1 (MQ=255)

                                          ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCACTGCGCCGGCAGgag                                >  1:714940/1‑76 (MQ=255)

                                          ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                                >  1:110688/1‑76 (MQ=255)

                                          ggCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgag                                >  1:2006879/1‑76 (MQ=255)

                                           gCTTCTTCACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGgaga                               >  1:4040556/1‑76 (MQ=255)

                                            nttcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                              >  1:2499512/2‑76 (MQ=255)

                                            gttcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                              <  2:2435122/75‑1 (MQ=255)

                                            gttcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGAt                              <  2:4604183/75‑1 (MQ=255)

                                              tcttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGa                            >  1:3275646/1‑76 (MQ=255)

                                               cttcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGAc                           >  2:1626419/1‑76 (MQ=255)

                                                 tcACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCg                         <  1:1208221/76‑1 (MQ=255)

                                                  cACTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCgg                        >  2:3309656/1‑76 (MQ=255)

                                                   aCTTATCAGCGAACTAAACGAGGGCAACACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCggc                       >  2:2203660/1‑75 (MQ=255)

                                                   aCTTATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGa                       <  1:4852866/76‑1 (MQ=255)

                                                     ttATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAcc                     <  2:3059397/76‑1 (MQ=255)

                                                      tATCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGAccc                    >  1:1053040/1‑76 (MQ=255)

                                                       aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa                   >  1:104585/1‑76 (MQ=255)

                                                       aTCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCa                   >  1:697629/1‑76 (MQ=255)

                                                        tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATgacccgacccct                  >  1:90835/1‑74 (MQ=255)

                                                        tCAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCAt                  <  1:659167/76‑1 (MQ=255)

                                                         cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATg                 >  2:3316541/1‑76 (MQ=255)

                                                         cAGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATg                 <  1:1011207/76‑1 (MQ=255)

                                                          aGCGAACTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCCTGa                >  2:2327295/1‑76 (MQ=255)

                                                              aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAACgg            >  1:1395706/1‑76 (MQ=255)

                                                              aaCTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAgg            <  1:2307045/76‑1 (MQ=255)

                                                                cTAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCGCGCAGGAGATGACCGGACCCATGAAAgggg          >  2:3023213/1‑76 (MQ=255)

                                                                 tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGCAAGGGgc         >  2:4475547/1‑76 (MQ=255)

                                                                 tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGgc         >  1:750069/1‑76 (MQ=255)

                                                                 tAAACGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGgc         <  1:1666398/76‑1 (MQ=255)

                                                                  aaaCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGgcg        >  1:3315424/1‑76 (MQ=255)

                                                                    aCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGgcgtg      <  1:4470184/76‑1 (MQ=255)

                                                                    aCGAGGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGgcgtg      >  1:1092200/1‑76 (MQ=255)

                                                                        aGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGGcgtgcgtg  <  1:1206809/76‑1 (MQ=255)

ACAAATTCGTCGCGCTCGTCCACACTGGCTGGATGTGCAACCGGCTTCTTCACTTATCAGCGAACTAAACGAAGGCAAAACACTGCTTCACGCCGGGCCGCCAATGCGCTGGCAGGAGATGACCGGACCCATGAAAGGGGCGTGCGTG  >  NC_000913/547622‑547769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`