BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,498,010	G→A	S9L (TCG→TTG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,498,010	0	G	A	98.9%	1093.7 / 0.2	268	S9L (TCG→TTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (141/124); total (143/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 279 total aligned reads displayed.

CGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGT  >  REL606/3497937‑3498082

 gggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATg                                                                       <  2:3620437/76‑1 (MQ=255)

  ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGa                                                                      <  1:4330787/76‑1 (MQ=255)

    tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAat                                                                    >  2:1507862/1‑76 (MQ=255)

    tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAat                                                                    <  2:6397012/76‑1 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGGAATGAatat                                                                  >  1:3139134/1‑76 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  >  1:1398072/1‑76 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  >  1:6883545/1‑76 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  >  1:3008913/1‑76 (MQ=255)

       tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATg                                                                 <  1:1461582/76‑1 (MQ=255)

       tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATg                                                                 >  1:5474335/1‑76 (MQ=255)

        ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGt                                                                <  1:6814953/76‑1 (MQ=255)

        ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGt                                                                >  2:2533265/1‑76 (MQ=255)

        ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGt                                                                >  1:4681330/1‑76 (MQ=255)

         tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTa                                                               <  2:4549639/76‑1 (MQ=255)

          ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaa                                                              <  2:4364483/76‑1 (MQ=255)

           cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGGGGGCAATGAATATGTaaa                                                             >  2:199671/1‑76 (MQ=255)

            aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAAc                                                            <  1:1830546/76‑1 (MQ=255)

            aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAAc                                                            <  2:2870400/76‑1 (MQ=255)

            aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAAc                                                            <  1:7632477/76‑1 (MQ=255)

             aTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACg                                                           >  1:6151113/1‑76 (MQ=255)

               cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATg                                                                       >  1:6367467/1‑62 (MQ=255)

               cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGa                                                         <  2:6914351/76‑1 (MQ=255)

                 tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATTTAAACGGAGc                                                       <  2:850635/76‑1 (MQ=255)

                 tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGc                                                       <  1:2644729/76‑1 (MQ=255)

                 tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGc                                                       >  2:800708/1‑76 (MQ=255)

                   aaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCAt                                                     <  2:5221351/76‑1 (MQ=255)

                     cATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTc                                                   <  2:113634/76‑1 (MQ=255)

                     cATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTc                                                   <  1:4628301/76‑1 (MQ=255)

                     cATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTc                                                   <  1:1992063/76‑1 (MQ=255)

                      aTCAGCTTGTAAGCGGTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGGCATGAATATGTAAACGGGGGATTCa                                                  >  2:6142688/1‑76 (MQ=255)

                      aTCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCa                                                  <  2:2107762/76‑1 (MQ=255)

                       tCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCat                                                 <  1:403966/76‑1 (MQ=255)

                         aGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCatat                                               <  1:740503/76‑1 (MQ=255)

                         aGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCatat                                               >  2:4314069/1‑76 (MQ=255)

                          gCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCAtata                                              >  1:6701321/1‑76 (MQ=255)

                          gCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCAtata                                              <  1:6186793/76‑1 (MQ=255)

                           cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAg                                             <  1:385016/76‑1 (MQ=255)

                            ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAgg                                            >  1:1002762/1‑76 (MQ=255)

                            ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAgg                                            >  2:19008/1‑76 (MQ=255)

                            ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAgg                                            >  1:924749/1‑76 (MQ=255)

                             tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                           >  1:7565397/1‑76 (MQ=255)

                             tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGc                                           >  1:2739640/1‑76 (MQ=255)

                              gTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCg                                          <  2:3354037/76‑1 (MQ=255)

                              gTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCg                                          >  1:5043383/1‑76 (MQ=255)

                               gaaGCGAGAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGt                                         <  2:5699244/75‑1 (MQ=255)

                                aaGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGtt                                        >  1:3607267/1‑76 (MQ=255)

                                 aGCGATAGAGTTCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGttt                                       <  1:7805592/76‑1 (MQ=255)

                                 aGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGttt                                       <  1:5979515/76‑1 (MQ=255)

                                 aGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGttt                                       >  2:7233622/1‑76 (MQ=255)

                                   cGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                     >  1:2296311/1‑76 (MQ=255)

                                   cGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                     >  2:2059079/1‑76 (MQ=255)

                                   cGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                     >  1:4201992/1‑76 (MQ=255)

                                   cGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                     >  1:5296394/1‑76 (MQ=255)

                                     aTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                                   <  1:1894455/76‑1 (MQ=255)

                                     aTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                                   <  2:4725857/76‑1 (MQ=255)

                                     aTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                                   >  1:2461440/1‑76 (MQ=255)

                                      tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                  >  2:6732773/1‑76 (MQ=255)

                                      tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                  >  2:238551/1‑76 (MQ=255)

                                       agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                 >  1:1749598/1‑76 (MQ=255)

                                        gagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaa                                <  2:6282138/76‑1 (MQ=255)

                                         agTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa                               <  2:3171625/76‑1 (MQ=255)

                                         agTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa                               >  2:7926447/1‑76 (MQ=255)

                                          gTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                              >  1:6925183/1‑76 (MQ=255)

                                          gTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                              <  1:2062164/76‑1 (MQ=255)

                                           tGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt                             >  2:1946751/1‑76 (MQ=255)

                                            gCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                            <  1:1680501/76‑1 (MQ=255)

                                             cTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                           <  2:5580218/76‑1 (MQ=255)

                                              ttAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                          >  1:4137514/1‑76 (MQ=255)

                                               tAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                         <  2:1250420/76‑1 (MQ=255)

                                               tAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                         >  1:3213419/1‑76 (MQ=255)

                                               tAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                         >  1:6258374/1‑76 (MQ=255)

                                               tAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                         >  1:4808466/1‑76 (MQ=255)

                                                aaGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                        >  1:5209464/1‑76 (MQ=255)

                                                 aGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                       <  1:2154257/76‑1 (MQ=255)

                                                 aGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                       >  1:5518776/1‑76 (MQ=255)

                                                  gggCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                      >  2:7647108/1‑76 (MQ=255)

                                                  gggCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                      <  2:4582791/76‑1 (MQ=255)

                                                  gggCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                      >  1:4782251/1‑76 (MQ=255)

                                                  gggCTTCTACGCTAAGCGGGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                      >  1:6972134/1‑76 (MQ=255)

                                                    gCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                    <  1:7371521/76‑1 (MQ=255)

                                                      ttCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                  >  2:6237544/1‑76 (MQ=255)

                                                      ttCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                  >  2:295432/1‑76 (MQ=255)

                                                        cTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa                >  2:1385109/1‑76 (MQ=255)

                                                          aCGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt              <  1:1189540/76‑1 (MQ=255)

                                                           cGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTa             >  1:2900149/1‑76 (MQ=255)

                                                             cTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa           >  2:5830059/1‑76 (MQ=255)

                                                             cTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa           <  2:7510309/76‑1 (MQ=255)

                                                             cTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa           >  1:1067486/1‑76 (MQ=255)

                                                             cTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa           >  1:1627171/1‑76 (MQ=255)

                                                               aaGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCACGCGATAGTAAAGc         >  1:7065021/1‑76 (MQ=255)

                                                               aaGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGc         >  1:3725765/1‑76 (MQ=255)

                                                                  cGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa      >  2:6028016/1‑76 (MQ=255)

                                                                   gTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCAc     >  2:1787473/1‑76 (MQ=255)

                                                                   gTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCAc     <  2:7363847/76‑1 (MQ=255)

                                                                    tGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg    >  2:5134538/1‑76 (MQ=255)

                                                                    tGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg    >  1:831965/1‑76 (MQ=255)

                                                                    tGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg    >  1:2431585/1‑76 (MQ=255)

                                                                      ggCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGt  <  1:7284478/76‑1 (MQ=255)

                                                                      ggCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGt  <  1:3795168/76‑1 (MQ=255)

                                                                      ggCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGt  <  1:3462739/76‑1 (MQ=255)

CGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGT  >  REL606/3497937‑3498082

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