BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,296	C→T	L198L (CTG→CTA)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,296	0	C	T	75.6%	160.6 / 61.3	90	L198L (CTG→CTA)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base C (11/11); new base T (30/38); total (41/49)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

Only 100 of 102 total aligned reads displayed.

AAATCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGGTG  >  NC_000913/1715249‑1715343

aaaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtag                                               <  1:3933356/50‑1 (MQ=255)

aaaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtag                                               <  1:2886594/50‑1 (MQ=255)

aaaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtag                                               <  1:1940648/50‑1 (MQ=255)

 aaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtagc                                              <  1:6136306/50‑1 (MQ=255)

 aaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtagc                                              <  1:6077853/50‑1 (MQ=255)

 aaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtagc                                              <  1:1391719/50‑1 (MQ=255)

 aaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtag                                               >  1:5450645/1‑49 (MQ=255)

 aaTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACcagc                                              <  1:5952338/50‑1 (MQ=255)

  aTCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtagca                                             <  1:7360054/50‑1 (MQ=255)

   tCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTTACtagcag                                            >  1:656273/1‑50 (MQ=255)

   tCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACtagcag                                            >  1:3242760/1‑50 (MQ=255)

     caccaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGGGCTgcagca                                          >  1:5409187/1‑50 (MQ=38)

     caccaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTagcagc                                           <  1:2217281/49‑1 (MQ=255)

     caccaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAgcagca                                          >  1:5141341/1‑50 (MQ=255)

     caccaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAgcagca                                          <  1:3557052/50‑1 (MQ=255)

     caccaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAgcagca                                          <  1:2184354/50‑1 (MQ=255)

      accaGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGAACAGCAGCAt                                         >  1:1516915/1‑50 (MQ=255)

        caGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCAttt                                       <  1:4358428/50‑1 (MQ=255)

        caGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCAttt                                       <  1:961964/50‑1 (MQ=255)

        caGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCAttt                                       >  1:2027498/1‑50 (MQ=255)

         aGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCAttt                                       >  1:5317019/1‑49 (MQ=255)

         aGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTc                                      <  1:93628/50‑1 (MQ=255)

         aGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTc                                      <  1:6872969/50‑1 (MQ=255)

         aGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTc                                      <  1:3184223/50‑1 (MQ=255)

         aGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCActt                                       >  1:6585766/1‑49 (MQ=255)

          gCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTc                                      >  1:7715563/1‑49 (MQ=255)

           cGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCaa                                    <  1:4132804/50‑1 (MQ=255)

           cGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCaa                                    <  1:780771/50‑1 (MQ=255)

           cGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCa                                     >  1:2003993/1‑49 (MQ=255)

            ggACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCaa                                    >  1:2463479/1‑49 (MQ=255)

            ggACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCaa                                    >  1:6617286/1‑49 (MQ=255)

             gACGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCCTTTCaaaa                                  >  1:3061216/1‑50 (MQ=255)

             gACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCaaaa                                  <  1:2977850/50‑1 (MQ=255)

             gACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCaaaa                                  >  1:5746745/1‑50 (MQ=255)

              aCGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAAc                                 <  1:6421950/50‑1 (MQ=255)

              aCGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAAc                                 >  1:2797248/1‑50 (MQ=255)

              aCGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAAc                                 <  1:344596/50‑1 (MQ=255)

               cGGCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACg                                <  1:3385884/50‑1 (MQ=255)

                ggCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACgg                               <  1:3436271/50‑1 (MQ=255)

                ggCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACgg                               >  1:5819979/1‑50 (MQ=255)

                 gCTGATATGCCAGGCTTCATCGGCGGTGTTTTGCAGCATTTCAAAACGGt                              >  1:1013337/1‑50 (MQ=37)

                 gCTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGt                              <  1:8185234/50‑1 (MQ=255)

                 gCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGt                              >  1:2932249/1‑50 (MQ=255)

                  cTGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTc                             <  1:4463294/50‑1 (MQ=255)

                   tGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCa                            >  1:1724387/1‑50 (MQ=255)

                   tGATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAAAGGTCa                            >  1:5688250/1‑50 (MQ=255)

                    gATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCAc                           >  1:1473170/1‑50 (MQ=255)

                    gATATGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACCGTccc                           >  1:4101521/1‑48 (MQ=255)

                     atatGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCa                            <  1:6195209/48‑1 (MQ=255)

                     atatGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCAcc                          >  1:3270925/1‑49 (MQ=255)

                     atatGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCAc                           <  1:8066689/49‑1 (MQ=255)

                     atatGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCAc                           <  1:1792138/49‑1 (MQ=255)

                     atatGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCAc                           <  1:4354751/49‑1 (MQ=255)

                      tatGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACgg                         <  1:2479352/50‑1 (MQ=255)

                       atGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACgg                         >  1:5247839/1‑49 (MQ=255)

                       atGCCAGGCTTCATCGGCGGTGACTAGCAACATTTCAAAACGGCCACCGc                        >  1:845517/1‑50 (MQ=255)

                        tGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCt                       >  1:8317311/1‑50 (MQ=255)

                        tGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCt                       <  1:6313444/50‑1 (MQ=255)

                        tGCCAGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCt                       <  1:1748895/50‑1 (MQ=255)

                        tGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCACGGCt                       <  1:5224605/50‑1 (MQ=255)

                           cAGGCTTCATCGGCGGTGACTAGCAGCATTTTAAAACGGTCACGGCTg                      >  1:907254/1‑48 (MQ=39)

                           cAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCACGGCTGTa                    <  1:2950652/50‑1 (MQ=255)

                            aGGCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAg                   <  1:2691257/50‑1 (MQ=255)

                             ggCTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGc                  <  1:1179028/50‑1 (MQ=255)

                               cTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCa                <  1:5314874/50‑1 (MQ=255)

                               cTTCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCa                <  1:3893686/50‑1 (MQ=255)

                                ttCATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAg               >  1:4886463/1‑50 (MQ=255)

                                 tCATCGGCGGTGACTAGCAGCATTTTAAAACGGTCACGGCTGTAGCCAGc              >  1:2162359/1‑50 (MQ=255)

                                  cATCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAg               <  1:3986611/48‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCt             <  1:7183890/49‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCt             <  1:3851213/49‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCt             <  1:6677998/49‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTc            <  1:2326755/50‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTc            <  1:3705158/50‑1 (MQ=255)

                                   aTCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTc            <  1:5156305/50‑1 (MQ=255)

                                   aTCGGCGGTGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTc            <  1:5345080/50‑1 (MQ=255)

                                    tCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcg           <  1:2361340/50‑1 (MQ=255)

                                    tCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcg           <  1:5869302/50‑1 (MQ=255)

                                    tCGGCGGTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcg           <  1:5755705/50‑1 (MQ=255)

                                     cggcggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgc          >  1:2713198/1‑50 (MQ=255)

                                      ggcggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgcg         >  1:1820140/1‑50 (MQ=255)

                                      ggcggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgcg         <  1:4253090/50‑1 (MQ=255)

                                      ggcggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgcg         >  1:710648/1‑50 (MQ=255)

                                      ggcggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgc          <  1:4076368/49‑1 (MQ=255)

                                      ggcggTGACTAGCAGCATTTCAAAACGGTCACGACTGTAGCCAGCTcgc          <  1:1001964/49‑1 (MQ=255)

                                      ggcggTGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTcgc          <  1:1816758/49‑1 (MQ=255)

                                        cggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCgcgc        >  1:312931/1‑49 (MQ=255)

                                        cggTGACTAGCAGCATTTC‑AAACGGTCACGGCTGTAGCCAGCTCgcgc        >  1:5855634/1‑48 (MQ=38)

                                         ggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCgcgc        >  1:6213783/1‑48 (MQ=255)

                                         ggTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCa      <  1:6664187/50‑1 (MQ=255)

                                         ggTGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCa      <  1:4961547/50‑1 (MQ=255)

                                          gTGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAg     >  1:1682499/1‑50 (MQ=255)

                                           tGACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAgg    >  1:8589450/1‑50 (MQ=255)

                                           tGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAgg    <  1:6630158/50‑1 (MQ=255)

                                            gACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgt   >  1:7885995/1‑50 (MQ=255)

                                            gACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgt   >  1:5821741/1‑50 (MQ=255)

                                            gACTAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgt   >  1:875000/1‑50 (MQ=255)

                                             aCCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgtg  >  1:4399692/1‑50 (MQ=255)

                                             aCCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgtg  >  1:3305223/1‑50 (MQ=255)

                                             aCCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGgt   <  1:7319301/49‑1 (MQ=255)

AAATCCACCAGCGGACGGCTGATATGCCAGGCTTCATCGGCGGTGACCAGCAGCATTTCAAAACGGTCACGGCTGTAGCCAGCTCGCGCCAGGTG  >  NC_000913/1715249‑1715343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`