BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,964,217	G→A	22.6%	N159N (AAC→AAT)	lgt ←	phosphatidylglycerol‑prolipoprotein diacylglyceryl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,964,217	0	G	A	22.6%	60.0 / 13.3	31	N159N (AAC→AAT)	lgt	phosphatidylglycerol‑prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand): ref base G (16/8); new base A (5/2); total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGAC  >  W3110S.gb/2964157‑2964281

ggttTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAATTCGGGTc                                                          >  1:1473896/1‑69 (MQ=255)

   ttGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAAc                                                       >  1:1688403/1‑69 (MQ=255)

   ttGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAAc                                                       >  1:2919493/1‑69 (MQ=255)

   ttGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGCAACAGCATGGCAAACGGGAAGTTCGGGTCAAc                                                       >  1:1313960/1‑69 (MQ=255)

              aaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCACCGCGGCCCcaca                                            >  1:2045316/1‑69 (MQ=255)

              aaaTATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCcaca                                            >  1:71481/1‑69 (MQ=255)

                  atCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTc                                        <  1:1152700/69‑1 (MQ=255)

                   tcttctGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTc                                                               <  1:1648632/45‑1 (MQ=255)

                         gTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGt                                   <  1:3098380/67‑1 (MQ=255)

                                aGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGt                                   <  1:1092745/60‑1 (MQ=255)

                                        aaCAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGcc                      >  1:392812/1‑65 (MQ=255)

                                        aaCAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGcc                      >  1:3015506/1‑65 (MQ=255)

                                         aCAGCATGGCAAACGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATaaa                            <  1:1439072/58‑1 (MQ=255)

                                            gCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTtaataa                             >  1:2786541/1‑54 (MQ=255)

                                            gCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTtaataa                             >  1:142972/1‑54 (MQ=255)

                                             cATGGCAAACGGGAAGTTCGGGTCAACGCTGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    <  1:649357/62‑1 (MQ=255)

                                             cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    <  1:97398/62‑1 (MQ=255)

                                             cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    <  1:2899169/62‑1 (MQ=255)

                                             cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGccc             >  1:1704513/1‑69 (MQ=255)

                                             cATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGccc             >  1:251649/1‑69 (MQ=255)

                                             cATGGCAAACGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGccc             <  1:1427279/69‑1 (MQ=255)

                                               tGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTtaata                              >  1:2928728/1‑50 (MQ=255)

                                               tGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTtaat                               >  1:1807417/1‑49 (MQ=255)

                                                ggCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCgg           <  1:809149/68‑1 (MQ=255)

                                                 gCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    >  1:3332610/1‑58 (MQ=255)

                                                 gCAAACGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    >  1:413473/1‑58 (MQ=255)

                                                 gCAAACGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    >  1:858554/1‑58 (MQ=255)

                                                     aCGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCaa     >  1:718486/1‑69 (MQ=255)

                                                      cGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    >  1:2297999/1‑53 (MQ=255)

                                                      cGGGAAATTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCa                    >  1:2441490/1‑53 (MQ=255)

                                                         gAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGAc  >  1:3168536/1‑68 (MQ=255)

GGTTTGCAGCAGCAAAATATCTTCTGTACGGGAGCCAGGGAACAGCATGGCAAACGGGAAGTTCGGGTCAACGCGGCCCCACAATTCACCGTTAATAAAGTTGCCCAGACGCCCGGCACCAAGAC  >  W3110S.gb/2964157‑2964281

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