BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,431,950	A→G	76.9%	V52A (GTT→GCT)	yjfZ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,431,950	0	A	G	76.9%	103.5 / 27.8	52	V52A (GTT→GCT)	yjfZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (5/7); new base G (15/25); total (20/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

AGCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCGCGTGCG  >  W3110S.gb/4431883‑4432012

aGCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             <  1:357964/71‑1 (MQ=255)

aGCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             <  1:2359543/71‑1 (MQ=255)

aGCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             >  1:1670515/1‑71 (MQ=255)

 gCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGt                                                            <  1:114301/71‑1 (MQ=255)

 gCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGt                                                            <  1:829767/71‑1 (MQ=255)

   aTTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTAc                                                          <  1:2229754/71‑1 (MQ=255)

    ttACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                                         <  1:1483917/71‑1 (MQ=255)

     tACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                                         <  1:1644111/70‑1 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             >  1:118145/1‑65 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             >  1:991867/1‑65 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTcagcag                                                             >  1:1476687/1‑65 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCa                                                              >  1:127349/1‑64 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCa                                                              >  1:1047676/1‑64 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTc                                                       >  1:1600908/1‑71 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTc                                                       >  1:2125744/1‑71 (MQ=255)

      acacCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTc                                                       >  1:998308/1‑71 (MQ=255)

              tcgCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                                         <  1:2189933/61‑1 (MQ=255)

                gCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGc                                              <  1:1656478/70‑1 (MQ=255)

                gCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGc                                              <  1:753742/70‑1 (MQ=255)

                    aCTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAAt                                           <  1:2536522/69‑1 (MQ=255)

                     cTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                                         >  1:1508152/1‑54 (MQ=255)

                     cTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCAt                                                 >  1:1544555/1‑62 (MQ=255)

                     cTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCAt                                                 >  1:596072/1‑62 (MQ=255)

                     cTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCAt                                                 >  1:2585342/1‑62 (MQ=255)

                     cTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCAt                                                 >  1:1772942/1‑62 (MQ=255)

                       tATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAAtt                                          <  1:1116/67‑1 (MQ=255)

                       tATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGa                                      >  1:1514994/1‑71 (MQ=255)

                       tATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGa                                      >  1:31712/1‑71 (MQ=255)

                       tATTGAATGCTTGCCAGTCTGAGAGATCAAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAAtt                                          <  1:139264/67‑1 (MQ=255)

                        aTTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATg                                       <  1:2350364/69‑1 (MQ=255)

                         ttGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGt                                                            <  1:567479/47‑1 (MQ=255)

                          tGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTa                                                           >  1:464271/1‑47 (MQ=255)

                            aaTGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATctct                                 <  1:205245/71‑1 (MQ=255)

                              tGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACg                                                         <  1:2126794/45‑1 (MQ=255)

                                cTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAAt                                           >  1:47272/1‑57 (MQ=255)

                                    ccAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGAt                         <  1:1188328/71‑1 (MQ=255)

                                        tCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGAtct                                   <  1:1666432/57‑1 (MQ=255)

                                        tCTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa                     <  1:2086760/71‑1 (MQ=255)

                                         cTGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGAt                                     <  1:1592644/54‑1 (MQ=255)

                                          tGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTa                                         >  1:2783168/1‑49 (MQ=255)

                                          tGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGAt                                     >  1:2381571/1‑53 (MQ=255)

                                          tGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa                     <  1:2017536/69‑1 (MQ=255)

                                          tGAGAGATCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGa                     >  1:2425610/1‑69 (MQ=255)

                                             gagaTCTAGTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTa                <  1:875933/71‑1 (MQ=255)

                                              agaTCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTaa               <  1:38244/71‑1 (MQ=255)

                                              agaTCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTaa               <  1:2545321/71‑1 (MQ=255)

                                                     gTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATg             <  1:753997/66‑1 (MQ=255)

                                                     gTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATg             <  1:2394540/66‑1 (MQ=255)

                                                     gTTTGTGGATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAgc         >  1:1985231/1‑70 (MQ=255)

                                                     gTTTGTGGATCTCAGAAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAgc         >  1:1321871/1‑70 (MQ=255)

                                                            gATCTCAGCAGTACGTCTGTCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCgcgtgcg  <  1:1373684/70‑1 (MQ=255)

                                                            gATCTCAGCAGTACGTCTGCCATCGCAATTATGAt                                     >  1:500236/1‑35 (MQ=255)

                                                            gATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCgcgtgcg  <  1:1383388/70‑1 (MQ=255)

                                                            gATCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCgagtgcg  <  1:2255001/70‑1 (MQ=255)

                                                            gATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCgcgtgcg  <  1:801411/70‑1 (MQ=255)

                                                             aTCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATg             <  1:291141/58‑1 (MQ=255)

                                                             aTCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATg             <  1:1661762/58‑1 (MQ=255)

                                                             aTCTCAGCAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCgcgtgcg  <  1:389086/69‑1 (MQ=255)

AGCATTACACCTCGTCGCTAACTTATTGAATGCTTGCCAGTCTGAGAGATCTAGTTTGTGGATCTCAACAGTACGTCTGCCATCGCAATTATGATCTCTTCCACGATTTGAAAGTAATGTAGCGCGTGCG  >  W3110S.gb/4431883‑4432012

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