BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,764	C→T	17.2%	R223R (CGC→CGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,764	0	C	T	17.2%	248.7 / 38.6	110	R223R (CGC→CGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (34/57); new base T (12/7); total (46/64)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 119 total aligned reads displayed.

TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGC  >  Exported/2570696‑2570831

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

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

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

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

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

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

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

        tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                            >  1:3519118/1‑70 (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:1826526/71‑1 (MQ=255)

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

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

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

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

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

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

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

                 aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                                  >  1:2652771/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)

                     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)

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

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

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

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

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

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

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

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

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

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

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

                                 gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                                     <  1:311883/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)

                                    gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                                <  1:617446/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)

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

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

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

                                       aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                            >  1:1253642/1‑71 (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)

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

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

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

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

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

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

                                            taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                       <  1:3749721/71‑1 (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:449845/49‑1 (MQ=255)

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

                                                cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                         <  1:2274309/49‑1 (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:3218372/71‑1 (MQ=255)

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

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

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

                                                        gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                     <  1:1560887/61‑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:917187/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGC  >  Exported/2570696‑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: `‑`