BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,603,951	G→A	27.6%	D108D (GAC→GAT)	yfeY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,603,951	0	G	A	27.6%	169.2 / 61.0	148	D108D (GAC→GAT)	yfeY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (70/37); new base A (20/21); total (90/58)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATC  >  minE/1603885‑1604014

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACGGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:3399786/69‑1 (MQ=255)

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:2702937/69‑1 (MQ=255)

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:374339/69‑1 (MQ=255)

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:1439058/69‑1 (MQ=255)

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:2271909/69‑1 (MQ=255)

gTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:3280339/69‑1 (MQ=255)

      gTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCg                                                         >  1:2301926/1‑69 (MQ=255)

      gTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCg                                                         <  1:1938854/69‑1 (MQ=255)

      gTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCg                                                         >  1:1254677/1‑69 (MQ=255)

      gTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCg                                                         >  1:1239574/1‑69 (MQ=255)

         gCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               >  1:616853/1‑60 (MQ=255)

         gCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               >  1:396410/1‑60 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCGCGTCAATGCGACTGATCGTGCCTTGa                                         >  1:811067/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:752011/47‑1 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:2172013/47‑1 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:1541726/47‑1 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:555809/47‑1 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTg                                               >  1:2069968/1‑63 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:1733714/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:3244648/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:2115908/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:1583096/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:659192/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:2532251/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:251986/1‑69 (MQ=255)

                      ccGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGa                                         >  1:2485753/1‑69 (MQ=255)

                        gATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTc                                                               <  1:2571140/45‑1 (MQ=255)

                        gATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATc                                       >  1:2196020/1‑69 (MQ=255)

                        gATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATc                                       >  1:3189556/1‑69 (MQ=255)

                        gATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATc                                       >  1:2785584/1‑69 (MQ=255)

                                    ggTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTg                                                     <  1:1731932/43‑1 (MQ=255)

                                    ggTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTg                                                     <  1:327430/43‑1 (MQ=255)

                                    ggTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCAc                           <  1:785022/69‑1 (MQ=255)

                                    ggTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCAc                           <  1:1175281/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:2917077/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:293604/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:284889/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:3011421/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:2616261/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:2743228/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:361191/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:741205/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:944347/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:957411/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:170721/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:1266595/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:1628071/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:134760/69‑1 (MQ=255)

                                      tgtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:1557178/69‑1 (MQ=255)

                                       gtCAGCAGGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:497693/68‑1 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:2281262/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:3390481/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:1334233/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:2243966/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGcc                                             >  1:1458015/1‑45 (MQ=38)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:2800402/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:3312860/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:1634481/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:1968906/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:558602/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:3008369/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:767778/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:1571084/1‑69 (MQ=255)

                                          agcagGAATATCGCTATCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCATCGc                     >  1:814895/1‑69 (MQ=255)

                                                aaTATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGtt               <  1:728748/69‑1 (MQ=255)

                                                aaTATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGtt               <  1:1787900/69‑1 (MQ=255)

                                                aaTATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGtt               <  1:405025/69‑1 (MQ=255)

                                                aaTATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGtt               <  1:3366589/69‑1 (MQ=255)

                                                 atatCGCTATCCAGCACGGCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGtt               <  1:3362697/68‑1 (MQ=255)

                                                  tatCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGt             >  1:1350190/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTc            >  1:2548694/1‑68 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:1549424/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:3391403/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:3147862/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:3325003/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:3163277/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:1978483/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:2061750/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:2823636/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:2192366/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:2683892/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:262214/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:920328/1‑69 (MQ=255)

                                                    tCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCg           >  1:2280152/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:1094564/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:751736/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:2422519/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:1417939/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:299840/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:2975874/1‑69 (MQ=255)

                                                         aTCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCtt      >  1:2181128/1‑69 (MQ=255)

                                                         aTCCAGCACATCAATGCGGCTGATCGTGCCTTGATCGCCGTTAATCACCa                         <  1:914276/50‑1 (MQ=39)

                                                            cAGCACGTCAATGCTACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat   >  1:2278527/1‑69 (MQ=255)

                                                            cAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTGATCACCATCGCCACGTTGTCGCCCTtcat   >  1:1391436/1‑69 (MQ=255)

                                                            cAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat   >  1:47780/1‑69 (MQ=255)

                                                            cAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat   >  1:1581317/1‑69 (MQ=255)

                                                            cAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat   >  1:2132611/1‑69 (MQ=255)

                                                            cAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcat   >  1:2354599/1‑69 (MQ=255)

                                                              gCACGTCAATGCGACTGATCGTGCCTTGATCTCCGTTAATCACCATCGCCACGTTGTCGCCCTtcatc  <  1:1081412/68‑1 (MQ=255)

                                                              gCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTtcatc  <  1:2389229/68‑1 (MQ=255)

GTAAAGGTCGCTAAACGGTGTACCGATTTTAACACCGGTGTCAGCAGGAATATCGCTATCCAGCACGTCAATGCGACTGATCGTGCCTTGATCGCCGTTAATCACCATCGCCACGTTGTCGCCCTTCATC  >  minE/1603885‑1604014

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