BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,438,124	T→C	W140R (TGG→CGG)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,438,124	0	T	C	92.7%	427.4 / 18.2	124	W140R (TGG→CGG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand): ref base T (1/8); new base C (48/67); total (49/75)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 124 total aligned reads displayed.

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAA  >  NC_000913/4438078‑4438171

cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTgg                                               <  1:3123294/49‑1 (MQ=255)

cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                               >  1:3699672/1‑49 (MQ=255)

 aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGa                                             <  2:5114005/50‑1 (MQ=255)

 aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                             >  2:213836/1‑50 (MQ=255)

  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                            <  1:3358572/50‑1 (MQ=255)

  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                            <  1:1265739/50‑1 (MQ=255)

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                           <  1:986463/50‑1 (MQ=255)

   aaaCCGTGAATGTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGata                                           <  2:4357505/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                          <  2:2831348/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                          >  2:2313423/1‑50 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                          <  1:4840112/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                          <  2:1688297/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                          <  1:73145/50‑1 (MQ=255)

     aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATAtt                                         <  2:3432860/50‑1 (MQ=255)

     aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                         >  1:3310523/1‑50 (MQ=255)

     aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                         >  1:2187703/1‑50 (MQ=255)

      ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                        <  2:2221678/50‑1 (MQ=255)

       cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                        >  2:2393443/1‑49 (MQ=255)

       cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                        >  2:878405/1‑49 (MQ=255)

       cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                       <  2:4890347/50‑1 (MQ=255)

         tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                     <  1:1286296/50‑1 (MQ=255)

         tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                     <  2:1688012/50‑1 (MQ=255)

          gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                      >  2:3191303/1‑48 (MQ=255)

           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                    <  2:4375369/49‑1 (MQ=255)

           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                    <  1:5362531/49‑1 (MQ=255)

           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                   <  2:3021194/50‑1 (MQ=255)

            aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                  <  2:1115774/50‑1 (MQ=255)

            aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                  <  1:5270905/50‑1 (MQ=255)

             tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                                 <  1:2701633/50‑1 (MQ=255)

               tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                                <  1:2369726/49‑1 (MQ=255)

                ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                               >  1:5279660/1‑49 (MQ=255)

                ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                               >  1:3594682/1‑49 (MQ=255)

                 cGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGaaaaa                              <  1:356355/49‑1 (MQ=255)

                 cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                             <  1:643520/50‑1 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAAc                            <  1:2837107/50‑1 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                             >  1:4345116/1‑49 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                             >  1:4434208/1‑49 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                            <  2:4006319/50‑1 (MQ=255)

                   ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                           <  2:4082724/50‑1 (MQ=255)

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                          <  2:4555/49‑1 (MQ=255)

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                          <  2:3074468/49‑1 (MQ=255)

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                         >  1:4359743/1‑50 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                         <  1:1559610/49‑1 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                         <  1:3164601/49‑1 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                        >  1:1208544/1‑50 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                        >  2:2587465/1‑50 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                        >  2:641982/1‑50 (MQ=255)

                      gAACAAAACAGAAGATTGTATGCCCTGGATATTCTACAGAAAAAACAGTc                        <  2:2421898/50‑1 (MQ=255)

                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                       >  2:57945/1‑50 (MQ=255)

                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                       >  2:2478357/1‑49 (MQ=255)

                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                      <  1:110200/50‑1 (MQ=255)

                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                      <  1:1595090/50‑1 (MQ=255)

                         cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                     >  2:3722735/1‑50 (MQ=255)

                         cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                     <  1:3517932/50‑1 (MQ=255)

                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                     >  2:1444197/1‑49 (MQ=255)

                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                     >  2:3502961/1‑49 (MQ=255)

                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                     >  1:1590294/1‑49 (MQ=255)

                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                    <  2:1244377/50‑1 (MQ=255)

                           aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                   <  2:1794802/50‑1 (MQ=255)

                           aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                   <  2:2332773/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa                 <  2:4301872/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa                 <  1:4072054/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa                 >  1:418599/1‑50 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa                 <  2:218297/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa                 <  2:584940/50‑1 (MQ=255)

                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc                >  1:2338686/1‑50 (MQ=255)

                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc                <  1:2799732/50‑1 (MQ=255)

                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               >  1:2714016/1‑50 (MQ=255)

                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               <  1:1972702/50‑1 (MQ=255)

                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               >  1:1812266/1‑50 (MQ=255)

                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               <  1:1499057/50‑1 (MQ=255)

                               aGACGATTGTATGCCCGAGATATTCTACAGAAAAAGCAGTCTTACATGcc               <  2:4091862/50‑1 (MQ=25)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               >  1:3042457/1‑49 (MQ=255)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc               >  1:2489242/1‑49 (MQ=255)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg              <  2:2248327/50‑1 (MQ=255)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg              <  2:4990019/50‑1 (MQ=255)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg              <  1:2113341/50‑1 (MQ=255)

                                 aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa             <  2:3187692/50‑1 (MQ=255)

                                 aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa             <  1:943800/50‑1 (MQ=255)

                                  cGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa             >  1:1236691/1‑49 (MQ=255)

                                  cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa            <  1:2113149/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa           <  2:2897881/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa           <  2:1895177/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa            >  1:5075214/1‑49 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGCCTGACATACCGaaa           <  1:638427/50‑1 (MQ=38)

                                    aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa          <  2:2900299/50‑1 (MQ=255)

                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg         >  2:3840264/1‑50 (MQ=255)

                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg         <  1:5310864/50‑1 (MQ=255)

                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg         >  1:1172452/1‑50 (MQ=255)

                                       gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg        <  1:3957358/49‑1 (MQ=255)

                                       gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa       >  2:3317684/1‑50 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa       <  2:5301247/49‑1 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa       <  1:1214102/49‑1 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt      >  2:3199972/1‑50 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt      <  2:2862197/50‑1 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt      >  1:4291359/1‑50 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt      >  1:1243467/1‑50 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt      <  2:2009155/50‑1 (MQ=255)

                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt     <  2:2084838/49‑1 (MQ=255)

                                            ccTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa  <  1:1035102/50‑1 (MQ=255)

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAA  >  NC_000913/4438078‑4438171

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