BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,097,424	A→T	L50Q (CTG→CAG)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,097,424	0	A	T	99.5%	756.4 / ‑3.6	182	L50Q (CTG→CAG)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (90/91); minor base G (0/1); total (90/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 183 total aligned reads displayed.

TGCAATGGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCAGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAACCCA  >  REL606/4097351‑4097496

 gcaatgGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTc                                                                       >  1:737382/1‑76 (MQ=255)

 gcaatgGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTc                                                                       >  1:875563/1‑76 (MQ=255)

   aatgGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTccc                                                                     <  2:1470968/76‑1 (MQ=255)

    atgGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTcccc                                                                    >  1:5105401/1‑76 (MQ=255)

     tgGTAACAGCGGGATTAAGATGCGCGCCGGAAGCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCa                                                                   <  2:6521080/76‑1 (MQ=255)

     tgGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCa                                                                   <  1:2531747/76‑1 (MQ=255)

       gTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCaaa                                                                 >  2:4058213/1‑76 (MQ=255)

         aaCAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATg                                                               <  2:2412897/76‑1 (MQ=255)

           cAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGac                                                             <  1:3610900/76‑1 (MQ=255)

            aGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGaca                                                            >  2:1447817/1‑76 (MQ=255)

            aGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCGGTCCCCAAATGaca                                                            <  2:4356612/76‑1 (MQ=255)

              cGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACt                                                          >  2:2212152/1‑76 (MQ=255)

                 gATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGAt                                                       <  1:5542816/76‑1 (MQ=255)

                  aTTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGAtt                                                      <  1:4094861/76‑1 (MQ=255)

                  aTTAAGATGCGCGCCGGAAACCCCCGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGAtt                                                      >  2:588953/1‑76 (MQ=255)

                   ttAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGAttt                                                     >  2:1566475/1‑76 (MQ=255)

                   ttAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGAttt                                                     >  1:1780476/1‑76 (MQ=255)

                     aaGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTcc                                                   >  1:1260821/1‑76 (MQ=255)

                     aaGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTcc                                                   >  1:2634844/1‑76 (MQ=255)

                     aaGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTcc                                                   <  2:2515431/76‑1 (MQ=255)

                     aaGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTcc                                                   >  1:5399795/1‑76 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  <  1:7039490/76‑1 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  >  1:7011754/1‑76 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  <  1:6595042/76‑1 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  <  1:7138483/76‑1 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  <  1:3774791/76‑1 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  <  2:1073415/76‑1 (MQ=255)

                      aGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTccc                                                  >  1:3393769/1‑76 (MQ=255)

                        aTGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCAc                                                <  2:3741129/76‑1 (MQ=255)

                         tGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACt                                               >  2:3493755/1‑76 (MQ=255)

                         tGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACt                                               >  2:799821/1‑76 (MQ=255)

                         tGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACt                                               >  1:3820956/1‑76 (MQ=255)

                          gcgcgcCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTg                                              <  1:4357681/76‑1 (MQ=255)

                           cgcgcCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGa                                             >  1:2268191/1‑76 (MQ=255)

                           cgcgcCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGa                                             >  2:6810658/1‑76 (MQ=255)

                             cgcCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGAcc                                           <  1:3900350/76‑1 (MQ=255)

                             cgcCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGAcc                                           <  2:2812124/76‑1 (MQ=255)

                               ccGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCaa                                         >  2:2062921/1‑76 (MQ=255)

                               ccGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCaa                                         >  2:4949940/1‑76 (MQ=255)

                               ccGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCaa                                         <  1:3230879/76‑1 (MQ=255)

                                 ggAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCaaaa                                       <  1:4882045/76‑1 (MQ=255)

                                 ggAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCaaaa                                       >  1:4161864/1‑76 (MQ=255)

                                  gAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAg                                      >  2:2990978/1‑76 (MQ=255)

                                  gAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAg                                      <  1:4810272/76‑1 (MQ=255)

                                  gAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAg                                      <  1:3594321/76‑1 (MQ=255)

                                   aaaCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGa                                     >  2:4167279/1‑76 (MQ=255)

                                   aaaCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGa                                     >  2:6650876/1‑76 (MQ=255)

                                    aaCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGAc                                    <  1:3508180/76‑1 (MQ=255)

                                    aaCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGAc                                    >  2:6544183/1‑76 (MQ=255)

                                     aCCCCTGCGGTCAGGTAGATGGCCATTGCTACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACg                                   <  2:6606086/76‑1 (MQ=255)

                                     aCCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACg                                   <  2:1540219/76‑1 (MQ=255)

                                      ccccTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGc                                  >  2:566947/1‑76 (MQ=255)

                                       cccTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCa                                 <  1:279154/76‑1 (MQ=255)

                                        ccTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCAc                                >  2:3916202/1‑76 (MQ=255)

                                        ccTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCAc                                <  2:2233200/76‑1 (MQ=255)

                                        ccTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACAAAAAGACGCAc                                <  2:140271/76‑1 (MQ=255)

                                         cTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCAcc                               <  2:1769707/76‑1 (MQ=255)

                                         cTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCAcc                               <  2:1141145/76‑1 (MQ=255)

                                         cTGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCAcc                               <  2:6733441/76‑1 (MQ=255)

                                          tGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCa                              >  2:4886356/1‑76 (MQ=255)

                                          tGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCa                              >  2:5538057/1‑76 (MQ=255)

                                          tGCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCa                              <  2:6033896/76‑1 (MQ=255)

                                           gCGGTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAg                             >  2:5470690/1‑76 (MQ=255)

                                             ggTCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCg                           <  2:4714476/76‑1 (MQ=255)

                                               tCAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGAc                         >  1:2473498/1‑76 (MQ=255)

                                               gcAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGAc                         <  2:5248982/75‑1 (MQ=255)

                                                cAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACt                        >  2:3226338/1‑76 (MQ=255)

                                                cAGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACt                        <  1:6053077/76‑1 (MQ=255)

                                                 aGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtt                       >  2:3334083/1‑76 (MQ=255)

                                                 aGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtt                       <  2:6293962/76‑1 (MQ=255)

                                                 aGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtt                       >  1:3446327/1‑76 (MQ=255)

                                                 aGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtt                       >  1:4508077/1‑76 (MQ=255)

                                                 aGGTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtt                       >  2:1364439/1‑76 (MQ=255)

                                                   gTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtttt                     >  1:5275330/1‑76 (MQ=255)

                                                   gTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtttt                     <  1:6136013/76‑1 (MQ=255)

                                                   gTAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACtttt                     <  2:1130077/76‑1 (MQ=255)

                                                    tAGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTa                    >  1:1949145/1‑76 (MQ=255)

                                                     aGATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAg                   <  2:5297929/76‑1 (MQ=255)

                                                      gATGGCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGt                  >  2:4304920/1‑76 (MQ=255)

                                                         ggCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGtgct               >  2:2142678/1‑76 (MQ=255)

                                                         ggCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGtgct               >  2:470551/1‑76 (MQ=255)

                                                         ggCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGtgct               >  1:3623320/1‑76 (MQ=255)

                                                          gCCATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGtgctg              <  2:4549909/76‑1 (MQ=255)

                                                           ccATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGtgctgc             >  2:4630784/1‑76 (MQ=255)

                                                            cATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgca            <  1:5740510/76‑1 (MQ=255)

                                                            cATTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgca            <  1:2878527/76‑1 (MQ=255)

                                                             aTTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaa           >  1:5644064/1‑76 (MQ=255)

                                                             aTTGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaa           >  1:5361230/1‑76 (MQ=255)

                                                              ttGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaac          >  2:6204037/1‑76 (MQ=255)

                                                              ttGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaac          <  1:2695360/76‑1 (MQ=255)

                                                               tGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacg         >  1:3054052/1‑76 (MQ=255)

                                                               tGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacg         <  2:5399795/76‑1 (MQ=255)

                                                               tGCCACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacg         <  1:5955546/76‑1 (MQ=255)

                                                                 ccACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacgca       <  1:2623111/76‑1 (MQ=255)

                                                                  cACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacgcaa      <  2:6251747/76‑1 (MQ=255)

                                                                  cACCCCCTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTgcaacgcaa      >  2:6438250/1‑76 (MQ=255)

                                                                     ccccTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAAccc   <  1:6322464/76‑1 (MQ=255)

                                                                     ccccTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAAccc   <  1:1277583/76‑1 (MQ=255)

                                                                      cccTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAACCca  >  2:2918905/1‑76 (MQ=255)

                                                                      cccTGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAACCca  <  2:6949168/76‑1 (MQ=255)

TGCAATGGTAACAGCGGGATTAAGATGCGCGCCGGAAACCCCTGCGGTCAGGTAGATGGCCATTGCCACCCCCAGTCCCCAAATGACACTGATTTCCCACTGACCAAAAGACGCACCAGCGACTTTTAGTGCTGCAACGCAACCCA  >  REL606/4097351‑4097496

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