BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,357,330	G→A	5.3%	V23M (GTG→ATG)	yjdK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,357,330	0	G	A	5.3%	438.1 / 8.8	132	V23M (GTG→ATG)	yjdK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (71/54); new base A (5/2); total (76/56)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 139 total aligned reads displayed.

ATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGAT  >  W3110S.gb/4357261‑4357400

attatGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAgt                                                                       >  1:3052214/1‑71 (MQ=255)

   atGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAgtgt                                                                     >  1:1835305/1‑70 (MQ=255)

   atGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAgtgt                                                                     <  1:532000/70‑1 (MQ=255)

    tGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTc                                                                   >  1:5124297/1‑71 (MQ=255)

           tAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                            >  1:6943109/1‑71 (MQ=255)

           tAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                            <  1:5314986/71‑1 (MQ=255)

            aaaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                            >  1:1726273/1‑70 (MQ=255)

             aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                          <  1:5116879/71‑1 (MQ=255)

             aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                          <  1:4168918/71‑1 (MQ=255)

             aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                           <  1:350967/70‑1 (MQ=255)

              aaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                           >  1:1577165/1‑69 (MQ=255)

              aaaCAAGTTCAATACTTATGTTGTTTCTTTTAATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                           >  1:3433550/1‑69 (MQ=255)

               aaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagattaag                                                        <  1:7346574/71‑1 (MQ=255)

               aaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCtt                                                                 <  1:451043/62‑1 (MQ=255)

                aCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                          <  1:5936349/68‑1 (MQ=255)

                 cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                           >  1:1008557/1‑66 (MQ=255)

                 cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                            >  1:3604005/1‑65 (MQ=255)

                 cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTtaagattaagat                                                      <  1:2739540/71‑1 (MQ=255)

                 cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATATCATACATCAATGATCttaatcttaagaa                                                      <  1:4302477/71‑2 (MQ=17)

                  aaGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTtaagattaagat                                                      >  1:3599066/1‑70 (MQ=255)

                   aGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                    >  1:2419598/1‑71 (MQ=255)

                   aGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                    <  1:4620995/71‑1 (MQ=255)

                   aGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                    >  1:5653615/1‑71 (MQ=255)

                    gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                   >  1:4873046/1‑71 (MQ=255)

                    gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                   >  1:7438407/1‑71 (MQ=255)

                    gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                   >  1:7447532/1‑71 (MQ=255)

                       cAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGa                                                >  1:4256943/1‑71 (MQ=255)

                        aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagattaa                                                         >  1:5582595/1‑61 (MQ=255)

                        aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGa                                                >  1:1814116/1‑70 (MQ=255)

                        aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGAt                                               >  1:4211756/1‑71 (MQ=255)

                        aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGAt                                               <  1:4178858/71‑1 (MQ=255)

                         aTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGa                                                >  1:3927905/1‑69 (MQ=255)

                         aTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATg                                              <  1:5331598/71‑1 (MQ=255)

                            cTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           <  1:5027524/71‑1 (MQ=255)

                             ttATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           >  1:2446765/1‑70 (MQ=255)

                             ttATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           >  1:7780165/1‑70 (MQ=255)

                               aTGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatga                                         <  1:2054130/70‑1 (MQ=255)

                               aTGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatga                                         <  1:574675/70‑1 (MQ=255)

                                 gttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTGTTAAGATTAAGATCa                                                    >  1:4364414/1‑57 (MQ=255)

                                  ttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgatgt                                      <  1:7076172/70‑3 (MQ=255)

                                  ttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgatat                                      <  1:5499642/70‑1 (MQ=255)

                                    gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgat                                        <  1:6679447/66‑1 (MQ=255)

                                    gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGaa                                   >  1:7554384/1‑71 (MQ=255)

                                    gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGaa                                   >  1:7557179/1‑71 (MQ=255)

                                      ttCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                  <  1:2330545/70‑1 (MQ=255)

                                        cTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           <  1:7379275/59‑1 (MQ=255)

                                         ttttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           >  1:5440388/1‑58 (MQ=255)

                                          tttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGAt                                               >  1:6130025/1‑53 (MQ=255)

                                          tttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                             >  1:6785309/1‑71 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCTtaagattaagat                                                      >  1:7197851/1‑45 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           >  1:2876694/1‑56 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                           >  1:1054335/1‑56 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtc                               <  1:3926448/68‑1 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtc                               <  1:1365705/68‑1 (MQ=255)

                                           ttGATTATCCATCATCTTATTCCTCAGTGTTCCTAAGATTAAGATCATTGATGTATGATATGAATTTctcc                            >  1:6893408/1‑71 (MQ=255)

                                            tGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAttt                                >  1:7894706/1‑66 (MQ=255)

                                            tGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTctcc                            <  1:2905483/70‑1 (MQ=255)

                                               ttATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAttt                                >  1:7874406/1‑63 (MQ=255)

                                               ttATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTa                        >  1:2451424/1‑71 (MQ=255)

                                               ttATCCATCATCTTATTCCTCAATGATCttaatcttaagaac                                                     <  1:233497/42‑3 (MQ=14)

                                                 aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTAt                       <  1:6327128/70‑1 (MQ=255)

                                                 aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTAt                       <  1:770864/70‑1 (MQ=255)

                                                 aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATc                      >  1:4003922/1‑71 (MQ=255)

                                                 aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                   >  1:5480597/1‑41 (MQ=2)

                                                 aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                   >  1:6793317/1‑41 (MQ=2)

                                                 aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                   >  1:3747075/1‑41 (MQ=2)

                                                  tCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                     <  1:2443517/71‑1 (MQ=255)

                                                    catcatCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGt                    <  1:349817/70‑1 (MQ=255)

                                                    catcatCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTg                   <  1:5480890/71‑1 (MQ=255)

                                                      tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                    <  1:5980845/36‑1 (MQ=255)

                                                      tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                  <  1:5413538/54‑1 (MQ=255)

                                                      tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                             <  1:36111/59‑1 (MQ=255)

                                                      tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                             <  1:7923413/59‑1 (MQ=255)

                                                       catcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                >  1:5813902/1‑71 (MQ=255)

                                                        atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTctcctc                          >  1:97069/1‑60 (MQ=255)

                                                        atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtg               <  1:6261208/71‑1 (MQ=255)

                                                        atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                >  1:2838986/1‑70 (MQ=255)

                                                        atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                <  1:6738139/70‑1 (MQ=255)

                                                        atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                >  1:2059021/1‑70 (MQ=255)

                                                         tcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTg                   >  1:6322325/1‑66 (MQ=255)

                                                         tcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTg                   >  1:5975218/1‑66 (MQ=255)

                                                          cTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGt                    <  1:6889891/64‑1 (MQ=255)

                                                           ttATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCctct                         >  1:6117283/1‑58 (MQ=255)

                                                           ttATACCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtgat             >  1:3909232/1‑70 (MQ=255)

                                                            tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATTTGAAt                                  >  1:6611989/1‑48 (MQ=255)

                                                            tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                  >  1:1919660/1‑48 (MQ=255)

                                                            tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTgg                  >  1:275983/1‑64 (MQ=255)

                                                            tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTgg                  >  1:7674056/1‑64 (MQ=255)

                                                             aTTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                     <  1:883418/60‑1 (MQ=255)

                                                               acatCAATGATCTTAATCTTAAGAACACTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata        >  1:7105651/4‑71 (MQ=18)

                                                                ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata        <  1:2408991/70‑1 (MQ=255)

                                                                ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata        <  1:4156986/70‑1 (MQ=255)

                                                                ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata        <  1:1002387/70‑1 (MQ=255)

                                                                ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata        <  1:5629002/70‑1 (MQ=255)

                                                                   cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                              >  1:2135179/1‑45 (MQ=255)

                                                                   cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATg      >  1:1854131/1‑69 (MQ=255)

                                                                    aGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATgg     <  1:3108018/69‑1 (MQ=255)

                                                                     gtgtTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                              <  1:3645466/43‑1 (MQ=255)

                                                                     gtgtTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                              <  1:5032439/43‑1 (MQ=255)

                                                                     gtgtTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGAt  <  1:1788338/71‑1 (MQ=255)

ATTATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGAT  >  W3110S.gb/4357261‑4357400

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