BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,570,758	C→A	16.5%	V221V (GTC→GTA) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,570,758	0	C	A	16.5%	322.8 / 32.3	115	V221V (GTC→GTA) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (42/54); new base A (12/7); total (54/61)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 115 total aligned reads displayed.

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTT  >  Exported/2570692‑2570826

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

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

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

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

                 gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                 <  1:1826526/71‑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:2380771/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)

                     aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                     >  1:413168/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)

                     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)

                        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)

                             aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                                <  1:5539252/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: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: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)

                                    tGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                   >  1:3481438/1‑50 (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:4456269/38‑1 (MQ=255)

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

                                        gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                           <  1:1806351/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: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: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:2880584/1‑70 (MQ=255)

                                             atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                     >  1:1906756/1‑71 (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: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: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)

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

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

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

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

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

GCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTT  >  Exported/2570692‑2570826

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