BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,477,935	G→A	81.0%	E360K (GAA→AAA)	ydbD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,477,935	0	G	A	81.0%	127.6 / 24.9	58	E360K (GAA→AAA)	ydbD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (4/7); new base A (20/27); total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGG  >  W3110S.gb/1477872‑1477998

gAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                            <  1:1762027/69‑1 (MQ=255)

gAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                            <  1:3191800/69‑1 (MQ=255)

  aataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtct                                                          <  1:3570077/69‑1 (MQ=255)

  aataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtct                                                          <  1:2420346/69‑1 (MQ=255)

    taatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTcc                                                        <  1:2476202/69‑1 (MQ=255)

     aatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTccc                                                       >  1:1878561/1‑69 (MQ=255)

     aatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTccc                                                       >  1:461266/1‑69 (MQ=255)

        gACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtctc                                                         <  1:1614974/64‑1 (MQ=255)

        gACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtctc                                                         <  1:1076950/64‑1 (MQ=255)

        gACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtctc                                                         <  1:79605/64‑1 (MQ=255)

           gcgAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGa                                                 <  1:2361974/69‑1 (MQ=255)

                 gAGTCGTCTTCCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                           <  1:3272275/69‑1 (MQ=255)

                 gAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                           <  1:1510230/69‑1 (MQ=255)

                 gAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                           <  1:866647/69‑1 (MQ=255)

                 gAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                           <  1:539943/69‑1 (MQ=255)

                  aGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                           <  1:818351/68‑1 (MQ=255)

                        cTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGAct                                    >  1:3435013/1‑69 (MQ=255)

                         tACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGActc                                   >  1:943153/1‑69 (MQ=255)

                         tACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGActc                                   >  1:3163868/1‑69 (MQ=255)

                         tACCTGATGGAGTCGGTCCGTTGAATCTTGAGCTACGTAAAAATCTcc                                                        >  1:3191609/1‑48 (MQ=255)

                              gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGc                               >  1:862468/1‑68 (MQ=255)

                                 ggTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            <  1:872428/68‑1 (MQ=255)

                                    gtCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAaga                        <  1:2236358/69‑1 (MQ=255)

                                    gtCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAaga                        <  1:1671198/69‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTg                                                   <  1:3199862/39‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTg                                                   <  1:19342/39‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGa                                                 <  1:694678/41‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGa                                                 <  1:2632521/41‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGa                                                 <  1:2800758/41‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTg                                 <  1:1406667/57‑1 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGAt                     >  1:620541/1‑69 (MQ=255)

                                       ggTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGAt                     >  1:3674372/1‑69 (MQ=255)

                                          ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTgg                  <  1:790287/69‑1 (MQ=255)

                                          ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTgg                  <  1:199986/69‑1 (MQ=255)

                                          ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATAGTCATGACTCTGGCGAAAAGAGATTgg                  <  1:1994545/69‑1 (MQ=255)

                                           cGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTgg                  >  1:2988066/1‑68 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            >  1:598800/1‑54 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            >  1:3114900/1‑54 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            >  1:199788/1‑54 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGa                             >  1:2983689/1‑53 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGatat             >  1:2910399/1‑69 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGatat             >  1:2068331/1‑69 (MQ=255)

                                               gAATCTTGCGCTACGTAAAAATCTCCCTCTCGATATCGTCATGACTCTGGCGAAAAgaga                      >  1:817198/1‑60 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            <  1:3612403/49‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            <  1:3100576/49‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            <  1:1531244/49‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGATATTGGGa                <  1:3728783/61‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                <  1:3133421/61‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                <  1:2978500/61‑1 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                >  1:2553280/1‑61 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                >  1:2166053/1‑61 (MQ=255)

                                                    ttGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                >  1:2076416/1‑61 (MQ=255)

                                                     tGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa        >  1:3290/1‑68 (MQ=255)

                                                     tGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                            <  1:3593442/48‑1 (MQ=255)

                                                     tGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa        >  1:852935/1‑68 (MQ=255)

                                                     tGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa        >  1:569486/1‑68 (MQ=255)

                                                     tGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa        >  1:1347309/1‑68 (MQ=255)

                                                          tACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAgg  <  1:1473973/69‑1 (MQ=255)

GAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGG  >  W3110S.gb/1477872‑1477998

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