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.3%	218.6 / ‑4.0	58	S9L (TCG→TTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (26/31); total (26/32)
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.73e-01

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGG  >  REL606/3497938‑3498081

gggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATg                                                                      >  1:238886/1‑76 (MQ=255)

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

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

   tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATACAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAat                                                                   >  2:804349/1‑76 (MQ=255)

   tCCTTTACAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTATGGGCTTCTACGCTAAGCGTGGGCAATGAat                                                                   <  1:1484106/76‑1 (MQ=255)

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

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

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

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

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

         ccAAGCCTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaa                                                             >  1:86539/1‑76 (MQ=255)

          cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaaa                                                            <  2:1219233/76‑1 (MQ=255)

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

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

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

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

               gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAg                                                       >  1:319520/1‑76 (MQ=255)

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

                 aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCa                                                     >  2:795242/1‑76 (MQ=255)

                 aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCa                                                     <  1:388865/76‑1 (MQ=255)

                   aCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCAtt                                                   <  2:482665/76‑1 (MQ=255)

                   aCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCAtt                                                   <  2:800192/76‑1 (MQ=255)

                      tCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATAGGTAAACGGAGCATTCat                                                >  1:172305/1‑76 (MQ=255)

                       cAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCata                                               >  2:550397/1‑76 (MQ=255)

                         gCTTGTAAGCGATATAGTGCTTAAAGGCGTCNACGCTAAGAGTGGGNAATGANTATGTAAACGGAGCATTCAtata                                             <  2:214166/76‑1 (MQ=255)

                          cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTGATATAg                                            <  2:436417/76‑1 (MQ=255)

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

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

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

                              tAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGTAATGAATATGTAAACGGAGCATTCATATAGGCGt                                        <  2:447291/76‑1 (MQ=255)

                              tAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGt                                        >  2:714884/1‑76 (MQ=255)

                               aaGCGATAGGGTGCTTAAGGGCTNCTACGNTAAGCTTGGGCAATNAACATGAAAACGGAGCATTCATATAGGCGtt                                       >  2:210138/1‑76 (MQ=255)

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

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

                                   gATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt                                   >  1:555366/1‑76 (MQ=255)

                                   gATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt                                   >  1:1227184/1‑76 (MQ=255)

                                   gATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt                                   <  2:1294108/76‑1 (MQ=255)

                                    aTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                                  <  2:1430959/76‑1 (MQ=255)

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

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

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

                                           gCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                           >  2:1266928/1‑76 (MQ=255)

                                           gCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                           >  1:916748/1‑76 (MQ=255)

                                            cTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                          >  1:1335636/1‑76 (MQ=255)

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

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

                                                   gCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                   >  2:443882/1‑76 (MQ=255)

                                                    cTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCg                  <  1:966557/76‑1 (MQ=255)

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

                                                      tcgaCGCTAAGCGNGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                <  1:210138/73‑1 (MQ=255)

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

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

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

                                                           gCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaa           <  2:933984/76‑1 (MQ=255)

                                                           gCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaa           >  2:96450/1‑76 (MQ=255)

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

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

                                                                    gggCAATGAATATGTAAACGGAGCATTCATATANGCGTTTCCTGAAAACTATTTCANGCGATAGTCANGCTCACgg  >  2:1362738/1‑76 (MQ=255)

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGG  >  REL606/3497938‑3498081

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