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.5%	264.4 / NA	67	S9L (TCG→TTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (32/34); total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGG  >  REL606/3497938‑3498081

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

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

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

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

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

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

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

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

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

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

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

          cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaaa                                                            >  2:186023/1‑76 (MQ=255)

          cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaaa                                                            >  2:809337/1‑76 (MQ=255)

          cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTaaa                                                            >  2:487665/1‑76 (MQ=255)

            aTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACg                                                          <  2:252569/76‑1 (MQ=255)

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

              cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGa                                                        >  1:1634482/1‑76 (MQ=255)

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

              cGTAAACATCAGCCTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGa                                                        >  1:798438/1‑76 (MQ=255)

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

                  aaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCAt                                                    >  2:974190/1‑76 (MQ=255)

                      tCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCat                                                >  1:769569/1‑76 (MQ=255)

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

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

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

                           ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAgn                                           <  1:1213507/76‑2 (MQ=255)

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

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

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

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

                                 gCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                     >  2:732820/1‑76 (MQ=255)

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

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:145428/76‑1 (MQ=255)

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:1296330/76‑1 (MQ=255)

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:1587229/76‑1 (MQ=255)

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:297864/76‑1 (MQ=255)

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:169322/76‑1 (MQ=255)

                                     tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                 <  2:920827/76‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                     ttCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                 <  2:991914/76‑1 (MQ=255)

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

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

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

                                                         aCGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt             <  1:954041/76‑1 (MQ=255)

                                                          cGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTa            <  1:889212/76‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                             tAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg         <  1:1103025/76‑1 (MQ=255)

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

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

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

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

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

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

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGG  >  REL606/3497938‑3498081

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