BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,602,757	T→A	21.9%	L260F (TTA→TTT)	trkH ←	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,602,757	0	T	A	21.9%	130.8 / 35.7	64	L260F (TTA→TTT)	trkH	potassium transporter
Reads supporting (aligned to +/- strand): ref base T (28/22); new base A (8/6); total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

AACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAA  >  W3110S.gb/3602692‑3602824

aaCTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAg                                                                >  1:3090583/1‑71 (MQ=255)

aaCTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAg                                                                >  1:2516926/1‑71 (MQ=255)

aaCTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAg                                                                >  1:1955165/1‑71 (MQ=255)

     tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                            <  1:1390253/70‑1 (MQ=255)

     tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGaa                                                            >  1:1488456/1‑70 (MQ=255)

       cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTa                                                                    >  1:3481297/1‑60 (MQ=255)

       cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAg                                                         <  1:621437/71‑1 (MQ=255)

        gCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGt                                                        <  1:117587/71‑1 (MQ=255)

         ccGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaa                                                           <  1:1513989/67‑1 (MQ=255)

           gATAAACATGCGAAATTCAGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:2200486/1‑69 (MQ=255)

            aTAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:3741828/1‑68 (MQ=255)

            aTAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:2754876/1‑68 (MQ=255)

            aTAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAg                                                         <  1:4305064/66‑1 (MQ=255)

            aTAAAAATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:4515383/1‑68 (MQ=255)

             tAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGa                                                   >  1:2877061/1‑71 (MQ=255)

              aaaCATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGAc                                                  >  1:3974648/1‑71 (MQ=255)

                 cATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGt                                               <  1:4386222/71‑1 (MQ=255)

                  aTGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCg                                                >  1:1984696/1‑69 (MQ=255)

                  aTGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCg                                                >  1:3415955/1‑69 (MQ=255)

                    gCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGt                                            >  1:1913404/1‑71 (MQ=255)

                     cGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:1665389/1‑59 (MQ=255)

                     agAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGt                                            >  1:3698853/2‑70 (MQ=255)

                      gAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTa                                          >  1:2377654/1‑71 (MQ=255)

                         aTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                            <  1:2276192/50‑1 (MQ=255)

                         aTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                            <  1:785670/50‑1 (MQ=255)

                         aTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCa                                                     >  1:2351225/1‑57 (MQ=255)

                         aTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAg                                       >  1:2878118/1‑71 (MQ=255)

                         aTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAg                                       >  1:4619946/1‑71 (MQ=255)

                          ttCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaa                                                           <  1:293759/50‑1 (MQ=255)

                          ttCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                            <  1:2129091/49‑1 (MQ=255)

                          ttCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGc                                      <  1:3570194/71‑1 (MQ=255)

                          ttCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGc                                      <  1:3551252/71‑1 (MQ=255)

                           tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:1990299/1‑53 (MQ=255)

                           tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       >  1:1309334/1‑53 (MQ=255)

                            ccGGATCGCGCCAATAAACCTTAAGACTACGCCAACTAAAc                                                                  >  1:69660/1‑41 (MQ=25)

                             cGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaaa                                                          >  1:3814053/1‑48 (MQ=255)

                             cGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaaa                                                          >  1:2000430/1‑48 (MQ=255)

                              ggATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                     <  1:3576878/36‑1 (MQ=255)

                              ggATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCg                                    >  1:2115701/1‑69 (MQ=255)

                               gATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                 <  1:1890008/71‑1 (MQ=255)

                                aTCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                            <  1:2217228/43‑1 (MQ=255)

                                aTCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGa                                                   >  1:2672301/1‑52 (MQ=255)

                                aTCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTAAAg                                       >  1:3763681/1‑64 (MQ=255)

                                 tCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgaga                                <  1:2840431/70‑1 (MQ=255)

                                  cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                               >  1:2966394/1‑70 (MQ=255)

                                      ccAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGt                                            >  1:4127132/1‑53 (MQ=255)

                                           aaaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                       <  1:1866978/37‑1 (MQ=255)

                                            aaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTAc                                         <  1:1915430/50‑1 (MQ=255)

                                            aaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTAc                                         <  1:1920834/50‑1 (MQ=255)

                                            aaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAgg                        >  1:1119980/1‑67 (MQ=255)

                                            aaCCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAg                     <  1:974222/70‑1 (MQ=255)

                                             aCCTTCAGACTACGCCCACTAAACAGTGAAAAGTg                                                       <  1:4184436/35‑1 (MQ=255)

                                               cTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAg                 >  1:3435503/1‑71 (MQ=255)

                                                  cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCaa              >  1:3472095/1‑71 (MQ=255)

                                                  cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCa               <  1:4153230/70‑1 (MQ=255)

                                                    gACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                 <  1:3081314/50‑1 (MQ=255)

                                                     aCTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGt                                            <  1:4152088/38‑1 (MQ=255)

                                                     aCTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGa           <  1:3939404/71‑1 (MQ=255)

                                                      cTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgaga                                >  1:3126349/1‑49 (MQ=255)

                                                             cACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                 <  1:482569/41‑1 (MQ=255)

                                                             cACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt    >  1:4326160/1‑70 (MQ=255)

                                                             cACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTa   <  1:4441328/71‑1 (MQ=255)

                                                             cACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTa   <  1:2337145/71‑1 (MQ=255)

                                                              aCTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATa         >  1:2885531/1‑64 (MQ=255)

                                                               cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt    <  1:1890182/68‑1 (MQ=255)

                                                                ttAACAGTGAAAACTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt    <  1:4404040/67‑1 (MQ=255)

                                                                 tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTaa  >  1:4363388/1‑68 (MQ=255)

AACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAA  >  W3110S.gb/3602692‑3602824

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