BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,870,877	A→C	20.0%	I268I (ATA→ATC)	yqcE →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,870,877	0	A	C	20.0%	78.4 / 17.5	50	I268I (ATA→ATC)	yqcE	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (12/28); new base C (2/8); total (14/36)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATCC  >  minE/1870817‑1870937

cGAAATGTACTGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:193491/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:1099031/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:2279593/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:1880078/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:1861809/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:2795324/67‑1 (MQ=255)

cGAAATGTACGGCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:3433590/67‑1 (MQ=255)

 gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                                                     <  1:3127912/69‑1 (MQ=255)

 gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                                                     <  1:1633831/69‑1 (MQ=255)

 gAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCcgcgcg                                                     <  1:120563/69‑1 (MQ=255)

    aTGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCtg                                                  <  1:132556/69‑1 (MQ=255)

        acggCATGTCACTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATCTTCCGCGCGTTATGc                                               >  1:3067380/3‑68 (MQ=255)

        acggCATGTCACTGGTGGCGGCGAGCTACATGAGGATTGTGATCAACAAAATCTTCCGCGCGTTATGc                                               >  1:71114/3‑68 (MQ=255)

           gCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGcc                                           >  1:843240/1‑69 (MQ=255)

           gCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGcc                                           >  1:2215134/1‑69 (MQ=255)

                tCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACtt                                      >  1:321762/1‑69 (MQ=255)

                     ggtggCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTggcgg                                 >  1:2115985/1‑69 (MQ=255)

                        ggcggcGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCAt                              <  1:869069/69‑1 (MQ=255)

                        ggcggcGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCAt                              <  1:380765/69‑1 (MQ=255)

                        ggcggcGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCAt                              <  1:1991811/69‑1 (MQ=255)

                        ggcggcGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCAt                              <  1:1482753/69‑1 (MQ=255)

                           ggcgAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAAt                           <  1:193108/69‑1 (MQ=255)

                           ggcgAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAAt                           <  1:2304623/69‑1 (MQ=255)

                           ggcgAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAAt                           <  1:1297578/69‑1 (MQ=255)

                           ggcgAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAAt                           <  1:1628980/69‑1 (MQ=255)

                               aGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATcacc                       <  1:2653783/69‑1 (MQ=255)

                                gCTACATGGGGATTGTGATCAACAAAATCTTCcgc                                                        <  1:1185319/35‑1 (MQ=255)

                                    cATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCaccacc                    >  1:609982/1‑67 (MQ=255)

                                    cATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCaccacc                    >  1:2617009/1‑67 (MQ=255)

                                      tGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTa                  <  1:2614565/67‑1 (MQ=255)

                                        gggATTGTGATCAACAAAATATTCCGCGCGCTGTGCgg                                             >  1:3206293/1‑38 (MQ=255)

                                        gggATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAg               >  1:2712380/1‑68 (MQ=255)

                                        gggATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAg               >  1:1632640/1‑68 (MQ=255)

                                        gggATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAg               >  1:2877465/1‑68 (MQ=255)

                                        gggATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAg               >  1:3123069/1‑68 (MQ=255)

                                                gATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:2183145/68‑1 (MQ=255)

                                                gATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:3028375/68‑1 (MQ=255)

                                                gATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:2659602/68‑1 (MQ=255)

                                                gATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:2212903/68‑1 (MQ=255)

                                                gATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:2095277/68‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:2481277/69‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:2400284/69‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:384510/69‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:1254531/69‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:834441/69‑1 (MQ=255)

                                                   caacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   >  1:1044598/1‑69 (MQ=255)

                                                    aacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATcc  <  1:1857035/69‑1 (MQ=255)

                                                    aacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATcc  <  1:3329635/69‑1 (MQ=255)

                                                    aacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATcc  <  1:3383254/69‑1 (MQ=255)

                                                    aacaaAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATc   <  1:2871605/68‑1 (MQ=255)

                                                            atTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGa       <  1:1852982/56‑1 (MQ=255)

CGAAATGTATGGCATGTCGCTGGTGGCGGCGAGCTACATGGGGATTGTGATCAACAAAATATTCCGCGCGCTGTGCGGCCCACTTGGCGGCATAATCACCACCTACAGCAAAGTGAAATCC  >  minE/1870817‑1870937

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: `‑`