BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	229,931	C→T	32.6%	R314R (CGC→CGT)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	229,931	0	C	T	32.6%	88.1 / 51.0	101	R314R (CGC→CGT)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base C (49/19); new base T (13/20); total (62/39)
Fisher's exact test for biased strand distribution p-value = 2.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 108 total aligned reads displayed.

GCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGCGATGCATTA  >  minE/229864‑229994

gCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCa                                                                <  1:1957667/69‑1 (MQ=255)

gCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCa                                                                <  1:1283204/69‑1 (MQ=255)

      cccTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCa                                                          <  1:1013512/69‑1 (MQ=255)

      cccTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCa                                                          <  1:739218/69‑1 (MQ=255)

      cccTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCa                                                          <  1:612295/69‑1 (MQ=255)

      cccTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCa                                                          <  1:548029/69‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGTCTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:1958129/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAATTCATCCGCATTTGCATcc                                                       >  1:618196/1‑69 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:2390783/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:833873/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:2664258/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:2580907/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:1045570/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGc                                                                 <  1:1214878/59‑1 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAAcc                                                       >  1:468634/1‑69 (MQ=255)

         tGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAAcc                                                       >  1:2764544/1‑69 (MQ=255)

          gTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAAcc                                                       >  1:2307403/1‑68 (MQ=255)

          gTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAAcc                                                       >  1:1029455/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAAcc                                                       >  1:25243/1‑67 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1450787/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1463848/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:93106/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1006629/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:847992/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:3407384/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:2305963/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1734264/1‑68 (MQ=255)

           ttGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1074896/1‑68 (MQ=255)

           ttGCCGAAAGGCATT‑AAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:1151146/1‑67 (MQ=255)

           ttGCCGAAAGGAATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGTATTTGCAACCt                                                      >  1:435326/1‑68 (MQ=255)

                    ggCATTAAAAGCGTGACTGGCAATTTCTCGCGTGCTGAAGTCATCCGCATTTGCa                                                          <  1:3052722/55‑1 (MQ=255)

                          aaaaGCGTGGCTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGatat                                       >  1:2442064/1‑68 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAAGCTCGAAGGCCGCGATATc                                      >  1:686164/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGa                                          >  1:3103918/1‑65 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1856398/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1326663/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:229143/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:2929776/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1420052/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:2115870/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:86583/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:2687956/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1163906/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:729236/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:630970/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1946301/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:3223386/1‑69 (MQ=255)

                          aaaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:1790772/1‑69 (MQ=255)

                           aaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCg                                     >  1:201683/1‑69 (MQ=255)

                           aaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCg                                     >  1:2043941/1‑69 (MQ=255)

                           aaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCg                                     >  1:2316302/1‑69 (MQ=255)

                           aaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCg                                     >  1:1991323/1‑69 (MQ=255)

                           aaaGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCg                                     >  1:1348553/1‑69 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACGTCGAAGGCCGCGATATCGCCCa                                 >  1:3345968/1‑69 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATc                                      >  1:2054037/1‑64 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGccc                                  >  1:651791/1‑68 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGGCCa                                 >  1:1276404/1‑69 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCa                                 >  1:3317977/1‑69 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCa                                 >  1:739051/1‑69 (MQ=255)

                               cGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCa                                 >  1:1969774/1‑69 (MQ=255)

                                         aaTTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCa                         <  1:2621019/67‑1 (MQ=255)

                                            ttCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:841597/1‑69 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:3304357/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:3366317/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:853615/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:899582/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:3072356/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:1168676/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2961439/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:153016/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:1691102/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:1869186/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:1924412/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2999752/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2437729/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:255717/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:125751/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:262811/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2655435/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2886394/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:2740800/50‑1 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGtt                   >  1:1486114/1‑69 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:2638315/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:1532375/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:1620916/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:317895/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGt                    >  1:425227/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACAGCGTCAGTCGt                    >  1:1566571/1‑68 (MQ=255)

                                             tctcGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGCCCGCGATATCGCCCACGGCGTCAg                        >  1:2414223/1‑64 (MQ=255)

                                                cgcgTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTac                 <  1:1182991/68‑1 (MQ=255)

                                                cgcgTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTac                 <  1:1203813/68‑1 (MQ=255)

                                                cgcgTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTac                 <  1:1696156/68‑1 (MQ=255)

                                                cgcgTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTac                 <  1:2249762/68‑1 (MQ=255)

                                                     ggtgAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGc           <  1:1772276/69‑1 (MQ=255)

                                                      gtgAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGc           <  1:316824/68‑1 (MQ=255)

                                                          aGTCATCCGTATTTGCAACCTCGAAGGTCGCGATATc                                      <  1:1403062/37‑1 (MQ=38)

                                                          aGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTac                 <  1:153284/58‑1 (MQ=255)

                                                              aTCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGCGATGCATTa  <  1:97083/69‑1 (MQ=255)

                                                              aTCAGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGCGATGCATTa  <  1:1076334/69‑1 (MQ=255)

                                                               tCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGCGATGCATTa  <  1:2169883/68‑1 (MQ=255)

GCAGCTCCCTGTTGCCGAAAGGCATTAAAAGCGTGACTGGCAATTTCTCGCGTGGTGAAGTCATCCGCATTTGCAACCTCGAAGGCCGCGATATCGCCCACGGCGTCAGTCGTTACAACAGCGATGCATTA  >  minE/229864‑229994

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