BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,160,663	C→T	D168N (GAT→AAT)	maf ←	Maf‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,160,663	0	C	T	98.4%	776.1 / 0.3	192	D168N (GAT→AAT)	maf	Maf‑like protein
Reads supporting (aligned to +/- strand): ref base C (3/0); new base T (87/102); total (90/102)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 194 total aligned reads displayed.

AGCGGGTTTTT‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGG  >  REL606/1160590‑1160735

aGCGGTTTTTT‑CCTTGCCCGACGTAACATCTGACACAGCGCGATTAGCTGTAAACCAACCAGCGTGTTAGGATTAc                                                                        <  2:1194259/76‑1 (MQ=255)

aGCGGGTTTTT‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTAc                                                                        >  2:5293755/1‑76 (MQ=255)

 aggggTTTTTCCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTAc                                                                        <  1:537764/74‑1 (MQ=255)

  cGGGTTTTT‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGc                                                                      >  1:2111751/1‑76 (MQ=255)

     gTTTTT‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGcccc                                                                   >  2:5195540/1‑76 (MQ=255)

      ttttt‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCct                                                                  >  1:4498443/1‑76 (MQ=255)

       tttc‑CCCTCCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctc                                                                 <  2:1816906/76‑1 (MQ=255)

       tttc‑CCCTCCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctc                                                                 <  1:5968685/76‑1 (MQ=255)

        ttt‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctct                                                                >  1:5791765/1‑76 (MQ=255)

        ttt‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctct                                                                >  1:2696984/1‑76 (MQ=255)

        ttt‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctct                                                                >  1:4974214/1‑76 (MQ=255)

        ttt‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCctct                                                                <  2:3489690/76‑1 (MQ=255)

         tt‑cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTa                                                               <  2:3112959/76‑1 (MQ=255)

          t‑cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTaa                                                              >  2:3956435/1‑76 (MQ=255)

            cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTTGGATTACGCCCCTCTaaa                                                             <  1:3941709/76‑1 (MQ=255)

            cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTaaa                                                             >  1:5190680/1‑76 (MQ=255)

            cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTaaa                                                             >  1:1054029/1‑76 (MQ=255)

            ccctccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCACTCTaaa                                                             <  1:5528321/76‑1 (MQ=255)

             cctccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCCGCGTGTTAGGATTACGCCCCTCTAAAc                                                            <  1:4376574/76‑1 (MQ=255)

             cctccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAAc                                                            <  1:2587594/76‑1 (MQ=255)

             cctccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAAc                                                            <  1:3931963/76‑1 (MQ=255)

              ctccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACg                                                           <  1:5873100/73‑1 (MQ=255)

                ccccGACGTACCATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCt                                                         <  2:3712348/75‑1 (MQ=255)

                 cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTc                                                        >  2:274373/1‑76 (MQ=255)

                 cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTc                                                        >  2:6183146/1‑76 (MQ=255)

                 cccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTc                                                        >  1:1835358/1‑76 (MQ=255)

                  ccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCa                                                       <  1:770981/76‑1 (MQ=255)

                  ccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCa                                                       >  1:2946910/1‑76 (MQ=255)

                   cGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCaa                                                      <  1:2076378/76‑1 (MQ=255)

                     aCGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAAc                                                    <  1:3274508/76‑1 (MQ=255)

                     aCGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAAc                                                    >  2:1456782/1‑76 (MQ=255)

                     aCGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAAc                                                    <  2:1084230/76‑1 (MQ=255)

                      cGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACa                                                   >  1:3000897/1‑76 (MQ=255)

                      cGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACa                                                   >  1:1521837/1‑76 (MQ=255)

                       gTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAg                                                  >  1:2134362/1‑76 (MQ=255)

                       gTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAg                                                  <  2:4616745/76‑1 (MQ=255)

                        tAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGc                                                 >  2:6060338/1‑76 (MQ=255)

                        tAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGc                                                 <  1:1415938/76‑1 (MQ=255)

                         aaCATCTGGCACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCg                                                <  1:932134/76‑1 (MQ=255)

                         aaCATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCg                                                >  2:2053731/1‑76 (MQ=255)

                          aCATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGt                                               >  1:3806118/1‑76 (MQ=255)

                           cATCTTACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTa                                              <  2:1228893/76‑1 (MQ=255)

                            ttCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTaa                                             <  1:4469531/75‑1 (MQ=255)

                            aTCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTaa                                             <  1:5050936/76‑1 (MQ=255)

                              cTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATg                                           >  1:946793/1‑76 (MQ=255)

                              cTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATg                                           <  2:4839110/76‑1 (MQ=255)

                              cTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATg                                           <  1:6241663/76‑1 (MQ=255)

                                gACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGcc                                         <  2:4879144/76‑1 (MQ=255)

                                 acacAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCa                                        >  2:2923380/1‑76 (MQ=255)

                                  cacaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCaa                                       <  2:4256112/76‑1 (MQ=255)

                                  cacaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCaa                                       <  2:5850858/76‑1 (MQ=255)

                                    caGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAAt                                     <  1:4362318/76‑1 (MQ=255)

                                    caGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAAt                                     <  2:5585445/76‑1 (MQ=255)

                                     aGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATc                                    >  2:1648873/1‑76 (MQ=255)

                                       cgcgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATctg                                  >  1:5075728/1‑74 (MQ=255)

                                       cgcgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCt                                  <  1:700052/76‑1 (MQ=255)

                                       cgcgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCt                                  >  1:4127427/1‑76 (MQ=255)

                                        gcgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCtt                                 <  2:4426084/76‑1 (MQ=255)

                                        gcgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCtt                                 <  2:2231064/76‑1 (MQ=255)

                                         cgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc                                <  1:2243435/76‑1 (MQ=255)

                                         cgATTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTc                                <  2:4320143/76‑1 (MQ=255)

                                           aTTAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAc                              >  1:1265622/1‑76 (MQ=255)

                                            ttAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAct                             <  2:3248352/76‑1 (MQ=255)

                                            ttAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAct                             >  1:218089/1‑76 (MQ=255)

                                            ttAGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCATcgga                             <  2:347847/76‑4 (MQ=255)

                                              aGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                           >  1:1090127/1‑76 (MQ=255)

                                              aGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                           >  1:864986/1‑76 (MQ=255)

                                              aGCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                           >  1:2329068/1‑76 (MQ=255)

                                               gCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCtt                          >  1:6304126/1‑76 (MQ=255)

                                               gCGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCtt                          <  2:435142/76‑1 (MQ=255)

                                                cGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTa                         >  1:16048/1‑76 (MQ=255)

                                                cGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTa                         <  1:5144122/76‑1 (MQ=255)

                                                cGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTa                         <  2:4396078/76‑1 (MQ=255)

                                                cGGTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTa                         <  1:264963/76‑1 (MQ=255)

                                                 ggTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTaa                        >  2:5378876/1‑76 (MQ=255)

                                                 ggTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTaa                        <  1:6081927/76‑1 (MQ=255)

                                                  gTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTaaa                       <  1:1646885/76‑1 (MQ=255)

                                                  gTAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTaaa                       <  1:1859688/76‑1 (MQ=255)

                                                   tAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAg                      <  2:5268755/76‑1 (MQ=255)

                                                   tAAACCAACCAGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAg                      >  1:1277651/1‑76 (MQ=255)

                                                   tAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAg                      >  2:2712178/1‑76 (MQ=255)

                                                     aaccaaccAGCGCGACAGCACCACGCCCCACTAAACGCTCAAACAACCTAATGCCAAACCCTTCACTCTTAAAGCt                    <  1:1956690/76‑1 (MQ=255)

                                                      accaaccaGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTa                   >  1:4838420/1‑76 (MQ=255)

                                                        caaccaGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTAcc                 <  2:5047500/76‑1 (MQ=255)

                                                          accaGCGTGTTAGGATTACGCCCCTCTAAACGCGCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCcg               <  2:704860/76‑1 (MQ=255)

                                                           ccaGCGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCcgc              >  2:1908457/1‑76 (MQ=255)

                                                               cGTGTTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAg          <  2:4189346/76‑1 (MQ=255)

                                                                 tgtTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGta        <  2:468045/76‑2 (MQ=255)

                                                                  gtTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGc       <  1:3546495/76‑1 (MQ=255)

                                                                  gtTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGc       >  2:4566585/1‑76 (MQ=255)

                                                                  gtTAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGc       <  1:1034313/76‑1 (MQ=255)

                                                                   ttAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCa      <  2:2307434/76‑1 (MQ=255)

                                                                   ttAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCa      <  1:2519073/76‑1 (MQ=255)

                                                                    tAGGATTTCGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAg     >  2:845175/1‑76 (MQ=255)

                                                                    tAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAg     >  2:3845024/1‑76 (MQ=255)

                                                                    tAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAg     >  2:341383/1‑76 (MQ=255)

                                                                    tAGGATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAg     >  1:682614/1‑76 (MQ=255)

                                                                      ggATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAggg   <  2:5171615/76‑1 (MQ=255)

                                                                      ggATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAggg   <  1:5321869/76‑1 (MQ=255)

                                                                       gATTACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAgggg  >  1:2005411/1‑76 (MQ=255)

AGCGGGTTTTT‑CCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGG  >  REL606/1160590‑1160735

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