BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	340,151	T→A	5.5%	A14A (GCT→GCA)	ybaS →	predicted glutaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	340,151	0	T	A	5.5%	365.9 / 7.2	123	A14A (GCT→GCA)	ybaS	predicted glutaminase
Reads supporting (aligned to +/- strand): ref base T (61/55); new base A (5/2); total (66/57)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTAC  >  Exported/340083‑340218

gCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   <  1:1524920/71‑1 (MQ=255)

gCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATa                                                                   >  1:39477/1‑71 (MQ=255)

 cATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTac                                                                  <  1:2251240/71‑1 (MQ=255)

  aTCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTaca                                                                 >  1:3505767/1‑71 (MQ=255)

      ataAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACAccc                                                              >  1:3623992/1‑70 (MQ=255)

        aaaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   >  1:798113/1‑63 (MQ=255)

        aaaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   >  1:3403000/1‑63 (MQ=255)

          aaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAttt                                                         <  1:963568/71‑1 (MQ=255)

            aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtca                                                       >  1:1291949/1‑71 (MQ=255)

            aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtca                                                       >  1:3025644/1‑71 (MQ=255)

            aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtca                                                       >  1:134994/1‑71 (MQ=255)

            aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtca                                                       >  1:3704817/1‑71 (MQ=255)

            aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtca                                                       <  1:1855518/71‑1 (MQ=255)

              gAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtcact                                                     >  1:3309208/1‑71 (MQ=255)

                gTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCaa                                                            >  1:1308473/1‑62 (MQ=255)

                gTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCaa                                                            >  1:2626245/1‑62 (MQ=255)

                gTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtcactc                                                    <  1:1684312/70‑1 (MQ=255)

                    aCAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACtt                                                <  1:983308/70‑1 (MQ=255)

                    aCAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACtt                                                <  1:1533294/70‑1 (MQ=255)

                      aaaaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTcactcact                                                 >  1:871260/1‑67 (MQ=255)

                      aaaaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACtt                                                <  1:3489372/68‑1 (MQ=255)

                      aaaaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACtt                                                <  1:2734039/68‑1 (MQ=255)

                       aaaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAAcg                                            >  1:3074717/1‑71 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCtt                                                                    >  1:350778/1‑46 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   >  1:1766152/1‑47 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   >  1:3539593/1‑47 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTaa                                              >  1:3485213/1‑68 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTaa                                              >  1:1318206/1‑68 (MQ=255)

                        aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAAcgg                                           <  1:2276246/71‑1 (MQ=255)

                         aGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACAccc                                                              <  1:1017337/51‑1 (MQ=255)

                          gATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAAcggcg                                         <  1:2122508/71‑1 (MQ=255)

                             gTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGa                                       >  1:1586805/1‑70 (MQ=255)

                               tAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTacac                                                                >  1:2508758/1‑43 (MQ=255)

                                aGATGCAAACAAATTACAGCAGGCAGTGGTTCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaa                                    <  1:2892118/70‑1 (MQ=255)

                                aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaa                                    >  1:2357046/1‑70 (MQ=255)

                                aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaa                                    >  1:3483878/1‑70 (MQ=255)

                                aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaa                                    >  1:3235864/1‑70 (MQ=255)

                                 gATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaaaa                                  <  1:3650925/71‑1 (MQ=255)

                                 gATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaaa                                   <  1:2520764/70‑1 (MQ=255)

                                  aTGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTa                                                                   >  1:2562324/1‑37 (MQ=255)

                                  aTGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACAccc                                                              >  1:486107/1‑42 (MQ=255)

                                  aTGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAAt                                 >  1:3614264/1‑71 (MQ=255)

                                  aTGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAAt                                 >  1:2914856/1‑71 (MQ=255)

                                    gCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGc                               >  1:1525982/1‑71 (MQ=255)

                                      aaacaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAAcg                                            >  1:2428428/1‑56 (MQ=255)

                                      aaacaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGa                                       <  1:2243922/61‑1 (MQ=255)

                                      aaacaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGa                                       <  1:2877519/61‑1 (MQ=255)

                                       aacaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAAcg                                            <  1:2167292/55‑1 (MQ=255)

                                       aacaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGcc                              <  1:781641/69‑1 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGTATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                           >  1:2560465/1‑71 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtcact                                                     >  1:2442141/1‑45 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtcac                                                      >  1:687092/1‑44 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTtc                                                        <  1:531625/42‑1 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGa                            >  1:2212065/1‑70 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGa                            >  1:3174903/1‑70 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                           >  1:627136/1‑71 (MQ=255)

                                        acaaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                           >  1:1622470/1‑71 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCg                             >  1:206463/1‑68 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCg                             >  1:916579/1‑68 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGa                            >  1:3535681/1‑69 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAtt                          <  1:1066246/71‑1 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAtt                          <  1:882980/71‑1 (MQ=255)

                                         caaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAtt                          <  1:2128453/71‑1 (MQ=255)

                                          aaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTa                         <  1:3715603/71‑1 (MQ=255)

                                          aaaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTa                         <  1:3719762/71‑1 (MQ=255)

                                           aaTTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTAc                        <  1:206063/71‑1 (MQ=255)

                                             ttACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAt                      <  1:1017477/71‑1 (MQ=255)

                                              tACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAtt                     >  1:809462/1‑71 (MQ=255)

                                               aCAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAtt                     >  1:2855514/1‑70 (MQ=255)

                                               aCAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTc                    >  1:1680012/1‑71 (MQ=255)

                                               aCAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTc                    >  1:1604457/1‑71 (MQ=255)

                                                cagcagGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTcc                   >  1:2288265/1‑71 (MQ=255)

                                                cagcagGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTcc                   >  1:2953773/1‑71 (MQ=255)

                                                  gcaggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCt                 <  1:1076105/71‑1 (MQ=255)

                                                  gcaggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCt                 <  1:1583022/71‑1 (MQ=255)

                                                  gcaggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCt                 <  1:1720156/71‑1 (MQ=255)

                                                   caggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCtt                <  1:2730055/71‑1 (MQ=255)

                                                   caggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCtt                <  1:3242789/71‑1 (MQ=255)

                                                    aggcagTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAt                      >  1:1680290/1‑64 (MQ=255)

                                                     ggcagTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTc              <  1:1530296/71‑1 (MQ=255)

                                                         gTGGATCAGGCTTACACCCAATTTCACTCACTTAAcggcg                                         <  1:3607745/40‑1 (MQ=255)

                                                         gTGGATCAGGCTTACACCCAATTTCACTCACTTAAcggcg                                         <  1:782985/40‑1 (MQ=255)

                                                          tGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAtt                     >  1:1647486/1‑59 (MQ=255)

                                                             aTCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTa                         <  1:2596524/52‑1 (MQ=255)

                                                             aTCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTa                         <  1:866861/52‑1 (MQ=255)

                                                             aTCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGaa       <  1:638152/70‑1 (MQ=255)

                                                             aTCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAAt      <  1:2696338/71‑1 (MQ=255)

                                                               cAGGCTTACCCCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATg     >  1:3693797/1‑70 (MQ=255)

                                                               cAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATg     >  1:3434290/1‑70 (MQ=255)

                                                               cAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCttt               >  1:767006/1‑60 (MQ=255)

                                                                aGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGt    <  1:1833134/70‑1 (MQ=255)

                                                                aGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGt    <  1:318154/70‑1 (MQ=255)

                                                                aGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTa   >  1:2285345/1‑71 (MQ=255)

                                                                 ggCTTACACCCAATTTCACTCACTTAACGGCGGACaaaa                                  <  1:2782135/39‑1 (MQ=255)

                                                                 ggCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTAc  <  1:3392863/71‑1 (MQ=255)

                                                                  gCTTACACCCAATTTCACTCACTTAACGGCGGACAAGATGCCGAt                           <  1:581038/45‑1 (MQ=255)

                                                                  gCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                           <  1:3319986/45‑1 (MQ=255)

                                                                  gCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                           <  1:49332/45‑1 (MQ=255)

                                                                  gCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTAc  >  1:721290/1‑70 (MQ=255)

                                                                  gCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTAc  >  1:313678/1‑70 (MQ=255)

GCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGTAC  >  Exported/340083‑340218

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