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.0%	179.1 / NA	51	S9L (TCG→TTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (24/26); total (24/27)
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:857970/1‑76 (MQ=255)

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

  gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGaa                                                                    >  1:863971/1‑76 (MQ=255)

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

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

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

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

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

        tCCAATCGTAAACATCAGCATGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTa                                                              >  2:519153/1‑76 (MQ=255)

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

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

             tCGTAAACATCAGCTTGTAAGCGATAGAGAGCTTAAGGGCTTCGACGCTAANCGTGAGCGATGAATATGTAAGCgg                                                         <  1:160447/76‑1 (MQ=255)

                tAAACATCAGCTTGTAAGCGTTAGAGTGCTTAAGGGCTTCTACGCAAAGCGTGGGCAATCAATATGTATACCGAGc                                                      >  1:33177/1‑76 (MQ=255)

                tAAACATCAGCTTGTAAGCGATAGAGTNCTTAAGGGCTTATACGCTAAGCGTGGGCAATGAATATGTAAACGGAGc                                                      >  1:1612672/1‑76 (MQ=255)

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

                 aaaCATCAGCTTGTAAGCGAGAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGACAATGAATATGTAAACGGAGCa                                                     >  2:148831/1‑76 (MQ=255)

                     aTCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCa                                                 >  1:835428/1‑76 (MQ=255)

                     aTCAGCTTGTAAGCCATCGCGTTCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGCATATGTAAACGGAGCATTCa                                                 <  1:84818/76‑1 (MQ=255)

                       cAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGANCATTCata                                               <  1:1756168/76‑1 (MQ=255)

                       cAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCata                                               <  1:1319927/76‑1 (MQ=255)

                          cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAATCGTGGGCAATGAATATGTAAACGGAGCATTCATATAg                                            <  2:1043368/76‑1 (MQ=255)

                           ttGTAAGCGATANAGNGCTTAAGGGCTTCTACGCNAAGCGTCGGCAATGAATATGTAAACGGAGCATTCATATAgg                                           <  1:1496084/76‑1 (MQ=255)

                           ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGCGCAATGAATATGTAAACGGAGCATTCATATAgg                                           <  1:1678194/76‑1 (MQ=255)

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

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

                                 gCGATAGAGTTCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATACATATAGGCGTTTc                                     >  1:146675/1‑76 (MQ=255)

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

                                    aTAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATCGGCGTTTCCTg                                  <  1:60456/76‑1 (MQ=255)

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

                                      agagagCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                >  1:138953/1‑76 (MQ=255)

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

                                          tGCTTAAGGGCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt                            <  1:1480405/76‑1 (MQ=255)

                                               aaGGGCTTCGACGCGAAGCGTACGCAATCAATATGTAAACAGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                       <  2:114480/76‑1 (MQ=255)

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

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

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

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

                                                  ggCTTCTACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAg                    <  1:1031996/76‑1 (MQ=255)

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

                                                      tCTACGCTAAGCGTGGGCAATGTATATGTAAACGGAGCATTCATATAGGCGTTTCGTGAAAACTATTTCAAGCGAt                <  2:146675/76‑1 (MQ=255)

                                                        tACGCTAAGCGTGGNCAATGAATATGTAAACGGAGCATTCATTTAGGCGTTTCCTGAAAACTATTTCAAGCGACAg              <  1:1570029/76‑1 (MQ=255)

                                                        tACGCTAAGCGTGGNCAATGAATATGTAAACGGAGCATTCATTTAGGCGTTTCCTGAAAAATAGTTCAAGCGATAg              <  1:1570053/76‑1 (MQ=255)

                                                        tACGCTAAGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAg              >  1:1099063/1‑76 (MQ=255)

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

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

                                                            cTAAGCGTGGGCAATGAATATGTCAACNGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa          >  2:5058/1‑76 (MQ=255)

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

                                                              aaGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGc        <  1:691314/76‑1 (MQ=255)

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

                                                               aGCGTGGGCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGATAACTATTTCAAGCGATAGTAAAGCt       >  1:285304/1‑76 (MQ=255)

                                                                gCGTGGGCAATGAATATGTAAACGGAGTATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTc      <  1:795727/76‑1 (MQ=255)

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

                                                                    gggCAATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg  <  1:186278/76‑1 (MQ=255)

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGG  >  REL606/3497938‑3498081

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