BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,393,339	A→G	10.1%	N30S (AAT→AGT)	yfbO →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,393,339	0	A	G	10.1%	261.1 / 19.3	99	N30S (AAT→AGT)	yfbO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (48/41); new base G (7/3); total (55/44)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 109 total aligned reads displayed.

TTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGA  >  W3110S.gb/2393270‑2393408

ttCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                       >  1:2520151/1‑70 (MQ=255)

   gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:1094110/71‑1 (MQ=255)

   gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:2927881/71‑1 (MQ=255)

   gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:2923193/71‑1 (MQ=255)

    tttGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:310644/70‑1 (MQ=255)

    tttGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCa                                                                  <  1:371363/71‑1 (MQ=255)

      tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 >  1:3315190/1‑70 (MQ=255)

      tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 >  1:617370/1‑70 (MQ=255)

      tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 >  1:135703/1‑70 (MQ=255)

      tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCaa                                                                 >  1:1047748/1‑70 (MQ=255)

      tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAAt                                                                >  1:3793331/1‑71 (MQ=255)

       gatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAAt                                                                <  1:1198176/70‑1 (MQ=255)

         tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                             <  1:1032394/71‑1 (MQ=255)

         tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                             <  1:647821/71‑1 (MQ=255)

         tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                             <  1:806931/71‑1 (MQ=255)

         tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                             <  1:3012290/71‑1 (MQ=255)

           tctatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGTTg                                                            <  1:718913/68‑1 (MQ=255)

           aaTATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                       <  1:1644705/59‑1 (MQ=255)

           aaTATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTg                                                            <  1:4003943/70‑1 (MQ=255)

             tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAata                                                                     >  1:1804787/1‑59 (MQ=255)

             tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:2316299/1‑69 (MQ=255)

             tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:2361113/1‑69 (MQ=255)

               tattatTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGt                                                        >  1:3497217/1‑70 (MQ=255)

               tattatTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGt                                                        >  1:192417/1‑70 (MQ=255)

                 ttattTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaa                                                     >  1:3568350/1‑71 (MQ=255)

                   attTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaat                                                    <  1:689232/70‑1 (MQ=255)

                       cAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATc                                                <  1:1111272/70‑1 (MQ=255)

                       cAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATc                                                <  1:3421091/70‑1 (MQ=255)

                        aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:233509/71‑1 (MQ=255)

                        aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:2311467/71‑1 (MQ=255)

                        aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:3864897/71‑1 (MQ=255)

                            aataatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATa                                          >  1:3998050/1‑71 (MQ=255)

                            aataatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATa                                          >  1:2379977/1‑71 (MQ=255)

                             ataatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:1377475/45‑1 (MQ=255)

                             ataatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                   <  1:3634370/45‑1 (MQ=255)

                              taatTATGACACCTTTAGAACTCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:3981492/1‑52 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 <  1:507027/46‑1 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 <  1:3799974/46‑1 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCaa                                                                 <  1:3893402/46‑1 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:2434559/1‑52 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:3357347/1‑52 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:2865959/1‑52 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATAct                                        >  1:1077149/1‑71 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATAct                                        >  1:870494/1‑71 (MQ=255)

                              taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATAct                                        >  1:3003572/1‑71 (MQ=255)

                               aatTATGACACTTTTAGAACGCATTACCCAACTTGTTaa                                                                       >  1:1910708/1‑39 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                       >  1:1931429/1‑39 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCaa                                                                 >  1:3260703/1‑45 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCaa                                                                 >  1:3135761/1‑45 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGAt                                                          >  1:778998/1‑52 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActc                                       >  1:641630/1‑71 (MQ=255)

                               aatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActc                                       >  1:141785/1‑71 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATg                                                               >  1:1152391/1‑46 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATg                                                               >  1:3416125/1‑46 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATg                                                               >  1:208345/1‑46 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATGGTGATGTTAATAATCCTGATActc                                       >  1:3746456/1‑70 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcc                                      >  1:2308497/1‑71 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcc                                      >  1:2999464/1‑71 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcc                                      >  1:1732599/1‑71 (MQ=255)

                                atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcc                                      >  1:1378269/1‑71 (MQ=255)

                                 ttATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                            <  1:381383/64‑1 (MQ=255)

                                 ttATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                            <  1:2211916/64‑1 (MQ=255)

                                  tATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaa                                                     >  1:586573/1‑54 (MQ=255)

                                   aTGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:3104581/1‑47 (MQ=255)

                                   aTGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCg                                   >  1:2010901/1‑71 (MQ=255)

                                          cTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:2421685/53‑1 (MQ=255)

                                          cTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:440132/53‑1 (MQ=255)

                                          cTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:1636374/53‑1 (MQ=255)

                                          cTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:2963101/53‑1 (MQ=255)

                                          cTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTg                            <  1:1748402/71‑1 (MQ=255)

                                           tttAGAACGCATTACCCAACTTTTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:3873024/52‑1 (MQ=255)

                                           tttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATc                                                <  1:2731452/50‑1 (MQ=255)

                                           tttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              <  1:3610901/52‑1 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                           >  1:1320879/1‑37 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaataa                                                  >  1:682558/1‑46 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaataa                                                  >  1:3033549/1‑46 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                              >  1:1114397/1‑50 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                            >  1:1134351/1‑52 (MQ=255)

                                             tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATa                                          >  1:433490/1‑54 (MQ=255)

                                               gAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTAt                        <  1:3641166/70‑1 (MQ=255)

                                                 aCGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTAtctc                     <  1:2508644/71‑1 (MQ=255)

                                                  cGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATctct                    >  1:2911998/1‑71 (MQ=255)

                                                  cGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATctct                    >  1:193397/1‑71 (MQ=255)

                                                       tACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTa                         >  1:916869/1‑61 (MQ=255)

                                                       tACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTa                         >  1:1216814/1‑61 (MQ=255)

                                                        aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGa              >  1:3829684/1‑71 (MQ=255)

                                                        aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGa              >  1:3319246/1‑71 (MQ=255)

                                                          ccAACTTGTTAATATCAATCGTGATGTTAATAATCCTGATACTCCTCGTCCATTGtt                          <  1:409863/57‑1 (MQ=255)

                                                           cAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAttt           >  1:142952/1‑71 (MQ=255)

                                                           cAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtt            >  1:1227918/1‑70 (MQ=255)

                                                              cTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt        <  1:3189034/71‑1 (MQ=255)

                                                              cTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt        <  1:156821/71‑1 (MQ=255)

                                                              cTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt        <  1:1407512/71‑1 (MQ=255)

                                                               ttGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt        >  1:3694031/1‑70 (MQ=255)

                                                                tGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTa      >  1:9103/1‑71 (MQ=255)

                                                                 gTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTa      <  1:1214627/70‑1 (MQ=255)

                                                                  ttAATATCAATTGTGATGTTAATAATCCTGATActcc                                      <  1:2952862/37‑1 (MQ=255)

                                                                    aaTATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGa  >  1:225331/1‑71 (MQ=255)

                                                                    aaTATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGa  >  1:1839168/1‑71 (MQ=255)

                                                                     atatCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTg   <  1:130008/69‑1 (MQ=255)

TTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGA  >  W3110S.gb/2393270‑2393408

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