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	99.1%	888.6 / ‑2.0	221	S9L (TCG→TTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (105/114); total (106/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 233 total aligned reads displayed.

CGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGT  >  REL606/3497937‑3498082

cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAAt                                                                        <  1:3875118/76‑1 (MQ=255)

cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAAt                                                                        >  1:1660124/1‑76 (MQ=255)

cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAAt                                                                        >  1:2479799/1‑76 (MQ=255)

   gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGaa                                                                     <  1:5647522/76‑1 (MQ=255)

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

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

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

    tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTAAGCTAAGCGTGGGCAAAGAat                                                                    >  2:790780/1‑76 (MQ=255)

     ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAata                                                                   <  2:3882664/76‑1 (MQ=255)

     ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAata                                                                   >  1:367696/1‑76 (MQ=255)

     ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAata                                                                   <  1:6055241/76‑1 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGGTAGAGTGCTTAAGGGGTTCTACGCTAAGCGTGGGCCATGAatat                                                                  >  2:5456298/1‑76 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  <  2:1166158/76‑1 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  <  2:1861341/76‑1 (MQ=255)

      cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAatat                                                                  <  1:3226564/76‑1 (MQ=255)

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

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

         tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTa                                                               >  2:5653287/1‑76 (MQ=255)

          ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaa                                                              >  2:2940539/1‑76 (MQ=255)

          ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaa                                                              >  1:3600936/1‑76 (MQ=255)

            aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAAc                                                            >  2:3917060/1‑76 (MQ=255)

              tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACgg                                                          <  2:275795/76‑1 (MQ=255)

              tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACgg                                                          <  2:1793750/76‑1 (MQ=255)

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

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

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

                    aCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCAtt                                                    >  2:6075025/1‑76 (MQ=255)

                      aTCAGCTTGTAAGCGATAGAGTGCTTAAGGGGTTCTACGCTAAGCGTGGGGAATGAATATGTAAACGGAGGAATCa                                                  >  2:3512351/1‑76 (MQ=255)

                      aTCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAACGTGGGCAATGAATATGTAAACGGAGCGTTCa                                                  >  2:4230882/1‑76 (MQ=255)

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

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

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

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

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

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

                           cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAg                                             >  2:4476549/1‑76 (MQ=255)

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

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

                             tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAgcc                                           >  2:887224/1‑74 (MQ=255)

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

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

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

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

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

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

                                aaGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGtt                                        <  1:3558544/76‑1 (MQ=255)

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

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

                                   cGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                     <  1:5770431/76‑1 (MQ=255)

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

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

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

                                       agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                 <  1:1617417/76‑1 (MQ=255)

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

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

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

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

                                            gCTTAAGGGCTTGTACGCTAACCGTGGGCAATGAATATGTAAAGGGCGCCTTCATATAGCCGATTCCTGAAAACTa                            >  2:1119143/1‑76 (MQ=255)

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

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

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

                                                aaGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                        <  1:2945392/76‑1 (MQ=255)

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

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

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

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

                                                aaGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                        <  1:1226911/76‑1 (MQ=255)

                                                 aGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                       <  1:1670870/76‑1 (MQ=255)

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

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

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

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

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

                                                   ggCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAg                     >  1:724798/1‑76 (MQ=255)

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

                                                       tCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                 <  1:213354/76‑1 (MQ=255)

                                                       tCTACGCTAAGCGGGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                 >  2:6296760/1‑76 (MQ=255)

                                                        cTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa                <  2:2054452/76‑1 (MQ=255)

                                                         tACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                 <  2:1909378/74‑1 (MQ=255)

                                                         tACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                 >  1:1909378/1‑74 (MQ=255)

                                                         tACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAg               <  1:697070/76‑1 (MQ=255)

                                                          aCGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt              >  1:101516/1‑76 (MQ=255)

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

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

                                                              tAAGGGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg          <  2:4130333/76‑1 (MQ=255)

                                                              tAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg          >  1:5207730/1‑76 (MQ=255)

                                                              tAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg          >  1:2682018/1‑76 (MQ=255)

                                                              tAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg          >  1:4735479/1‑76 (MQ=255)

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

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

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

                                                                aGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt        >  1:3975389/1‑76 (MQ=255)

                                                                 gCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTc       <  1:5152042/76‑1 (MQ=255)

                                                                  cGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa      <  1:4595199/76‑1 (MQ=255)

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

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

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

                                                                    tGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg    <  2:5911343/76‑1 (MQ=255)

                                                                     gggCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg   >  2:3324309/1‑76 (MQ=255)

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

CGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGT  >  REL606/3497937‑3498082

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