BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,571,115	C→A	8.4%	P340P (CCC→CCA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,571,115	0	C	A	8.4%	369.3 / 9.2	116	P340P (CCC→CCA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (43/63); new base A (8/2); total (51/65)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGC  >  Exported/2571046‑2571185

cGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTAccc                                                                        >  1:2953177/1‑70 (MQ=255)

 gTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGa                                                                      <  1:1970903/71‑1 (MQ=255)

  ggCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGAGGATATCGATGCCTTACCCGa                                                                      <  1:4087652/69‑1 (MQ=255)

    cccGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACgg                                                                   >  1:591326/1‑71 (MQ=255)

     ccGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                                                  <  1:4589198/71‑1 (MQ=255)

     ccGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGc                                                                  <  1:3488973/71‑1 (MQ=255)

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

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

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

                 gtgtGCTGGGGACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCg                                                                       >  1:483428/1‑54 (MQ=255)

                 gtgtGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCg                                                                       >  1:515689/1‑54 (MQ=255)

                  tgtgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCAt                                                                >  1:5487429/1‑60 (MQ=255)

                  tgtgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCAt                                                                >  1:25842/1‑60 (MQ=255)

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

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

                    tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCa                                                                 <  1:2336986/57‑1 (MQ=255)

                    tgCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCa                                                                 <  1:3621813/57‑1 (MQ=255)

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

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

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

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

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

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

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

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

                            aCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTa                                           <  1:3724870/71‑1 (MQ=255)

                            aCCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTa                                           <  1:5313952/71‑1 (MQ=255)

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

                                ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATc                                                               <  1:4433594/47‑1 (MQ=255)

                                ttGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTcc                                       >  1:2639266/1‑71 (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:5111350/1‑55 (MQ=255)

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

                                  gAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAgg                                                     >  1:1955201/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)

                                   aaCTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                    <  1:2608988/71‑1 (MQ=255)

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

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

                                      tGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCa                                  >  1:4962606/1‑70 (MQ=255)

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

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

                                        gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                <  1:3070300/70‑1 (MQ=255)

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

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

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

                                               aCGATGATATCGATGCCTTACCCGACGGCATCTGCGCCCAgg                                                     >  1:5106830/1‑42 (MQ=255)

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

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

                                               aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                         <  1:62860/70‑1 (MQ=255)

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

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

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

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

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

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

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

                                                     atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                  >  1:3922071/1‑71 (MQ=255)

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

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

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

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

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

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

                                                          gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt             >  1:3172167/1‑71 (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:2597122/3‑70 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                               cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa        <  1:1087588/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:2442065/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

CGTGCCCGGCACGGACGGTGTGCTGGGAACCATTGAACTGGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGC  >  Exported/2571046‑2571185

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