BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,603,818	A→G	24.1%	Y285H (TAC→CAC)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,603,818	0	A	G	24.1%	110.8 / 34.9	65	Y285H (TAC→CAC)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (32/17); new base G (5/10); total (37/28)
Fisher's exact test for biased strand distribution p-value = 3.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AACCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGTT  >  minE/2603752‑2603885

aaCCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTagag                                                                  >  1:3546896/1‑70 (MQ=255)

  ccGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAAt                                                               >  1:660789/1‑71 (MQ=255)

  ccGAGAATTTCATGATGCACGTTGTTGATGGCAAACAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAAt                                                               >  1:2069704/1‑71 (MQ=255)

    gagaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATCTTGGCATCTGGGCGTAGAGAATGt                                                             <  1:1125236/71‑1 (MQ=255)

     agaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCAACTGGTCGTAGAGAATGTa                                                            <  1:3599896/71‑1 (MQ=255)

      gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTa                                                            >  1:850067/1‑70 (MQ=255)

       aaTTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                           >  1:2570853/1‑70 (MQ=255)

         tttCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTa                                                            <  1:2435869/67‑1 (MQ=255)

           tCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaaa                                                         >  1:2934064/1‑68 (MQ=255)

           tCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaaaa                                                         >  1:769233/1‑68 (MQ=255)

           tCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaaaa                                                         >  1:491545/1‑68 (MQ=255)

            cATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                      <  1:2872141/70‑1 (MQ=255)

              tgatgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAAcaccac                                                   <  1:674069/71‑1 (MQ=255)

                 tgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGc                                                >  1:3093169/1‑71 (MQ=255)

                 tgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGc                                                <  1:2608946/71‑1 (MQ=255)

                  gCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCt                                               >  1:1056477/1‑71 (MQ=255)

                    aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaa                                                           <  1:1351936/57‑1 (MQ=255)

                    aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaa                                                           <  1:1786601/57‑1 (MQ=255)

                    aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaa                                                           <  1:231278/57‑1 (MQ=255)

                    aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCtt                                              >  1:1880749/1‑70 (MQ=255)

                     cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTc                                             >  1:3100600/1‑70 (MQ=255)

                     cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTc                                             >  1:1453095/1‑70 (MQ=255)

                      gttgttGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaaa                                                         >  1:1282695/1‑57 (MQ=255)

                      gttgttGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGt                                                             >  1:3399089/1‑53 (MQ=255)

                        tgttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaaaa                                                         >  1:2942599/1‑55 (MQ=255)

                        tgttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaaaa                                                         >  1:230738/1‑55 (MQ=255)

                         gttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGt                                                             >  1:3507454/1‑50 (MQ=255)

                          ttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCAt                                       <  1:1706617/71‑1 (MQ=255)

                            gATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                      >  1:2145689/1‑54 (MQ=255)

                                 cAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACa                                <  1:2828364/71‑1 (MQ=255)

                                  aaaGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTa                                                            >  1:3184626/1‑42 (MQ=255)

                                  aaaGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAAcacc                                                     >  1:468926/1‑49 (MQ=255)

                                   aagaagTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaa                                                           <  1:2070797/42‑1 (MQ=255)

                                   aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                              <  1:1663592/71‑1 (MQ=255)

                                   aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACAACCGCTTCGa                                           >  1:3816979/1‑58 (MQ=255)

                                    agaagTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAg                             >  1:798383/1‑71 (MQ=255)

                                    agaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAg                             >  1:1688710/1‑71 (MQ=255)

                                    agaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAg                             >  1:1489928/1‑71 (MQ=255)

                                    agaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCAACGCtt                                              >  1:1105653/1‑54 (MQ=255)

                                      aagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAGAACACCACCGCTTCGAGCATCAGGACAAAGGc                           >  1:3755215/1‑71 (MQ=255)

                                      aagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGa                                  >  1:2485400/1‑64 (MQ=255)

                                       agTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCg                          <  1:2793061/71‑1 (MQ=255)

                                        gTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGa                         <  1:1330669/71‑1 (MQ=255)

                                         ttCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCCAc                        <  1:2481611/71‑1 (MQ=255)

                                         ttCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCCAc                        <  1:2477073/71‑1 (MQ=255)

                                          tCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACa                                <  1:131630/62‑1 (MQ=255)

                                          tCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAgg                            <  1:3364640/66‑1 (MQ=255)

                                           cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACa                       >  1:894439/1‑70 (MQ=255)

                                            aGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACaaa                     >  1:3446734/1‑71 (MQ=255)

                                            aGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACaaa                     >  1:1400368/1‑71 (MQ=255)

                                             gCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACaaa                     >  1:3084756/1‑70 (MQ=255)

                                               gATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAg                                          >  1:2956320/1‑47 (MQ=255)

                                                aTGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtt                 <  1:923544/71‑1 (MQ=255)

                                                  gTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAttt                >  1:2894875/1‑70 (MQ=255)

                                                   ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa              >  1:2431651/1‑71 (MQ=255)

                                                   ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa              <  1:1923088/71‑1 (MQ=255)

                                                   ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa              <  1:1712887/71‑1 (MQ=255)

                                                     ggCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCa                                     <  1:2640573/46‑1 (MQ=255)

                                                     ggCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtt                 >  1:2405945/1‑66 (MQ=255)

                                                     ggCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtt                 >  1:900343/1‑66 (MQ=255)

                                                           tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCg                          <  1:3175830/51‑1 (MQ=255)

                                                           tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCg                          <  1:1209576/51‑1 (MQ=255)

                                                           tGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGccc      >  1:3829035/1‑71 (MQ=255)

                                                             ggCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGcc       <  1:3178159/68‑1 (MQ=255)

                                                             ggCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCgg    <  1:2444561/71‑1 (MQ=255)

                                                               cGTTGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGtt  <  1:3006500/71‑1 (MQ=255)

                                                                  agagAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGccc      >  1:1158590/1‑64 (MQ=255)

AACCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGTT  >  minE/2603752‑2603885

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