BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,571,124	C→T	14.6%	I343I (ATC→ATT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,571,124	0	C	T	14.6%	283.1 / 33.0	116	I343I (ATC→ATT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (41/58); new base T (14/3); total (55/61)
Fisher's exact test for biased strand distribution p-value = 2.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAA  >  Exported/2571054‑2571194

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:3013951/71‑1 (MQ=255)

gCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:4581699/71‑1 (MQ=255)

 cACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCt                                                                       <  1:810724/71‑1 (MQ=255)

    ggacggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:320691/67‑1 (MQ=255)

    ggacggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:3828242/67‑1 (MQ=255)

    ggacggTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:3972112/67‑1 (MQ=255)

          tgtgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCa                                                                <  1:3359755/69‑1 (MQ=255)

          tgtgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:1255654/1‑71 (MQ=255)

            tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGcc                                                                  <  1:4639024/65‑1 (MQ=255)

            tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGt                                                             <  1:2945067/70‑1 (MQ=255)

                gggAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCgg                                                        >  1:262481/1‑71 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      <  1:747772/71‑1 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      <  1:2570848/71‑1 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      <  1:1136149/71‑1 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      >  1:211201/1‑71 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      >  1:217817/1‑71 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      <  1:1229629/71‑1 (MQ=255)

                  gAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGCtt                                                      <  1:5537316/71‑1 (MQ=255)

                   aaCCATTGAACTGGACCCTAACGACGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGtt                                                      <  1:4378663/70‑1 (MQ=255)

                     ccATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAg                                                               >  1:4850571/1‑59 (MQ=255)

                     ccATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTAc                                                   >  1:2240849/1‑71 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                                >  1:2639266/1‑71 (MQ=255)

                        ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                                <  1:4493236/71‑1 (MQ=255)

                         tGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                                        <  1:2012615/46‑1 (MQ=255)

                         tGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCa                                                                >  1:2810396/1‑54 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTa                                                                      >  1:3791913/1‑47 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:2158537/1‑55 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:1955201/1‑55 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:5199416/1‑55 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:5111350/1‑55 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                              >  1:1043643/1‑55 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                                  >  1:3339208/1‑67 (MQ=255)

                          gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCg                                               <  1:4911662/70‑1 (MQ=255)

                              tGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTa                                                    <  1:1882072/61‑1 (MQ=255)

                              tGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                                  <  1:3219290/63‑1 (MQ=255)

                              tGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                                  <  1:5065744/63‑1 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                        <  1:2367287/71‑1 (MQ=255)

                                gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                        <  1:3539778/71‑1 (MQ=255)

                                   ccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                             <  1:3077075/63‑1 (MQ=255)

                                     tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGt                                                             >  1:4982549/1‑45 (MQ=255)

                                      aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                  >  1:880349/1‑71 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTa                                                    >  1:2660279/1‑52 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTc                                       <  1:53430/65‑1 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                 >  1:18021/1‑71 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                 >  1:4710826/1‑71 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                 >  1:589619/1‑71 (MQ=255)

                                       aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                 >  1:1023888/1‑71 (MQ=255)

                                           tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                             >  1:4241635/1‑71 (MQ=255)

                                             atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                  >  1:200360/1‑64 (MQ=255)

                                              tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                          <  1:3286039/71‑1 (MQ=255)

                                              tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                           <  1:5151667/70‑1 (MQ=255)

                                               atCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACt                                                  <  1:5679794/46‑1 (MQ=255)

                                                tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:3017175/70‑1 (MQ=255)

                                                tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:623872/70‑1 (MQ=255)

                                                 cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                       >  1:222200/1‑71 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                      >  1:4753995/1‑71 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                      >  1:3726211/1‑71 (MQ=255)

                                                  gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                      <  1:2769342/71‑1 (MQ=255)

                                                   aTGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATg                     >  1:5314680/1‑71 (MQ=255)

                                                    tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                             <  1:1158953/62‑1 (MQ=255)

                                                    tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                    >  1:4482994/1‑71 (MQ=255)

                                                      ctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                   >  1:2556208/3‑70 (MQ=255)

                                                      ctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                   >  1:2597122/3‑70 (MQ=255)

                                                      ctttACCAGACGGTATTTACGACCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCg                   >  1:4595646/3‑70 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:1480555/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:4756545/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:442177/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:2004031/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:1555635/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:3948811/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:2265140/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:3007160/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:2489517/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:1087588/71‑1 (MQ=255)

                                                       cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:3091865/71‑1 (MQ=255)

                                                        ttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 >  1:3570946/1‑68 (MQ=255)

                                                        ttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 >  1:61714/1‑68 (MQ=255)

                                                         tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 >  1:2958551/1‑67 (MQ=255)

                                                         tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 <  1:2708292/69‑3 (MQ=255)

                                                         tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 >  1:2442065/1‑67 (MQ=255)

                                                         tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                 <  1:4370586/69‑3 (MQ=255)

                                                          aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg              >  1:4422962/1‑71 (MQ=255)

                                                           ccagACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                      >  1:3602603/4‑62 (MQ=255)

                                                            ccGACGGCATCTACGCCCAGGTGTCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:3640913/58‑1 (MQ=255)

                                                            ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:101163/58‑1 (MQ=255)

                                                            ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:3453689/58‑1 (MQ=255)

                                                            ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                         <  1:4489800/58‑1 (MQ=255)

                                                             cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc           >  1:3225277/1‑71 (MQ=255)

                                                             agACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc           <  1:1571938/70‑1 (MQ=255)

                                                              gACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc           >  1:2934453/1‑70 (MQ=255)

                                                                cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat        >  1:5093293/1‑69 (MQ=255)

                                                                cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat        >  1:1634113/1‑69 (MQ=255)

                                                                cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                 <  1:1933922/62‑1 (MQ=255)

                                                                    aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct       >  1:177315/1‑68 (MQ=255)

                                                                    aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct       >  1:1785017/1‑68 (MQ=255)

                                                                     tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                      <  1:2807924/52‑1 (MQ=255)

                                                                     tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg    >  1:3844818/1‑70 (MQ=255)

                                                                      cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa               >  1:3105923/1‑58 (MQ=255)

                                                                      cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa  >  1:5194102/1‑71 (MQ=255)

                                                                      cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa  >  1:3046096/1‑71 (MQ=255)

GCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAA  >  Exported/2571054‑2571194

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`