BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,562,802	G→A	6.5%	S425S (TCC→TCT)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,562,802	0	G	A	6.5%	365.5 / 9.5	122	S425S (TCC→TCT)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand): ref base G (67/47); new base A (3/5); total (70/52)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 125 total aligned reads displayed.

GAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAG  >  Exported/2562733‑2562873

gAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCgg                                                                        <  1:5551658/71‑1 (MQ=255)

  aCGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGaa                                                                      <  1:1743839/71‑1 (MQ=255)

  aCGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGaa                                                                      <  1:1619930/71‑1 (MQ=255)

  aCGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGaa                                                                      >  1:3912729/1‑71 (MQ=255)

   cGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAAt                                                                     <  1:4564681/71‑1 (MQ=255)

       gagaAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGtt                                                                 >  1:5157284/1‑71 (MQ=255)

        agaagaAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTc                                                                <  1:3103375/71‑1 (MQ=255)

          aagaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGa                                                              <  1:1586393/71‑1 (MQ=255)

           agaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCg                                                               >  1:4599888/1‑69 (MQ=255)

           agaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGa                                                              >  1:4212015/1‑70 (MQ=255)

            gaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:3084253/1‑71 (MQ=255)

            gaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:290423/1‑71 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:5415613/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:1409973/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:4212770/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:3844628/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:2822310/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:2597269/71‑1 (MQ=255)

                    aTGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                    <  1:835765/71‑1 (MQ=255)

                      gCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATaa                                                         <  1:2481695/64‑1 (MQ=255)

                       cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAAt                                                 >  1:3893123/1‑71 (MQ=255)

                       cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAAt                                                 >  1:5169307/1‑71 (MQ=255)

                       cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAAt                                                 >  1:106148/1‑71 (MQ=255)

                           aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaa                                                             >  1:5344643/1‑55 (MQ=255)

                           aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaa                                                             >  1:593341/1‑55 (MQ=255)

                           aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAg                                                       >  1:3647968/1‑61 (MQ=255)

                           aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGc                                               >  1:4808018/1‑69 (MQ=255)

                            ccGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGc                                                     >  1:5466577/1‑62 (MQ=255)

                            ccGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAg                                            <  1:4606298/71‑1 (MQ=255)

                            ccGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAg                                            >  1:1020275/1‑71 (MQ=255)

                             cGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATa                                                          <  1:4622033/56‑1 (MQ=255)

                             cGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATa                                                          <  1:5176153/56‑1 (MQ=255)

                                tgtgCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCg                                                               <  1:1141503/48‑1 (MQ=255)

                                  tgCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGa                                                              <  1:5610371/47‑1 (MQ=255)

                                  tgCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGa                                                              <  1:1259051/47‑1 (MQ=255)

                                  tgCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATaaa                                                        <  1:4947671/53‑1 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:2550696/1‑47 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:4698748/1‑47 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:1587077/1‑47 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGaaa                                                            >  1:1418383/1‑47 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAAt                                                           <  1:5028110/48‑1 (MQ=255)

                                    cgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATAc                                    >  1:5556386/1‑71 (MQ=255)

                                      ccAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATAcc                                   >  1:4946521/1‑70 (MQ=255)

                                      ccAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATAcc                                   >  1:3398620/1‑70 (MQ=255)

                                      ccAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCg                                  <  1:4700014/71‑1 (MQ=255)

                                       cAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGccaa                                 <  1:4206624/71‑3 (MQ=255)

                                       cAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGccaa                                 <  1:5280274/71‑3 (MQ=255)

                                        aGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGAt                                >  1:2159907/1‑71 (MQ=255)

                                         gTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                               <  1:1936672/71‑1 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATa                                                          <  1:4429701/43‑1 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATa                                                          <  1:1021423/43‑1 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATa                                                          <  1:3836005/43‑1 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAgg                                                      <  1:4159643/47‑1 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                               >  1:5103044/1‑70 (MQ=255)

                                          tGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAAT‑AAGGCGAATGCCAGCACCATACCGATGGa                             <  1:1114507/71‑1 (MQ=255)

                                           accagcgccagcgTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAATCGa                             >  1:4741748/2‑71 (MQ=255)

                                            ccagcgccagcgTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAATCGa                             >  1:118170/1‑70 (MQ=255)

                                            ccagcgccagcgTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAATCGa                             <  1:408708/70‑1 (MQ=255)

                                            ccagcgccagcgTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAATCGa                             >  1:769714/1‑70 (MQ=255)

                                            ccAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATa                                     >  1:1750329/1‑62 (MQ=255)

                                            ccAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            >  1:5545325/1‑71 (MQ=255)

                                            ccAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            <  1:3234846/71‑1 (MQ=255)

                                            ccAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            >  1:1110052/1‑71 (MQ=255)

                                             cAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            >  1:5599599/1‑70 (MQ=255)

                                             cAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            >  1:2693699/1‑70 (MQ=255)

                                             cAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGCCAATCGa                             <  1:2241158/69‑1 (MQ=255)

                                                cgcCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                            >  1:4506506/1‑67 (MQ=255)

                                                 gcCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGa                        >  1:4009300/1‑70 (MQ=255)

                                                  ccAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTa                          <  1:5275838/67‑1 (MQ=255)

                                                    aGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCa                                       >  1:3025819/1‑52 (MQ=255)

                                                    aGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCa                                       >  1:1596421/1‑52 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:1405628/1‑40 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:5662538/1‑40 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:2072745/1‑40 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGa                                                   <  1:2713171/39‑1 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATggg                    >  1:248751/1‑70 (MQ=255)

                                                     gtgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATggg                    >  1:4043806/1‑70 (MQ=255)

                                                      tgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCa                  <  1:1048842/71‑1 (MQ=255)

                                                      tgtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCa                  <  1:5484630/71‑1 (MQ=255)

                                                       gtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCgtt                                <  1:33239/56‑3 (MQ=255)

                                                        ttGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                              >  1:4503950/1‑57 (MQ=255)

                                                         tgatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCaga                <  1:2691275/70‑1 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:4021223/1‑35 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:1249381/1‑35 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGaa                                                  >  1:4942369/1‑35 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCagag               <  1:892243/70‑1 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGc              >  1:5530277/1‑71 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGc              >  1:741882/1‑71 (MQ=255)

                                                          gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGc              >  1:3373119/1‑71 (MQ=255)

                                                          gatgaCAGT‑CGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGc              >  1:3243766/1‑70 (MQ=255)

                                                             gaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCAt                                      >  1:3298977/1‑44 (MQ=255)

                                                               cAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCag                 >  1:1691005/1‑63 (MQ=255)

                                                                 gTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGa                        >  1:1273389/1‑54 (MQ=255)

                                                                 gTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGa                        >  1:1311283/1‑54 (MQ=255)

                                                                 gTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCACACCGATGGAGTAGa                        >  1:5237428/1‑54 (MQ=255)

                                                                    cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCttt     >  1:1495555/1‑70 (MQ=255)

                                                                    cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCttt     >  1:5252222/1‑70 (MQ=255)

                                                                    cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCttt     <  1:2064984/70‑1 (MQ=255)

                                                                    cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCttt     >  1:4434644/1‑70 (MQ=255)

                                                                    cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCttt     >  1:3395007/1‑70 (MQ=255)

GAACGGCGAGAAGAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAG  >  Exported/2562733‑2562873

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