BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,980	T→C	H19R (CAC→CGC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,980	0	T	C	98.8%	325.6 / ‑6.3	80	H19R (CAC→CGC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (36/43); total (37/43)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAA  >  REL606/3497910‑3498051

cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCtt                                                                     <  1:1939670/75‑1 (MQ=255)

 aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTa                                                                    <  2:854801/75‑1 (MQ=255)

 aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTa                                                                    >  2:2030820/1‑75 (MQ=255)

 aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTa                                                                    >  2:74574/1‑75 (MQ=255)

    cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAgg                                                                 >  1:572565/1‑75 (MQ=255)

    cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAgg                                                                 <  2:4592/75‑1 (MQ=255)

      ttCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGc                                                               <  2:851588/75‑1 (MQ=255)

       tCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCt                                                              >  1:26280/1‑75 (MQ=255)

       tCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCt                                                              <  2:438361/75‑1 (MQ=255)

        cATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCtt                                                             >  1:195293/1‑75 (MQ=255)

          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTACGGGCTTCt                                                           >  2:740391/1‑75 (MQ=255)

          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCt                                                           <  1:481224/75‑1 (MQ=255)

          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCt                                                           >  2:473958/1‑75 (MQ=255)

          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCt                                                           <  2:1825433/75‑1 (MQ=255)

              tATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTActc                                                       >  1:1710231/1‑73 (MQ=255)

              tATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGc                                                       <  1:943037/75‑1 (MQ=255)

                ttGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTa                                                     >  2:320415/1‑75 (MQ=255)

                 tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTaa                                                    <  2:1113113/75‑1 (MQ=255)

                 tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTaa                                                    <  2:925493/75‑1 (MQ=255)

                 tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTaa                                                    <  2:213982/75‑1 (MQ=255)

                  ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAg                                                   >  2:1382859/1‑75 (MQ=255)

                  ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAg                                                   <  1:1041570/75‑1 (MQ=255)

                    cgacgaCCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCACTTAAGGGCTTCTACGCTAAGCg                                                 <  2:1387948/75‑1 (MQ=255)

                     gacgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGt                                                >  2:1031114/1‑75 (MQ=255)

                      acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTg                                               <  2:1139012/75‑1 (MQ=255)

                      acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTg                                               >  2:140498/1‑75 (MQ=255)

                        gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTggg                                             <  1:320415/75‑1 (MQ=255)

                        gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTggg                                             >  2:1242422/1‑75 (MQ=255)

                         acCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGc                                            <  1:1991369/75‑1 (MQ=255)

                          ccGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCg                                           >  2:1256063/1‑75 (MQ=255)

                           cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGa                                          >  2:2064345/1‑75 (MQ=255)

                           cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGa                                          <  2:540635/75‑1 (MQ=255)

                             ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATg                                        <  1:286779/75‑1 (MQ=255)

                             ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATg                                        <  1:286773/75‑1 (MQ=255)

                              gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGa                                       <  2:1291/75‑1 (MQ=255)

                               tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGaa                                      >  1:560269/1‑75 (MQ=255)

                                ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAat                                     <  2:418865/75‑1 (MQ=255)

                                 cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAata                                    <  2:1197199/75‑1 (MQ=255)

                                  tttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAatat                                   <  2:365010/75‑1 (MQ=255)

                                     ccAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTa                                <  1:1878385/75‑1 (MQ=255)

                                      cAATCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaa                               <  2:1397966/75‑1 (MQ=255)

                                       aaTCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                              >  1:1878270/1‑75 (MQ=255)

                                       aaTCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                              >  1:935822/1‑75 (MQ=255)

                                       aaTCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                              >  2:381163/1‑75 (MQ=255)

                                        aTCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAAc                             >  1:898600/1‑75 (MQ=255)

                                        aTCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAAc                             >  1:1231496/1‑75 (MQ=255)

                                         tCGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACg                            <  1:564119/75‑1 (MQ=255)

                                          cGTAAACATCATCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACgg                           <  2:673191/75‑1 (MQ=255)

                                          cGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAGGCGTGGGCGATGAATATGTAAACgg                           <  1:1527799/75‑1 (MQ=255)

                                          cGTAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACgg                           <  1:668959/75‑1 (MQ=255)

                                            tAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAg                         <  1:1371085/75‑1 (MQ=255)

                                            tAAACATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAg                         >  2:3394/1‑75 (MQ=255)

                                             aaaCATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGc                        >  2:1327158/1‑75 (MQ=255)

                                             aaaCATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGc                        >  2:1257148/1‑75 (MQ=255)

                                             aaaCATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGc                        <  1:316300/75‑1 (MQ=255)

                                             aaaCATCAAGTTGGAAGTGATAGAGGGGTTAAAGGGTTCTTAGGTAAACCTGGGCGATGAAGATGTAAACGGAGc                        <  1:316323/75‑1 (MQ=255)

                                               aCATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCAt                      <  1:1612957/75‑1 (MQ=255)

                                               aCATCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCAt                      <  2:1951170/75‑1 (MQ=255)

                                                 aTCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTc                    >  1:1598088/1‑75 (MQ=255)

                                                 aTCAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTc                    >  1:378569/1‑75 (MQ=255)

                                                   cAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCat                  >  1:387166/1‑75 (MQ=255)

                                                   cAGCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCat                  <  1:1793893/75‑1 (MQ=255)

                                                     gcctGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                <  2:1594220/72‑1 (MQ=255)

                                                     gCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                >  2:431714/1‑75 (MQ=255)

                                                     gCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                >  2:971873/1‑75 (MQ=255)

                                                     gCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                <  2:1731939/75‑1 (MQ=255)

                                                     gCTTGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                <  1:870126/75‑1 (MQ=255)

                                                        tGTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAgg             >  2:1514309/1‑75 (MQ=255)

                                                         gTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc            <  2:1595731/75‑1 (MQ=255)

                                                         gTAAGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc            <  1:1675309/75‑1 (MQ=255)

                                                           aaGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGt          >  2:1520685/1‑75 (MQ=255)

                                                           aaGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGt          >  1:1975999/1‑75 (MQ=255)

                                                            aGCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGtt         >  2:1345541/1‑75 (MQ=255)

                                                             gCGATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGttt        >  1:843790/1‑75 (MQ=255)

                                                               gATAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc      >  1:477550/1‑75 (MQ=255)

                                                                aTAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt     <  1:2052435/75‑1 (MQ=255)

                                                                aTAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt     <  1:1777999/75‑1 (MQ=255)

                                                                aTAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt     >  1:1177987/1‑75 (MQ=255)

                                                                 tAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg    <  1:1327158/75‑1 (MQ=255)

                                                                 tAGAGCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg    <  1:1234871/75‑1 (MQ=255)

                                                                   gagCGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa  >  2:57031/1‑75 (MQ=255)

CAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAA  >  REL606/3497910‑3498051

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