BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,130,079	G→A	42.6%	T61T (ACC→ACT)	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,130,079	0	G	A	42.6%	25.2 / 55.8	78	T61T (ACC→ACT)	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base G (11/33); new base A (18/15); total (30/48)
Fisher's exact test for biased strand distribution p-value = 9.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

ATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATT  >  minE/1130014‑1130141

aTCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGt                                                              <  1:3006189/68‑1 (MQ=255)

aTCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGt                                                              <  1:3114614/68‑1 (MQ=255)

aTCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGt                                                              <  1:1416903/68‑1 (MQ=255)

          cTAACAGGTGCGCACAGGAGTGACGAATGATCTCAAGACCTTCTTCGTCTTTGGCGGTAATGATCGAc                                                    <  1:97268/68‑1 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCTGTAATGATCGACAg                                                  >  1:2630619/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGAc                                                    <  1:533976/67‑1 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACa                                                   >  1:412120/1‑68 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:184766/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:465974/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:1929244/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:481166/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:2038510/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:2082464/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:2253872/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:289536/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:275343/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:2761546/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:1827960/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:982138/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:861770/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:712053/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:686795/1‑69 (MQ=255)

           tAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATAGACa                                                   >  1:3271236/1‑68 (MQ=255)

           tAAAAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCAGTAATGATCGACAg                                                  >  1:1742633/1‑69 (MQ=255)

                                aCGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCa                                  >  1:1901784/1‑64 (MQ=255)

                                aCGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCa                                  >  1:3000191/1‑64 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     <  1:359555/68‑1 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     >  1:530881/1‑68 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     <  1:2924383/68‑1 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     >  1:2676976/1‑68 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     >  1:954498/1‑68 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAg                     >  1:1754411/1‑68 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGc                    >  1:2975730/1‑69 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGc                    >  1:2968996/1‑69 (MQ=255)

                                         ctcCAGACCTTCTTCGTCTTTGG‑GGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGc                    >  1:125557/1‑68 (MQ=255)

                                            cAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 >  1:1485284/1‑69 (MQ=255)

                                            cAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 >  1:2907481/1‑69 (MQ=255)

                                              gACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:811271/67‑1 (MQ=255)

                                              gACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:2570314/67‑1 (MQ=255)

                                              gACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:820920/67‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:3081579/69‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:983363/69‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:2545923/69‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:1208065/69‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:1239397/69‑1 (MQ=255)

                                               aCCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAAc              <  1:1315307/69‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:661307/64‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:2001935/64‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:2050952/64‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:760305/64‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:2137080/64‑1 (MQ=255)

                                                 cttcttCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATc                 <  1:2413480/64‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:1228309/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:810244/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:761051/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:89426/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:744197/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:1143720/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:1451978/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:1622552/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:47646/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:2513972/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:341475/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:2216356/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:2056504/69‑1 (MQ=255)

                                                   tcttcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:2207532/69‑1 (MQ=255)

                                                   tcctcGTCTTTGGCAGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGt          <  1:632698/66‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:1080867/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:2658964/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:1122832/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:2353986/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:310003/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:3363051/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:2055217/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:1332150/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:1993700/69‑1 (MQ=255)

                                                           tttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:1523228/69‑1 (MQ=255)

                                                            ttGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCAtt  <  1:2394040/68‑1 (MQ=255)

ATCGCGTGCCCTAACAGGTGCGCACAGGAGTGACGAATGATCTCCAGACCTTCTTCGTCTTTGGCGGTAATGATCGACAGTTGTGCGTCGTTTTCAATCAGATCGCAAGCATCAACCAGTTCGCCATT  >  minE/1130014‑1130141

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