BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,855,194	A→T	35.5%	I90N (ATC→AAC)	ygbA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,855,194	0	A	T	35.5%	42.5 / 48.2	48	I90N (ATC→AAC)	ygbA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (21/10); new base T (4/13); total (25/23)
Fisher's exact test for biased strand distribution p-value = 5.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

ACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTG  >  W3110S.gb/2855130‑2855260

aCGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCg                                                                 <  1:793615/68‑1 (MQ=255)

        tttCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCg                                                                 >  1:3265564/1‑60 (MQ=255)

            cGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                                     <  1:3097175/68‑1 (MQ=255)

             ggCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                                     <  1:2835584/67‑1 (MQ=255)

              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 <  1:862992/54‑1 (MQ=255)

              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 <  1:2673838/54‑1 (MQ=255)

              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:298684/1‑54 (MQ=255)

              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:430564/1‑54 (MQ=255)

              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:1778694/1‑54 (MQ=255)

                aGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                                                <  1:816171/69‑1 (MQ=255)

                aGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                                                <  1:3261248/69‑1 (MQ=255)

                   tCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtctc                                             >  1:1616848/1‑69 (MQ=255)

                            ggACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccggc                                     >  1:621967/1‑68 (MQ=255)

                             gACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                                              <  1:3480415/42‑1 (MQ=255)

                              aCGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCa                                  >  1:871509/1‑69 (MQ=255)

                                   tGTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCa                             <  1:2637601/69‑1 (MQ=255)

                                   tGTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCa                             <  1:2845198/69‑1 (MQ=255)

                                    gTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCggcccggccc                                   >  1:1832626/1‑62 (MQ=255)

                                      cATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCAt                            <  1:2895966/67‑1 (MQ=255)

                                      cATCGATGAGATGACGCACGGTTAAGATCGGATGCCGCCACAGCATTCTCGGCCCGGCCCAACGCATaa                          >  1:131253/1‑69 (MQ=255)

                                       aTCGATGAGATGACGCACGGTTAAGATCGGATGACGGCACAGCATTCTCggcccggccc                                   >  1:429238/1‑59 (MQ=255)

                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:370288/1‑63 (MQ=255)

                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:371049/1‑63 (MQ=255)

                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:2736970/1‑63 (MQ=255)

                                             gagatcaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:1326228/1‑63 (MQ=255)

                                              agatgaCGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2828775/69‑1 (MQ=255)

                                              agatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2784073/69‑1 (MQ=255)

                                              agatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:1210727/69‑1 (MQ=255)

                                               gatgaCGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2026483/68‑1 (MQ=255)

                                                   aCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCAtctc             >  1:1198865/1‑69 (MQ=255)

                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:2618488/69‑1 (MQ=255)

                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1349684/69‑1 (MQ=255)

                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:790530/69‑1 (MQ=255)

                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1439811/69‑1 (MQ=255)

                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1548255/69‑1 (MQ=255)

                                                        cGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  >  1:2636434/1‑59 (MQ=255)

                                                          gTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      <  1:44252/69‑1 (MQ=255)

                                                          gTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc                                   <  1:1847964/40‑1 (MQ=255)

                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:1687835/1‑69 (MQ=255)

                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:729302/1‑69 (MQ=255)

                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:733577/1‑69 (MQ=255)

                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:2561935/1‑69 (MQ=255)

                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGc                              >  1:602898/1‑41 (MQ=255)

                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:525077/1‑67 (MQ=255)

                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:1367089/1‑67 (MQ=255)

                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:3134851/1‑67 (MQ=255)

                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:2774758/1‑67 (MQ=255)

                                                               gATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTg  <  1:2517998/68‑1 (MQ=255)

                                                                aTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCttt   >  1:2556003/1‑66 (MQ=255)

                                                                aTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTg  <  1:1212662/67‑1 (MQ=255)

ACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTG  >  W3110S.gb/2855130‑2855260

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