BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,872,654	T→C	29.7%	E60E (GAA→GAG)	yjfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,872,654	0	T	C	29.7%	78.0 / 38.7	84	E60E (GAA→GAG)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (38/21); new base C (12/13); total (50/34)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGC  >  minE/2872589‑2872718

ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttc                                                                >  1:2068978/1‑68 (MQ=255)

ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttc                                                                >  1:2602144/1‑68 (MQ=255)

ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttc                                                                >  1:1497576/1‑68 (MQ=255)

ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttc                                                                >  1:2243733/1‑68 (MQ=255)

ttGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttc                                                                >  1:557969/1‑68 (MQ=255)

        aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                       >  1:978889/1‑69 (MQ=255)

        aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                       >  1:1840055/1‑69 (MQ=255)

        aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                       >  1:2973566/1‑69 (MQ=255)

        aCTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                       >  1:2644370/1‑69 (MQ=255)

               ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGa                                                <  1:1910827/69‑1 (MQ=255)

               ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGa                                                >  1:1573297/1‑69 (MQ=255)

               ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGa                                                <  1:1450629/69‑1 (MQ=255)

               ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGa                                                <  1:785040/69‑1 (MQ=255)

               ggTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGa                                                >  1:2244868/1‑69 (MQ=255)

                gTTCCCGCATCAGAATTACGTTGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                               >  1:2019678/1‑69 (MQ=255)

                gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                               >  1:2129101/1‑69 (MQ=255)

                gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                               >  1:260286/1‑69 (MQ=255)

                gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                               >  1:910494/1‑69 (MQ=255)

                gTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                               >  1:3177886/1‑69 (MQ=255)

                   cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTg                                            >  1:3276982/1‑69 (MQ=255)

                   cccGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTg                                            >  1:1013031/1‑69 (MQ=255)

                           aGATTTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:551314/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:2699303/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:2797695/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:2350969/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:2387105/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:3203445/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:367921/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:1047474/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:1772890/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                             <  1:759742/44‑1 (MQ=255)

                           aGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAttt                                                         <  1:868050/48‑1 (MQ=255)

                                 aCGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTg                              >  1:2717039/1‑69 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                      >  1:1341653/1‑69 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                      >  1:1337689/1‑69 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                      >  1:1242241/1‑69 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                      >  1:1190805/1‑69 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAt                       >  1:291542/1‑68 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:490493/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:2565638/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:2414603/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:928428/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:2725384/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:3077685/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:2150991/69‑1 (MQ=255)

                                         gtaCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:1191190/69‑1 (MQ=255)

                                         gtaCGTTGAGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:629691/69‑1 (MQ=255)

                                         gtaAGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                      >  1:93374/1‑69 (MQ=255)

                                          taCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:1647279/68‑1 (MQ=255)

                                          taCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGtt                      <  1:1963330/68‑1 (MQ=255)

                                          taCGTTGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCCCAGTGGCGAAGtt                      <  1:242846/68‑1 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTc                                  >  1:242706/1‑54 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGaaa                        >  1:3286930/1‑64 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAt                       >  1:510367/1‑65 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:3372267/1‑69 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:2110359/1‑69 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:2181839/1‑69 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:2709407/1‑69 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:1057561/1‑69 (MQ=255)

                                            cGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtg                   >  1:3265471/1‑69 (MQ=255)

                                              ttGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGTTAtga                  >  1:3006730/1‑68 (MQ=255)

                                              ttGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGTTAtga                  >  1:1045536/1‑68 (MQ=255)

                                              ttGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGTTAtga                  >  1:450233/1‑68 (MQ=255)

                                              ttGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGTTAtga                  >  1:2177438/1‑68 (MQ=255)

                                              ttGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGTGCTTCACTCTCAGTGGCGAAGTTAtga                  >  1:2597171/1‑68 (MQ=255)

                                               tGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgatt                <  1:863004/69‑1 (MQ=255)

                                               tGCGCCTTGCGCATCTGTCTCTTCATTTAATGCCTGACGAGCTTCACTCTCAGTGGCTAAGtt                      <  1:2232629/63‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:117831/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2320917/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:1428602/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:1226709/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:59680/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:1179133/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2133914/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2355656/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2585095/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2923173/68‑1 (MQ=255)

                                                         gcATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:3013352/68‑1 (MQ=255)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:3001534/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:469495/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:1286925/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:312256/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:2317325/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAATGCCTGACGTGCTTCACt                                     >  1:2199070/1‑35 (MQ=37)

                                                            tCTGTCTCTTCATTTAAAGCCTGACGTGCTTCACt                                     >  1:2772719/1‑35 (MQ=25)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:66405/63‑1 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:246409/63‑1 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:2720772/63‑1 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:1208835/63‑1 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt       <  1:1698789/63‑1 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:732857/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:1747726/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:2901691/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:1283187/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:2908568/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:2318403/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:1644384/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:944398/1‑68 (MQ=255)

                                                              tGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGc  >  1:3240994/1‑68 (MQ=255)

TTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCACGC  >  minE/2872589‑2872718

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