BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,494,791	C→T	34.5%	P92S (CCA→TCA) ‡	ibpB →	heat shock chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,494,791	0	C	T	34.5%	78.6 / 58.3	104	P92S (CCA→TCA) ‡	ibpB	heat shock chaperone
Reads supporting (aligned to +/- strand): ref base C (27/41); new base T (18/18); total (45/59)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGTTT  >  minE/2494724‑2494858

cGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGcc                                                                    >  1:1162132/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

         cACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATCAATCAGCCATTTAGc                                                             >  1:3045955/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                            gaagaGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGaa                                        <  1:2449165/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:1988830/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                         ggCTGCATCAAGGGCTTATGAATCAGTCTTTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGtct                           >  1:2517534/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:2996810/1‑69 (MQ=255)

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

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

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

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

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

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

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

                                                      gCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCtt              >  1:308875/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

CGTAAAAGGCACGCCGGAGCAGCCAAAAGAAGAGAAAAAATGGCTGCATCAAGGGCTTATGAATCAGCCATTTAGCCTGAGCTTTACGCTGGCTGAAAATATGGAAGTCTCTGGCGCAACCTTCGTAAACGGTTT  >  minE/2494724‑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: `‑`