BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,779	C→T	20.0%	G228G (GGC→GGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,779	0	C	T	20.0%	217.0 / 43.8	105	G228G (GGC→GGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (45/39); new base T (17/4); total (62/43)
Fisher's exact test for biased strand distribution p-value = 2.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 114 total aligned reads displayed.

GCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCT  >  Exported/2570709‑2570849

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

gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                                        <  1:2443825/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)

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

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

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

   tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTAc                                                                     >  1:2671339/1‑71 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

                        cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                                        <  1:5421041/47‑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)

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

                          aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                              <  1:3208415/71‑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:2733471/1‑70 (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:4119278/1‑71 (MQ=255)

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

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

                                  ccTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgtc                                               >  1:3467639/1‑60 (MQ=255)

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

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

                                       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)

                                           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)

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

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

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

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

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

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

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

                                                     cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                     >  1:4819656/1‑69 (MQ=255)

                                                        gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                        >  1:974992/1‑47 (MQ=255)

                                                        gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                        >  1:3157193/1‑47 (MQ=255)

                                                        gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                        >  1:2029958/1‑47 (MQ=255)

                                                        gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                        >  1:2616414/1‑47 (MQ=255)

                                                        gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                        >  1:3886713/1‑47 (MQ=255)

                                                         cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt               >  1:2187208/1‑71 (MQ=255)

                                                         cgcCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGt               >  1:4854762/1‑71 (MQ=255)

                                                         cgcCGAGCCATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTatg                                   >  1:5221741/1‑49 (MQ=255)

                                                          gcagAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc              >  1:5095792/4‑71 (MQ=255)

                                                             gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg           <  1:4523482/71‑1 (MQ=255)

                                                             gAGCAATGGCTACGTTACTCAGGTACTGAACCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTcc             >  1:2462052/1‑69 (MQ=255)

                                                              aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                             >  1:3905850/1‑52 (MQ=255)

                                                              aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                             >  1:13360/1‑52 (MQ=255)

                                                              aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                             >  1:1658442/1‑52 (MQ=255)

                                                               gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg           <  1:3771478/69‑1 (MQ=255)

                                                               gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg           <  1:2791965/69‑1 (MQ=255)

                                                                  aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg      <  1:4950043/71‑1 (MQ=255)

                                                                  aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg      >  1:544441/1‑71 (MQ=255)

                                                                  aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc              <  1:2249279/63‑1 (MQ=255)

                                                                  aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc              <  1:5461260/63‑1 (MQ=255)

                                                                   tGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATgg     >  1:622311/1‑71 (MQ=255)

                                                                   tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa        <  1:198856/68‑1 (MQ=255)

                                                                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg            >  1:951436/1‑63 (MQ=255)

                                                                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg            >  1:4733665/1‑63 (MQ=255)

                                                                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg     >  1:2486290/1‑70 (MQ=255)

                                                                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt    >  1:5129572/1‑71 (MQ=255)

                                                                    ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt    >  1:4296935/1‑71 (MQ=255)

                                                                     gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc   >  1:3392954/1‑71 (MQ=255)

                                                                     gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc   <  1:2328269/71‑1 (MQ=255)

GCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCT  >  Exported/2570709‑2570849

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