BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,761	T→G	16.8%	I222M (ATT→ATG) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,761	0	T	G	16.8%	216.9 / 49.4	113	I222M (ATT→ATG) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (40/54); new base G (12/7); total (52/61)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 123 total aligned reads displayed.

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGC  >  minE/2570692‑2570831

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

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

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

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

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

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

      gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                                 >  1:3657184/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)

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

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

                aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                        <  1:5068480/70‑1 (MQ=255)

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

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

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

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

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

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

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

                    tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTAc                                                   >  1:2671339/1‑71 (MQ=255)

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

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

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

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

                      agagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                                 <  1:5429040/71‑1 (MQ=255)

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

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

                            cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                            <  1:989349/54‑1 (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)

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

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

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

                                  aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                     >  1:129133/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                          aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgt                              <  1:212453/70‑2 (MQ=255)

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

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

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

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

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

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

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

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

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

                                            aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                           >  1:2733471/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:3749721/71‑1 (MQ=255)

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

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

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

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

                                                    cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                         <  1:2222215/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)

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

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

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

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

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

                                                            gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg           >  1:3098634/1‑71 (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)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:2016074/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:4652072/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:3065483/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:1277911/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:3667828/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:96965/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:3303062/71‑1 (MQ=255)

                                                                    ttCGCGCGCCGAGCAAAGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg   <  1:1593534/71‑1 (MQ=255)

                                                                     tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc  >  1:4405786/1‑71 (MQ=255)

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGC  >  minE/2570692‑2570831

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