BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,494,793	A→T	29.9%	P92P (CCA→CCT) ‡	ibpB →	heat shock chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,494,793	0	A	T	29.9%	146.0 / 77.5	117	P92P (CCA→CCT) ‡	ibpB	heat shock chaperone
Reads supporting (aligned to +/- strand): ref base A (34/48); new base T (17/18); total (51/66)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 119 total aligned reads displayed.

GTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGTTT  >  minE/2494725‑2494858

gTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCa                                                                   <  1:220605/69‑1 (MQ=255)

   aaaGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCAttt                                                                >  1:1986627/1‑69 (MQ=255)

   aaaGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCAttt                                                                >  1:870435/1‑69 (MQ=255)

   aaaGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCAttt                                                                >  1:1982735/1‑69 (MQ=255)

   aaaGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCAttt                                                                >  1:2082126/1‑69 (MQ=255)

    aacgCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:2287291/4‑69 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAATAAATGGCTGCATCAAGGGCTTATGAATCAGTCtttt                                                                >  1:318181/1‑68 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:1292014/1‑69 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:1871553/1‑69 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:3261174/1‑69 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:1767732/1‑69 (MQ=255)

    aaGGCACGCCAGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGTCTTTTa                                                               >  1:1749676/1‑69 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGTGCTTATGAATCAGCCATTTAGc                                                             >  1:2455675/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:2051386/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:1929642/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:1971425/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:1819120/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:2548000/1‑67 (MQ=255)

        cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGc                                                             >  1:678706/1‑67 (MQ=255)

              ggAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt                                                     <  1:2876480/69‑1 (MQ=255)

              ggAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt                                                     <  1:3154645/69‑1 (MQ=255)

              ggAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt                                                     <  1:675807/69‑1 (MQ=255)

              ggAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt                                                     <  1:1085025/69‑1 (MQ=255)

                agcagcCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTa                                                   <  1:3122648/69‑1 (MQ=255)

                agcagcCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTa                                                   <  1:2188071/69‑1 (MQ=255)

                 gcagcCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTAc                                                  <  1:2336644/69‑1 (MQ=255)

                           gaagaGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:2449165/69‑1 (MQ=255)

                           gaagaGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:2504033/69‑1 (MQ=255)

                           gaagaGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:2989272/69‑1 (MQ=255)

                           gaagaGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:342369/69‑1 (MQ=255)

                             agagAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaaaa                                      <  1:2582752/69‑1 (MQ=255)

                             agagAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaaaa                                      <  1:2486514/69‑1 (MQ=255)

                             agagAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaaaa                                      <  1:1998408/69‑1 (MQ=255)

                             agagAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaaaa                                      <  1:1524034/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:2005966/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:1118793/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:426702/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:3077079/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:2849848/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:1607418/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:2507401/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:1793146/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:2118025/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:1966911/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:1988830/69‑1 (MQ=255)

                                  aaaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:215762/69‑1 (MQ=255)

                                   aaaaTGGCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:2206189/68‑1 (MQ=255)

                                   aaaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:1406379/61‑1 (MQ=255)

                                   aaaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:2555138/69‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:743492/68‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:3011402/68‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:1722884/68‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:1646398/68‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:1514493/68‑1 (MQ=255)

                                    aaaTGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATACGGa                                <  1:681531/68‑1 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:91347/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:2833235/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:2808173/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:3278842/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:1880699/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:561347/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:1865704/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:2277341/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTAAGCTGGCTGAAAATATGGAAGtct                           >  1:155321/1‑69 (MQ=255)

                                        ggCTGCATCAAGGGCTTATGAATCAGTCTTTTACCCTGAGCTTTACGCTGGCTGAAAATATTGAAGtct                           >  1:1109997/1‑69 (MQ=255)

                                         gatGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtctc                          >  1:2623219/3‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAGTATGGAAGtctc                          >  1:1734632/1‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtctc                          >  1:3407242/1‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtctc                          >  1:3226373/1‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtctc                          >  1:2996810/1‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtctc                          >  1:2255173/1‑69 (MQ=255)

                                         gCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtc                            >  1:2957338/1‑67 (MQ=255)

                                           tGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCtt                                                     <  1:2168938/40‑1 (MQ=255)

                                             cATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGa                                <  1:3251061/59‑1 (MQ=255)

                                                     gCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCtt              >  1:308875/1‑69 (MQ=255)

                                                     gCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCtt              >  1:1099472/1‑69 (MQ=255)

                                                     gCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCtt              >  1:682116/1‑69 (MQ=255)

                                                     gCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCtt              >  1:608276/1‑69 (MQ=255)

                                                        tATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATgg                                 <  1:3024759/47‑1 (MQ=38)

                                                              tCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACgg     >  1:2816090/1‑69 (MQ=255)

                                                              tCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACg      <  1:2729262/68‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:2846180/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:903439/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:789361/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:1241229/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:2665652/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:2189428/69‑1 (MQ=255)

                                                                aGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   <  1:570035/69‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:418461/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:2356562/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:498324/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:1138188/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:2783555/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:2795081/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGttt  <  1:13730/68‑1 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   >  1:252352/1‑67 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   >  1:1012790/1‑67 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   >  1:1854323/1‑67 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   >  1:1414061/1‑67 (MQ=255)

                                                                  ccATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGtt   >  1:2872001/1‑67 (MQ=255)

GTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGTTT  >  minE/2494725‑2494858

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