BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,759,336	C→T	T45I (ACA→ATA)	yfjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,759,336	0	C	T	95.3%	258.2 / 2.9	86	T45I (ACA→ATA)	yfjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/4); new base T (38/44); total (38/48)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

ACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTT  >  W3110S.gb/2759273‑2759396

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:102255/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:433265/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:462585/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:6843/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:3917918/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:822906/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:2273056/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:3767863/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:2661703/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:3730595/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:3648655/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:907410/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:3225486/70‑1 (MQ=255)

aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                        <  1:310669/70‑1 (MQ=255)

  tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAatat                                                     >  1:2194579/1‑71 (MQ=255)

  tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAatat                                                     >  1:2189678/1‑71 (MQ=255)

    ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGATAATATg                                                    <  1:9057/70‑1 (MQ=255)

    ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                    <  1:3715831/70‑1 (MQ=255)

    ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                    <  1:3477890/70‑1 (MQ=255)

    ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                    <  1:116809/70‑1 (MQ=255)

    ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                    <  1:2026324/70‑1 (MQ=255)

        aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGAt                                                 <  1:3948152/69‑1 (MQ=255)

           gAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                          <  1:646830/57‑1 (MQ=255)

           gAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                          <  1:846750/57‑1 (MQ=255)

              cATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAtt                                         >  1:3994452/1‑71 (MQ=255)

              cATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAtt                                         >  1:3829577/1‑71 (MQ=255)

              cATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAtt                                         >  1:3802552/1‑71 (MQ=255)

              cATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAtt                                         >  1:841929/1‑71 (MQ=255)

              cATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAtt                                         >  1:846686/1‑71 (MQ=255)

                tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                        <  1:2199277/54‑1 (MQ=255)

                tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                        <  1:1697716/54‑1 (MQ=255)

                tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                        <  1:1402112/54‑1 (MQ=255)

                            gATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           >  1:3391179/1‑71 (MQ=255)

                            gATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           >  1:3573282/1‑71 (MQ=255)

                            gATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAACCATTCTTTCAACCCAggg                           >  1:1028713/1‑71 (MQ=255)

                             aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           <  1:2860424/70‑1 (MQ=255)

                             aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           <  1:1327570/70‑1 (MQ=255)

                             aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           <  1:3340224/70‑1 (MQ=255)

                             aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           <  1:1097085/70‑1 (MQ=255)

                             aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAggg                           <  1:1092321/70‑1 (MQ=255)

                              ttCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                    <  1:1974107/44‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:84927/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:3313485/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:3117457/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:1165623/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:1452423/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:2035359/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:2966170/65‑1 (MQ=255)

                                       ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:2554520/65‑1 (MQ=255)

                                        tGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                      <  1:2113189/64‑1 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:2314548/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:954942/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:1908434/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:2166479/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:3022343/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGActc                    >  1:2970408/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACACAGGGTTGActc                    >  1:2707050/1‑65 (MQ=255)

                                         gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACACAGGGTTGActc                    >  1:2702227/1‑65 (MQ=255)

                                          aTCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGc              <  1:355861/70‑1 (MQ=255)

                                          aTCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGc              <  1:1193302/70‑1 (MQ=255)

                                          aTCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGc              <  1:2217889/70‑1 (MQ=255)

                                               cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            <  1:2864031/67‑1 (MQ=255)

                                               cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            <  1:3741446/67‑1 (MQ=255)

                                               cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            <  1:1612345/67‑1 (MQ=255)

                                               cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            <  1:842179/67‑1 (MQ=255)

                                                aTGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            <  1:202530/66‑1 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:3427428/1‑67 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:324853/1‑67 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:298519/1‑67 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:2062128/1‑67 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:3433313/1‑67 (MQ=255)

                                                 tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtata          >  1:3710452/1‑67 (MQ=255)

                                                     ccAATTGATATAGAAAATATGGATAACCCATTCTTTCAAcc                                >  1:4112916/1‑41 (MQ=255)

                                                     ccAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCtt  >  1:3131102/1‑71 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:3972995/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:3954690/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:199806/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:1671378/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:3481180/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:2629306/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:3737869/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:1411942/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:2760651/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta            >  1:3075569/1‑56 (MQ=255)

                                                          tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCt             >  1:3999888/1‑55 (MQ=255)

                                                          agATATAGAAAATATGGATAACCCATTCTATCAAccc                               >  1:3494878/2‑37 (MQ=38)

ACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTT  >  W3110S.gb/2759273‑2759396

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