BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,756	G→A	16.2%	V221I (GTC→ATC) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,756	0	G	A	16.2%	275.4 / 37.9	117	V221I (GTC→ATC) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (44/54); new base A (12/7); total (56/61)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 119 total aligned reads displayed.

CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTG  >  Exported/2570686‑2570822

cGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTg                                                                    <  1:2334190/71‑1 (MQ=255)

cGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTg                                                                    <  1:4627215/71‑1 (MQ=255)

 ggCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   >  1:2578877/1‑71 (MQ=255)

 ggCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   >  1:1530045/1‑71 (MQ=255)

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

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

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

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

          tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                          >  1:508860/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)

                    tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGATCAGACTGTCATTCGCGCGCCGAGCaa                                                 <  1:4146578/70‑1 (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:4609709/71‑1 (MQ=255)

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

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

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

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

                         ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   <  1:2248416/47‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                  cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:989349/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:766332/1‑71 (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:2963613/1‑50 (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)

                                           gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                            <  1:4658593/52‑3 (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:617446/38‑1 (MQ=255)

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

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

                                                 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)

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

                                                  aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3725302/38‑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:2733471/1‑70 (MQ=255)

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

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

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

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

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

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

                                                      taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa              >  1:4119278/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: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: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:917187/1‑67 (MQ=255)

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

CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTG  >  Exported/2570686‑2570822

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