BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,758,182	C→T	45.5%	V52I (GTT→ATT)	def ←	peptide deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,758,182	0	C	T	45.5%	34.1 / 136.3	101	V52I (GTT→ATT)	def	peptide deformylase
Reads supporting (aligned to +/- strand): ref base C (31/24); new base T (28/18); total (59/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 103 total aligned reads displayed.

ATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATACGCT  >  minE/2758115‑2758250

attAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTg                                                                   <  1:1776041/71‑1 (MQ=255)

  tAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTggg                                                                 >  1:1546686/1‑71 (MQ=255)

    acacTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGtt                                                               >  1:1138177/1‑71 (MQ=255)

      acTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTgc                                                             <  1:859072/71‑1 (MQ=255)

         aGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTgccgc                                                          >  1:3075603/1‑71 (MQ=255)

         aGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTgccgc                                                          >  1:905834/1‑71 (MQ=255)

          gCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTgccgc                                                          >  1:828709/1‑70 (MQ=255)

          gCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTgccgc                                                          >  1:725301/1‑70 (MQ=255)

             gTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAgg                                                      <  1:3240539/71‑1 (MQ=255)

              ttCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGc                                                     >  1:3834764/1‑71 (MQ=255)

              ttCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGc                                                     >  1:1396536/1‑71 (MQ=255)

              ttCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGc                                                     >  1:2602549/1‑71 (MQ=255)

                cGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGcc                                                    <  1:2626312/70‑1 (MQ=255)

                cGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGcc                                                    <  1:3195254/70‑1 (MQ=255)

                cGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGc                                                     >  1:1975790/1‑69 (MQ=255)

                  tCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTgccgcc                                                         <  1:1995952/63‑1 (MQ=255)

                  tCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTgg                                                                  <  1:2473014/54‑1 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTgccgc                                                          >  1:1301534/1‑60 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTg                                                              >  1:1247222/1‑56 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATa                                                <  1:3652259/70‑1 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATa                                                >  1:2278147/1‑70 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATAc                                               >  1:2763939/1‑71 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATa                                                >  1:3178308/1‑70 (MQ=255)

                    aCGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATAc                                               >  1:3561108/1‑71 (MQ=255)

                      ggTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCTGCCAGGCCAATACct                                             >  1:43947/1‑71 (MQ=255)

                       gTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACctt                                            >  1:2922375/1‑71 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAAc                                                                      >  1:584000/1‑44 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAAc                                                                      >  1:4087794/1‑44 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTgccgc                                                          <  1:1709615/56‑1 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACcttc                                           >  1:2599641/1‑71 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACcttc                                           >  1:915423/1‑71 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACcttc                                           <  1:3599417/71‑1 (MQ=255)

                        ttttCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACcttc                                           <  1:1403808/71‑1 (MQ=255)

                         tttCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACcttct                                          >  1:276385/1‑71 (MQ=255)

                           tCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCttcttc                                        >  1:1635912/1‑71 (MQ=255)

                                aaCATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGt                                   <  1:7585/71‑1 (MQ=255)

                                aaCATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAAGCCAATACCTTCTTCTGCGt                                   <  1:4082126/71‑1 (MQ=255)

                                 aCATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGt                                   >  1:839274/1‑70 (MQ=255)

                                 aCATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTa                                  <  1:3259674/71‑1 (MQ=255)

                                 aCATCAATAACAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTa                                  <  1:3506761/71‑1 (MQ=255)

                                  cATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTa                                  >  1:2719985/1‑70 (MQ=255)

                                  cATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTa                                  >  1:2009934/1‑70 (MQ=255)

                                   aTCAATAACAATGATACGTTGATGGATATCAATCTGGGtt                                                               <  1:1481601/40‑1 (MQ=255)

                                   aTCAATAACAATGATACGTTGATGGATATCAACCTggg                                                                 >  1:4076803/1‑38 (MQ=255)

                                   aTCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCaa                                                  >  1:2237006/1‑53 (MQ=255)

                                     cAATAACAATGATACGTTGATGGATATCAATCTGGGTTgc                                                             >  1:3311712/1‑40 (MQ=255)

                                      aaTAACAATGATACGTTGATGGATATCAATCTGGGTTgc                                                             >  1:946057/1‑39 (MQ=255)

                                      aaTAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCa                                                        >  1:2877082/1‑44 (MQ=255)

                                      aaTAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCaa                                                  >  1:1986142/1‑50 (MQ=255)

                                          aCAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGtctc                         >  1:3745120/1‑71 (MQ=255)

                                          aCAATGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGtctc                         >  1:1240841/1‑71 (MQ=255)

                                          aCAATGATACGTTGATGGATATCAATATGGGTTGCCGCCAGGCCaa                                                  >  1:3548136/1‑46 (MQ=37)

                                           cAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCg                        <  1:1648910/71‑1 (MQ=255)

                                            aaTGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGa                       <  1:2827397/71‑1 (MQ=255)

                                             aTGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGaa                      <  1:1456540/71‑1 (MQ=255)

                                             aTGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGa                       <  1:3726127/70‑1 (MQ=255)

                                             aTGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTAAATCGTCTCGaa                      <  1:3257299/71‑1 (MQ=255)

                                              tGATACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGc                                     <  1:1426503/55‑1 (MQ=255)

                                              cgATCCGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGaa                      <  1:2236917/69‑1 (MQ=255)

                                               gATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGc                                     >  1:494118/1‑54 (MQ=255)

                                                aTACGTTGATGGATATCAATCTGGGTTGCCGCCAgg                                                      >  1:3361532/1‑36 (MQ=255)

                                                aTACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGtc                           <  1:918035/63‑1 (MQ=255)

                                                 tACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACata                  >  1:3941134/1‑71 (MQ=255)

                                                 tACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACata                  >  1:3979517/1‑71 (MQ=255)

                                                 tACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACata                  >  1:3367003/1‑71 (MQ=255)

                                                 tACGTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACat                   <  1:939955/70‑1 (MQ=255)

                                                  aCGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACata                  >  1:2015528/1‑70 (MQ=255)

                                                    gTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAAt                                                 <  1:2817173/37‑1 (MQ=255)

                                                    gTTGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatca               <  1:829533/71‑1 (MQ=255)

                                                    gTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatca               <  1:3304238/71‑1 (MQ=255)

                                                     ttGATGGATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGaa                      <  1:1321770/63‑1 (MQ=255)

                                                       gATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatcatc             <  1:4009198/70‑1 (MQ=255)

                                                       gATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatcatc             <  1:791961/70‑1 (MQ=255)

                                                       gATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatcatc             <  1:2365657/70‑1 (MQ=255)

                                                       gATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATatcatc             <  1:2283975/70‑1 (MQ=255)

                                                         tGGATATCAATCTGGGTTGCCGCCAGGCCAATACcttct                                          <  1:2698730/39‑1 (MQ=255)

                                                         tGGATATCAATCTGGGTTGCCGCCAGGCCAATACcttct                                          <  1:2847917/39‑1 (MQ=255)

                                                          ggATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACa                    <  1:529291/60‑1 (MQ=255)

                                                           gATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGaa                      <  1:714427/57‑1 (MQ=255)

                                                           gATATCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATcgacga        <  1:2625884/71‑1 (MQ=255)

                                                           gATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATcgacga        >  1:2076122/1‑71 (MQ=255)

                                                            atatCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACa                    >  1:646231/1‑58 (MQ=255)

                                                            atatCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATc                              >  1:3792963/1‑48 (MQ=255)

                                                            atatCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATc                              >  1:1239389/1‑48 (MQ=255)

                                                            atatCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATcgacga        <  1:482456/70‑1 (MQ=255)

                                                             tatCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATa      >  1:452646/1‑71 (MQ=255)

                                                             tatCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATa      >  1:1780483/1‑71 (MQ=255)

                                                             tatCAATCTGGGTTGCCGCCAGGCCAATACCTTCTACTGCGTAc                                 >  1:1707308/1‑44 (MQ=38)

                                                             tatCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATa      >  1:3377856/1‑71 (MQ=255)

                                                              atCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATcg            <  1:2164598/64‑1 (MQ=255)

                                                               tCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGtctc                         >  1:3095467/1‑50 (MQ=255)

                                                               tCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATcgac          >  1:3594866/1‑65 (MQ=255)

                                                               tCAATCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATACg    >  1:1064430/1‑71 (MQ=255)

                                                               tCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACata                  >  1:625954/1‑57 (MQ=255)

                                                               tCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACa                    >  1:1671397/1‑55 (MQ=255)

                                                                cAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTa                                  >  1:3997882/1‑40 (MQ=255)

                                                                 aaCCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACa                    >  1:1466568/1‑53 (MQ=255)

                                                                 aaCCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATACGc   >  1:3370727/1‑70 (MQ=255)

                                                                 aaCCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATACGCt  >  1:427761/1‑71 (MQ=255)

                                                                  aCCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCg                                    <  1:670919/36‑1 (MQ=255)

ATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTTGATGGATATCAACCTGGGTTGCCGCCAGGCCAATACCTTCTTCTGCGTACATCGTCTCGAACATATCATCGACGATACGCT  >  minE/2758115‑2758250

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