BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,787,280	T→C	18.1%	E384E (GAA→GAG)	lysC ←	aspartokinase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,787,280	0	T	C	18.1%	214.1 / 32.4	137	E384E (GAA→GAG)	lysC	aspartokinase III
Reads supporting (aligned to +/- strand): ref base T (72/40); new base C (10/15); total (82/55)
Fisher's exact test for biased strand distribution p-value = 4.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 155 total aligned reads displayed.

AATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGC  >  minE/2787212‑2787347

aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACct                                                                     <  1:1830102/69‑1 (MQ=255)

aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACct                                                                     <  1:1082357/69‑1 (MQ=255)

aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACct                                                                     <  1:1136035/69‑1 (MQ=255)

aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACct                                                                     <  1:2468865/69‑1 (MQ=255)

    ccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACcttct                                                                  >  1:81238/1‑68 (MQ=255)

    ccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACcttct                                                                  >  1:1884289/1‑68 (MQ=255)

    ccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACcttct                                                                  >  1:2268888/1‑68 (MQ=255)

    ccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACcttct                                                                  >  1:235641/1‑68 (MQ=255)

    ccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACcttct                                                                  >  1:2812301/1‑68 (MQ=255)

      tctTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCttcttc                                                                <  1:2547890/68‑1 (MQ=255)

      tctTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCttcttc                                                                <  1:145721/68‑1 (MQ=255)

          tGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccac                                                             <  1:2408839/67‑1 (MQ=255)

          tGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccac                                                             <  1:1964918/67‑1 (MQ=255)

               aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccacctccac                                                       >  1:2525843/1‑68 (MQ=255)

               aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccacctccac                                                       >  1:1428143/1‑68 (MQ=255)

               aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccacctccac                                                       >  1:549368/1‑68 (MQ=255)

               aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTtccacctccac                                                       >  1:591884/1‑68 (MQ=255)

               aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCC‑‑ACCTTCTTCCACCTCCAccc                                                     >  1:1597250/1‑68 (MQ=255)

                  ccgcAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGa                                                   <  1:2513418/69‑1 (MQ=255)

                  ccgcAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGa                                                   <  1:889062/69‑1 (MQ=255)

                   cgcAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGa                                                   <  1:2060710/68‑1 (MQ=255)

                          ttttGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTg                                            >  1:515505/1‑68 (MQ=255)

                          ttttGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTg                                            >  1:2041178/1‑68 (MQ=255)

                          ttttGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTg                                            >  1:1475482/1‑68 (MQ=255)

                          ttttGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTg                                            >  1:1517898/1‑68 (MQ=255)

                                 aGGTCATTGCGAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGc                                    >  1:2540323/1‑69 (MQ=255)

                                 aGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGc                                    >  1:2985699/1‑69 (MQ=255)

                                 aGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGc                                    >  1:1569511/1‑69 (MQ=255)

                                 aGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGc                                    >  1:411332/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGCCCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:775679/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCt                                   >  1:3003775/1‑68 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:2572898/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:185423/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:299436/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:2262701/1‑69 (MQ=255)

                                   gTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                  >  1:2120262/1‑69 (MQ=255)

                                    tCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:684250/1‑69 (MQ=255)

                                    ctaTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:2056072/67‑1 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:3390317/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:500947/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:3237710/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:2343884/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:160749/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:1198412/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGACCAGCGCCAGACCCTCTTCCACTTCCACCCGGCACAGTGCGGAAAGCTcc                                 >  1:2567838/1‑68 (MQ=255)

                                     cATTGCCAATCAACGCGAACAGCGCCAGACCCTCTTCCACTTCCACCCTGCACAGTGCGGAAAGc                                    >  1:2886429/1‑65 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:1786324/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:2026230/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:2940983/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:108077/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:421504/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:3420537/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:1471597/1‑67 (MQ=255)

                                      aTTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 >  1:581626/1‑67 (MQ=255)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:3033242/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:3040999/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:3045232/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:3051691/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:383455/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:2375105/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:1227574/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:2019265/38‑2 (MQ=38)

                                        tGCCAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:2142453/38‑2 (MQ=38)

                                        tGCAAATCAACGCGACCAGCGCCAGACCCTCTtccact                                                            <  1:444112/38‑2 (MQ=37)

                                         gCCAATCAACGCGACCAGCGCCAGTCCCTCTtccact                                                            <  1:3331994/37‑2 (MQ=37)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTtccacctcc                                                         >  1:3295520/1‑40 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCGGACACAGTGCGGAAAGCTCCATc                              >  1:3396095/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:1741599/64‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:3119356/64‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:2816873/64‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:2759282/64‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                 <  1:1419991/64‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcag                            <  1:2344062/69‑1 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:1723530/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:1026693/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:1466714/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:356235/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:3273916/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:1831596/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:1272413/1‑67 (MQ=255)

                                         gCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                              >  1:2445442/1‑67 (MQ=255)

                                          ccAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcag                            <  1:669256/68‑1 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCagag                       >  1:2047748/1‑67 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAgaga                      >  1:1356105/1‑68 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGAt                     <  1:3268743/69‑1 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGAt                     >  1:3247061/1‑69 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGAt                     >  1:3093616/1‑69 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGAt                     >  1:2109768/1‑69 (MQ=255)

                                                aaCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGAt                     <  1:1653860/69‑1 (MQ=255)

                                                           cgcCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCaa          <  1:2919267/69‑1 (MQ=255)

                                                           cgcCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCaa          <  1:313930/69‑1 (MQ=255)

                                                           cgcCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCaa          <  1:1272668/69‑1 (MQ=255)

                                                            gcCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCaa          <  1:1015025/68‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGTTTGCGTCAGCAACGt       <  1:1158849/69‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt       <  1:638635/69‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt       <  1:298792/69‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt       <  1:2691278/69‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt       <  1:697706/69‑1 (MQ=255)

                                                              cAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt       <  1:1803064/69‑1 (MQ=255)

                                                                   cttcttCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGc  <  1:3195094/69‑1 (MQ=255)

AATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGC  >  minE/2787212‑2787347

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