BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,759	A→G	16.9%	I222V (ATT→GTT) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,759	0	A	G	16.9%	250.9 / 38.9	112	I222V (ATT→GTT) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (39/54); new base G (12/7); total (51/61)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCC  >  minE/2570692‑2570829

gCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                     <  1:2452704/71‑1 (MQ=255)

gCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                     <  1:3401954/71‑1 (MQ=255)

 cAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcg                                                                    <  1:2453502/71‑1 (MQ=255)

    tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                 >  1:5307078/1‑71 (MQ=255)

    tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                 >  1:508860/1‑71 (MQ=255)

    tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                 >  1:19330/1‑71 (MQ=255)

      gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                               >  1:5277902/1‑71 (MQ=255)

          aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgcg                                                                <  1:3348057/66‑1 (MQ=255)

          aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                                   <  1:1117837/63‑1 (MQ=255)

          aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAg                                                           >  1:4646603/1‑71 (MQ=255)

           ttGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:3267752/71‑1 (MQ=255)

            tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcg                                                                    >  1:2182206/1‑60 (MQ=255)

            tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          >  1:3519118/1‑70 (MQ=255)

              tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGATCAGACTGTCATTCGCGCGCCGAGCaa                                                        <  1:4146578/70‑1 (MQ=255)

              tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                        <  1:2134918/70‑1 (MQ=255)

                aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:3576077/66‑1 (MQ=255)

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:4609709/71‑1 (MQ=255)

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:4654431/71‑1 (MQ=255)

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:2380771/71‑1 (MQ=255)

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:4739520/71‑1 (MQ=255)

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:2443825/71‑1 (MQ=255)

                   ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                       >  1:750550/1‑66 (MQ=255)

                   ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                  >  1:2094604/1‑71 (MQ=255)

                     aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                        >  1:3434175/1‑47 (MQ=255)

                     aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                >  1:3452316/1‑71 (MQ=255)

                     aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                >  1:2652771/1‑71 (MQ=255)

                     aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                >  1:5260611/1‑71 (MQ=255)

                     aaGAGCGCAACTTACTGAAGCAAAATAAAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                >  1:4828358/1‑71 (MQ=255)

                       gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcg                                                                  >  1:1735797/1‑51 (MQ=255)

                        agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                             >  1:4617417/1‑71 (MQ=255)

                         gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTAc                                            >  1:217300/1‑71 (MQ=255)

                            cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:989349/54‑1 (MQ=255)

                            cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:2017278/54‑1 (MQ=255)

                            cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                           >  1:2620716/1‑69 (MQ=255)

                             aaCTTACTGAAGCAAAATATAACCTTGAGTAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                   <  1:1854577/60‑1 (MQ=255)

                             aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                   <  1:2430738/60‑1 (MQ=255)

                             aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                   <  1:5539252/60‑1 (MQ=255)

                              aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCa                                                         <  1:2932552/53‑1 (MQ=255)

                              aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgt                                       <  1:1579372/71‑2 (MQ=255)

                              aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                       >  1:5504789/1‑71 (MQ=255)

                              aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                       >  1:766332/1‑71 (MQ=255)

                               cTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                        >  1:2418141/1‑53 (MQ=255)

                                ttACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:2555669/50‑1 (MQ=255)

                                  aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                        >  1:3004292/1‑50 (MQ=255)

                                  aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCATGTa                                     <  1:3341062/69‑1 (MQ=255)

                                   cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                     <  1:1242373/68‑1 (MQ=255)

                                    tGAAGCAAAATATAACCTTGAGCAGACTGTCATTcg                                                                    <  1:4180285/36‑1 (MQ=255)

                                    tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                      >  1:3481438/1‑50 (MQ=255)

                                     gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                                   <  1:311883/52‑3 (MQ=255)

                                     gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                                   <  1:4658593/52‑3 (MQ=255)

                                      aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                               >  1:5325897/1‑71 (MQ=255)

                                      aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                               >  1:2660324/1‑70 (MQ=255)

                                      aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                               >  1:5628654/1‑70 (MQ=255)

                                        gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                              <  1:1806351/38‑1 (MQ=255)

                                        gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                              <  1:3694242/38‑1 (MQ=255)

                                        gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                              <  1:4456269/38‑1 (MQ=255)

                                         cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                    <  1:5421041/47‑1 (MQ=255)

                                         cAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCg                             <  1:930966/70‑1 (MQ=255)

                                          aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGAc                                            >  1:4144447/1‑54 (MQ=255)

                                           aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGtt                                              <  1:3532462/51‑1 (MQ=255)

                                           aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                       <  1:2680923/42‑1 (MQ=12)

                                           aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                       <  1:1039785/42‑1 (MQ=12)

                                           aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                          >  1:1253642/1‑71 (MQ=255)

                                           aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                          <  1:3208415/71‑1 (MQ=255)

                                            aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:3718608/38‑1 (MQ=255)

                                            aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:3725302/38‑1 (MQ=255)

                                            aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                         <  1:1698509/71‑1 (MQ=255)

                                            aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTa                                     >  1:3456040/1‑59 (MQ=255)

                                            aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                         >  1:2880584/1‑70 (MQ=255)

                                            aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                         >  1:2733471/1‑70 (MQ=255)

                                            aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                         >  1:2610651/1‑70 (MQ=255)

                                             atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                        >  1:1906756/1‑71 (MQ=255)

                                              tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:629321/36‑1 (MQ=255)

                                              tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                          <  1:3082228/36‑1 (MQ=255)

                                              tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                               <  1:2890290/47‑1 (MQ=255)

                                              tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAgg                       >  1:3649775/1‑71 (MQ=255)

                                                taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                     >  1:1397815/1‑71 (MQ=255)

                                                taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                     <  1:3749721/71‑1 (MQ=255)

                                                taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                     >  1:4119278/1‑71 (MQ=255)

                                                taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTa                     >  1:5628176/1‑71 (MQ=255)

                                                  aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                   <  1:2085988/71‑1 (MQ=255)

                                                    cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                       <  1:449845/49‑1 (MQ=255)

                                                    cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                       <  1:2222215/49‑1 (MQ=255)

                                                    cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                       <  1:2274309/49‑1 (MQ=255)

                                                    cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa                 >  1:1271195/1‑71 (MQ=255)

                                                     ttGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGa                                 >  1:562770/1‑54 (MQ=255)

                                                      tGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg               >  1:2174325/1‑71 (MQ=255)

                                                       gAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                 >  1:2387129/1‑52 (MQ=255)

                                                        aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc              <  1:1028155/70‑1 (MQ=255)

                                                           aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt          <  1:428766/71‑1 (MQ=255)

                                                           aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt          <  1:3218372/71‑1 (MQ=255)

                                                            gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                            >  1:359309/1‑52 (MQ=255)

                                                            gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                   <  1:1560887/61‑1 (MQ=255)

                                                            gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg         >  1:3098634/1‑71 (MQ=255)

                                                            gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg         <  1:2670348/71‑1 (MQ=255)

                                                              cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc           >  1:5506930/1‑67 (MQ=255)

                                                              cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc           >  1:3899776/1‑67 (MQ=255)

                                                              cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc           >  1:917187/1‑67 (MQ=255)

                                                                 tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt     <  1:124460/70‑1 (MQ=255)

                                                                   aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc  >  1:2726988/1‑71 (MQ=255)

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCC  >  minE/2570692‑2570829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`