BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,432,742	T→C	5.5%	K37K (AAA→AAG)	trxA ←	thioredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,432,742	0	T	C	5.5%	280.8 / 6.0	90	K37K (AAA→AAG)	trxA	thioredoxin 1
Reads supporting (aligned to +/- strand): ref base T (55/30); new base C (3/2); total (58/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 108 total aligned reads displayed.

TGTTCAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGT  >  Exported/2432672‑2432810

tGTTCAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtt                                                                      >  1:1333984/1‑71 (MQ=255)

tGTTCAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtt                                                                      >  1:6409295/1‑71 (MQ=255)

  ttCAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtttt                                                                    <  1:4164896/71‑1 (MQ=255)

    cAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtttt                                                                    >  1:2881425/1‑69 (MQ=255)

    cAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTg                                                                   <  1:1058218/70‑1 (MQ=255)

     aGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTg                                                                   <  1:3234939/69‑1 (MQ=255)

     aGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTg                                                                   <  1:1823576/69‑1 (MQ=255)

        tttGCAACGGTCAGTTTGCGCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACg                                                               >  1:3810735/1‑70 (MQ=255)

        tttGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACg                                                               >  1:2895492/1‑70 (MQ=255)

        tttGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACg                                                               >  1:523211/1‑70 (MQ=255)

            cAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGAcc                                                           <  1:4971117/70‑1 (MQ=255)

              aCGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGca                                                        >  1:4693778/1‑71 (MQ=255)

              aCGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGca                                                        >  1:3594657/1‑71 (MQ=255)

              aCGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGca                                                        >  1:384752/1‑71 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:2247126/1‑70 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:2517074/1‑70 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:1778734/1‑70 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:1761980/1‑70 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:4512591/1‑70 (MQ=255)

                ggTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcac                                                       >  1:153735/1‑70 (MQ=255)

                 gTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtttt                                                                    <  1:4519503/56‑1 (MQ=255)

                  tCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGcaccac                                                    <  1:1471248/71‑1 (MQ=255)

                    aGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCActc                                                  >  1:6797327/1‑71 (MQ=255)

                    aGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCAct                                                   <  1:5376043/70‑1 (MQ=255)

                        tGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGcc                                              >  1:938998/1‑71 (MQ=255)

                             tGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGa                                          >  1:5139522/1‑70 (MQ=255)

                             tGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGa                                          >  1:1662793/1‑70 (MQ=255)

                             tGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACAACTCTGCCCAGa                                          >  1:4557844/1‑70 (MQ=255)

                              gATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCAtt                                                                      <  1:5164318/41‑1 (MQ=255)

                               atatTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGaaa                                        <  1:5482392/70‑1 (MQ=255)

                               atatTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAAt                                       >  1:6322139/1‑71 (MQ=255)

                                tatTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATc                                      <  1:1605802/71‑1 (MQ=255)

                                tatTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATc                                      <  1:1923800/71‑1 (MQ=255)

                                 atTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCTCACCACTCTGCCCAGaa                                         <  1:2435066/67‑1 (MQ=255)

                                  ttCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcga                                    >  1:6600835/1‑71 (MQ=255)

                                   tCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcgac                                   >  1:5278249/1‑71 (MQ=255)

                                   tCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcgac                                   >  1:3516058/1‑71 (MQ=255)

                                    cGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCa                                                                 >  1:4868466/1‑40 (MQ=255)

                                    cGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCa                                                                 >  1:442608/1‑40 (MQ=255)

                                    cGTCAGCGATTTCATCCAGAATCGCGGCGATCATTTTGCa                                                                 >  1:2991461/1‑40 (MQ=255)

                                     gTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcgacg                                  <  1:5412400/70‑1 (MQ=255)

                                        aGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               <  1:1948317/70‑1 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCg                                                          >  1:2152335/1‑40 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCg                                                          >  1:6399522/1‑40 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               >  1:1820562/1‑67 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               >  1:35852/1‑67 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               >  1:3402237/1‑67 (MQ=255)

                                           gATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCg                           <  1:6070039/71‑1 (MQ=255)

                                            aTTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGc                          <  1:3291053/71‑1 (MQ=255)

                                            aTTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGc                          <  1:2532160/71‑1 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGc                                                         >  1:1990091/1‑39 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGc                                                         >  1:6131623/1‑39 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         >  1:1699938/1‑71 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         >  1:5424877/1‑71 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         >  1:3597575/1‑71 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         >  1:3037692/1‑71 (MQ=255)

                                             tttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         >  1:3125715/1‑71 (MQ=255)

                                              ttcttcCAGGATCGGGGCAATCATCTTGCAAGGACCGcacc                                                      <  1:6233763/41‑1 (MQ=19)

                                              ttcttcCAGGATCGGGGCAATCATCTTGCAAGGACCGcacc                                                      <  1:3396147/41‑1 (MQ=19)

                                              ttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcga                                    >  1:6500999/1‑59 (MQ=255)

                                              ttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGccc                        >  1:6837776/1‑71 (MQ=255)

                                              ttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGccc                        >  1:161159/1‑71 (MQ=255)

                                              ttCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGccc                        >  1:410147/1‑71 (MQ=255)

                                                cttCCAGGATCGGGGCAATCATCTTGCAAGGACCGcaccac                                                    >  1:6848461/3‑41 (MQ=19)

                                                cttCCAGGATCGGGGCAATCATCTTGCAAGGACCGcaccac                                                    >  1:1499214/3‑41 (MQ=19)

                                                cATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                       >  1:3450008/1‑70 (MQ=255)

                                                cATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                       <  1:6783380/70‑1 (MQ=255)

                                                cATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                       <  1:1672537/70‑1 (MQ=255)

                                                 aTCCAGAATCGTGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcga                                    <  1:1533356/56‑1 (MQ=255)

                                                 aTCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATcga                                    <  1:3649687/56‑1 (MQ=255)

                                                     aGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAg                                >  1:2097998/1‑56 (MQ=255)

                                                        aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt              <  1:4746724/71‑1 (MQ=255)

                                                        aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt              <  1:3176627/71‑1 (MQ=255)

                                                        aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt              <  1:5267096/71‑1 (MQ=255)

                                                        aTCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCttt              <  1:6723745/71‑1 (MQ=255)

                                                         tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGaa                                         >  1:3561665/1‑43 (MQ=255)

                                                         tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAg                                >  1:4895707/1‑52 (MQ=255)

                                                         tCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAg                                >  1:5120585/1‑52 (MQ=255)

                                                         tCGGGGCGATCATTATGCACGGACCGCACCACTCTGCCCAGAAATcgacg                                  >  1:1633515/1‑50 (MQ=255)

                                                               cGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:1558568/1‑66 (MQ=255)

                                                                gATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACAt      >  1:5142259/1‑71 (MQ=255)

                                                                 aTCATTTTGCACGGACCGCACCACTCTGCCCAGAATTCGACGAGGATCGCCCCGTCCGCTTTGAGTACATc     <  1:759273/71‑1 (MQ=255)

                                                                  tCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATcc    <  1:48951/71‑1 (MQ=255)

                                                                   cATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCg   >  1:649373/1‑71 (MQ=255)

                                                                    aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt  >  1:1503031/1‑71 (MQ=255)

                                                                    aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt  <  1:3303300/71‑1 (MQ=255)

                                                                    aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt  <  1:4080890/71‑1 (MQ=255)

                                                                    aTTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCgt  >  1:373663/1‑71 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               <  1:4697657/41‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               <  1:6503732/41‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAgg                               <  1:3281149/41‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                       <  1:4992384/49‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcccc                       <  1:1585809/49‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         <  1:6570328/47‑1 (MQ=255)

                                                                     ttttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGcc                         <  1:1864014/47‑1 (MQ=255)

                                                                      tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:5547456/1‑59 (MQ=255)

                                                                      tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:3349765/1‑59 (MQ=255)

                                                                      tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:3094848/1‑59 (MQ=255)

                                                                      tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:5638739/1‑59 (MQ=255)

                                                                      tttGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGa            >  1:6249832/1‑59 (MQ=255)

TGTTCAGTTTTGCAACGGTCAGTTTGCCCTGATATTCGTCAGCGATTTCATCCAGAATCGGGGCGATCATTTTGCACGGACCGCACCACTCTGCCCAGAAATCGACGAGGATCGCCCCGTCCGCTTTGAGTACATCCGT  >  Exported/2432672‑2432810

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