BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,498,004	G→A	T11M (ACG→ATG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,498,004	0	G	A	96.9%	245.9 / ‑1.1	65	T11M (ACG→ATG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (30/33); total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAG  >  REL606/3497931‑3498074

gacgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATg                                                                      >  2:1389235/1‑76 (MQ=255)

    acCGGGTCCTTTCCAATCGTCAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCg                                                                  <  1:1081511/76‑1 (MQ=255)

    acCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCg                                                                  <  1:1658621/76‑1 (MQ=255)

     ccGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGa                                                                 >  2:784113/1‑76 (MQ=255)

       gggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATg                                                               <  2:458361/76‑1 (MQ=255)

        ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGa                                                              <  1:1015104/76‑1 (MQ=255)

        ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGa                                                              >  1:1264668/1‑76 (MQ=255)

         gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGaa                                                             <  2:971491/76‑1 (MQ=255)

          tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAat                                                            >  1:325798/1‑76 (MQ=255)

          tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAat                                                            <  1:839801/76‑1 (MQ=255)

          tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAat                                                            >  1:517603/1‑76 (MQ=255)

              ttCCAATCGTAAACATCAGTTTGCAAGCGATAGAGGGATTAAGGGCTTCTGCGCTAAGCATGGGCGATGAATATGt                                                        <  2:394647/76‑1 (MQ=255)

                ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTaa                                                      >  1:359871/1‑76 (MQ=255)

                 cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTaaa                                                     <  2:1085410/76‑1 (MQ=255)

                  aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAAc                                                    <  1:80310/76‑1 (MQ=255)

                    tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACgg                                                  <  2:382801/76‑1 (MQ=255)

                    tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACgg                                                  >  2:470098/1‑76 (MQ=255)

                      gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAg                                                >  1:458451/1‑76 (MQ=255)

                      gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAg                                                >  1:500991/1‑76 (MQ=255)

                      gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAg                                                >  1:1215218/1‑76 (MQ=255)

                        aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCa                                              <  1:1757539/76‑1 (MQ=255)

                        aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCa                                              >  1:208134/1‑76 (MQ=255)

                        aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCa                                              >  1:1776531/1‑76 (MQ=255)

                            aTCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCa                                          >  1:1600118/1‑76 (MQ=255)

                               aGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCatat                                       <  2:1727456/76‑1 (MQ=255)

                               aGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCatat                                       <  2:1418262/76‑1 (MQ=255)

                                 cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAg                                     <  2:1250749/76‑1 (MQ=255)

                                 cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAg                                     <  1:1103278/76‑1 (MQ=255)

                                 cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAg                                     >  1:1346124/1‑76 (MQ=255)

                                  ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAgg                                    <  1:1046100/76‑1 (MQ=255)

                                  ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAgg                                    <  2:1134339/76‑1 (MQ=255)

                                   tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                   <  2:558343/76‑1 (MQ=255)

                                   tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                   >  2:687080/1‑76 (MQ=255)

                                    gTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCg                                  <  2:65255/76‑1 (MQ=255)

                                       aGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGttt                               <  1:1335998/76‑1 (MQ=255)

                                        gCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                              <  2:1707389/76‑1 (MQ=255)

                                          gATAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt                            <  1:320419/76‑1 (MQ=255)

                                           aTAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTccag                           >  1:642538/1‑74 (MQ=255)

                                           aTAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                           >  1:175718/1‑76 (MQ=255)

                                           aTAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                           >  2:1020951/1‑76 (MQ=255)

                                            tAGAGTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                          <  2:202371/76‑1 (MQ=255)

                                             agagTGCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                         >  1:125425/1‑76 (MQ=255)

                                                  gCTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                    <  1:29396/76‑1 (MQ=255)

                                                   cTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                   >  2:833870/1‑76 (MQ=255)

                                                   cTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                   >  1:1596742/1‑76 (MQ=255)

                                                   cTTAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                   >  2:646975/1‑76 (MQ=255)

                                                     tAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                 <  1:1000469/76‑1 (MQ=255)

                                                     tAAGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                 >  1:837990/1‑76 (MQ=255)

                                                       aGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa               <  1:1627605/76‑1 (MQ=255)

                                                       aGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa               <  1:1692039/76‑1 (MQ=255)

                                                       aGGGCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa               <  2:359871/76‑1 (MQ=255)

                                                        gggCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa              <  1:862940/76‑1 (MQ=255)

                                                         ggCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAg             <  2:1513103/76‑1 (MQ=255)

                                                          gCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc            <  1:1806258/76‑1 (MQ=255)

                                                          gCTTCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc            >  1:1068396/1‑76 (MQ=255)

                                                            ttCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa          >  2:1145758/1‑76 (MQ=255)

                                                            ttCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa          >  1:802714/1‑76 (MQ=255)

                                                             tCTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt         >  2:1006496/1‑76 (MQ=255)

                                                              cTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa        >  1:711036/1‑76 (MQ=255)

                                                              cTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa        <  1:863509/76‑1 (MQ=255)

                                                              cTACGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa        <  1:889755/76‑1 (MQ=255)

                                                                 cGCTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTa     <  1:602145/76‑1 (MQ=255)

                                                                   cTAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa   >  2:854379/1‑76 (MQ=255)

                                                                    tAAGCTTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGCTTCCTGAAAACTATTTCAAGCGATAGTaaac  >  2:356691/1‑75 (MQ=255)

                                                                    tAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg  >  1:390612/1‑76 (MQ=255)

                                                                    tAAGCATGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg  <  1:1506783/76‑1 (MQ=255)

GACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAG  >  REL606/3497931‑3498074

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