BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,341,007	A→G	14.8%	S143S (TCA→TCG)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,341,007	0	A	G	14.8%	150.0 / 12.8	87	S143S (TCA→TCG)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base A (18/56); new base G (8/5); total (26/61)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGGCGCATGGAACTCCTCC  >  minE/2340942‑2341064

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGatcatc                                                         >  1:2324372/1‑68 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGatcatc                                                         >  1:2489119/1‑68 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGatcatc                                                         >  1:2501549/1‑68 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGatcatc                                                         >  1:1708916/1‑68 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGatcatc                                                         >  1:3001563/1‑68 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:209574/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:2388198/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1920755/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:246245/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:2631288/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1769569/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:2656791/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1544855/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1415193/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:3151730/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1263200/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:639825/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:1219433/1‑69 (MQ=255)

gCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCatct                                                        >  1:727098/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:1412003/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:2835538/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:3078014/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:2177183/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:207628/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:1899755/1‑69 (MQ=255)

      acCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTcc                                                  >  1:1113866/1‑69 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTGGGATCGTCTATCTCCGGt                                               <  1:2329762/69‑1 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGt                                               <  1:3271365/69‑1 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGt                                               <  1:2637868/69‑1 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCgg                                                >  1:2147883/1‑68 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGt                                               >  1:1326353/1‑69 (MQ=255)

         gTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGt                                               <  1:501309/69‑1 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:1109575/1‑69 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:2203226/1‑69 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:2253972/1‑69 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:2280850/1‑69 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:827370/1‑69 (MQ=255)

            cTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTcc                                            >  1:1775233/1‑69 (MQ=255)

                   ggCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGtt                                      <  1:1220509/68‑1 (MQ=255)

                   ggCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGtt                                      <  1:459651/68‑1 (MQ=255)

                   ggCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGtt                                      <  1:1841526/68‑1 (MQ=255)

                    gCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGtt                                      <  1:2100998/67‑1 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTAtctc                                                   >  1:910611/1‑53 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCg                                                 >  1:1216014/1‑55 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTg                                     >  1:1382315/1‑67 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTGc                                    >  1:3348421/1‑68 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTGc                                    >  1:2567359/1‑68 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTGc                                    >  1:1956582/1‑68 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTGCt                                   >  1:399201/1‑69 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCGTCTATCTCCGGTTCCCACGTTGCt                                   >  1:527834/1‑69 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCt                                   <  1:3085775/69‑1 (MQ=255)

                     ccAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCt                                   <  1:1050038/69‑1 (MQ=255)

                     acAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCt                                   <  1:2054554/68‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:3118451/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2939143/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2699117/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2676098/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1283591/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:668334/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2184697/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:260581/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:201256/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2443834/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2472958/69‑1 (MQ=255)

                           ccTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:721530/69‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:1820759/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:350802/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:875131/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:2632917/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:2861800/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:605816/49‑1 (MQ=255)

                                   cttcCTGTTCGATGACATCCGATTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:2568646/49‑1 (MQ=255)

                                    ttcCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCAc                                         <  1:1278710/48‑1 (MQ=255)

                                               tGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGGCGCATGGa          <  1:2207510/68‑1 (MQ=255)

                                               tGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGGCGCATGGa          <  1:3299958/68‑1 (MQ=255)

                                                      cGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGGCGCATGGAACtcctcc  >  1:2819573/1‑69 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:530648/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:789833/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:715158/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:913726/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:648386/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2670837/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2034538/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2339334/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1889417/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2525795/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2232544/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1784282/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2644907/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2177286/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:2681732/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1446871/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:286359/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1360099/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1321094/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:3243725/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:3384581/38‑1 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGgcg                >  1:335026/1‑51 (MQ=255)

                                                          tCGGATCATCTATCTCCGGTTCCCACGTT‑CTATcgac                             <  1:2553835/37‑1 (MQ=38)

                                                             gATCATCTATCTCCGGTTCCCACGTTGCTATcgac                             <  1:1762436/35‑1 (MQ=255)

GCCGACACCGTACTGTTCGGGCCAGAACCTGAATTCTTCCTGTTCGATGACATCCGTTTCGGATCATCTATCTCCGGTTCCCACGTTGCTATCGACGATATCGAAGGCGCATGGAACTCCTCC  >  minE/2340942‑2341064

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