BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,275,544	T→C	D135G (GAT→GGT)	yjcO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,275,544	0	T	C	98.3%	239.6 / ‑4.5	60	D135G (GAT→GGT)	yjcO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (18/41); total (18/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTAATGCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGA  >  REL606/4275472‑4275615

ttAATGCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCtt                                                                      >  1:1519674/1‑76 (MQ=255)

 tAATGCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTc                                                                     <  1:50634/76‑1 (MQ=255)

   aTGCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCg                                                                   <  2:815278/76‑1 (MQ=255)

     gCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAg                                                                 >  1:245728/1‑76 (MQ=255)

     gCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAg                                                                 >  1:756078/1‑76 (MQ=255)

     gCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAg                                                                 <  2:1052530/76‑1 (MQ=255)

     gCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAg                                                                 <  2:55813/76‑1 (MQ=255)

        cACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCg                                                              <  1:610965/76‑1 (MQ=255)

        cACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCg                                                              <  1:1563547/76‑1 (MQ=255)

         aCGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGt                                                             >  1:40773/1‑76 (MQ=255)

           gCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGttt                                                           >  1:1181447/1‑76 (MQ=255)

            ccGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGtttt                                                          <  1:302319/76‑1 (MQ=255)

             cGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTc                                                         <  1:437532/76‑1 (MQ=255)

             cGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTc                                                         <  2:1234182/76‑1 (MQ=255)

              gTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCt                                                        <  1:968171/76‑1 (MQ=255)

               ttGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTa                                                       >  2:1077532/1‑76 (MQ=255)

                tGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTaa                                                      >  1:1696094/1‑76 (MQ=255)

                 ggCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAAc                                                     <  1:663316/76‑1 (MQ=255)

                 ggCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAAc                                                     <  1:303582/76‑1 (MQ=255)

                 ggCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAAc                                                     <  2:1232131/76‑1 (MQ=255)

                 ggCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAAc                                                     <  1:100512/76‑1 (MQ=255)

                  gCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACa                                                    <  1:34626/76‑1 (MQ=255)

                   cGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAg                                                   <  1:509840/76‑1 (MQ=255)

                   cGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAg                                                   <  1:1236173/76‑1 (MQ=255)

                     tAAATCAAACCAAGCAACATTTGGGCGTCGACCGCAACGTCGTTCTCCAGACCTTCCCAGCAGTTATCAAACAGca                                                 >  2:1837907/1‑75 (MQ=255)

                      aaaTCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGa                                                >  1:24072/1‑76 (MQ=255)

                       aaTCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGaa                                               <  1:1702730/76‑1 (MQ=255)

                        aTCAAACCAAGCAGCATTTTGGCATCGACGGCAGAGTCGTGCTCCAGGCCTTCCGAGGCGTTTTCTAACAGCGaaa                                              <  1:1118997/76‑1 (MQ=255)

                        aTCAAACCAAGCAGCATTTAGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGaaa                                              <  1:1118995/76‑1 (MQ=255)

                          cAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAAtt                                            <  1:1397370/76‑1 (MQ=255)

                               cAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCtt                                       <  2:1293260/76‑1 (MQ=255)

                                   cagcaTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGa                                   >  2:1693131/1‑76 (MQ=255)

                                   cagcaTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTGTTCTAACAGCGCAAATGCCTTTGGa                                   >  1:963587/1‑76 (MQ=255)

                                    agcaTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGAt                                  <  1:868728/76‑1 (MQ=255)

                                    agcaTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGAt                                  >  2:1811923/1‑76 (MQ=255)

                                      caTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATaa                                <  2:1357532/76‑1 (MQ=255)

                                       aTTTGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAAt                               <  1:908401/76‑1 (MQ=255)

                                        tttGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATc                              >  1:562272/1‑76 (MQ=255)

                                          tGGGCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCg                            >  2:876689/1‑76 (MQ=255)

                                           gggCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCgg                           <  2:990172/76‑1 (MQ=255)

                                           gggCATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCgg                           <  1:1607572/76‑1 (MQ=255)

                                              cATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGttt                        <  2:1379165/76‑1 (MQ=255)

                                              cATCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGttt                        <  2:1690701/76‑1 (MQ=255)

                                               aTCGACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTa                       >  2:451926/1‑76 (MQ=255)

                                                tCGGCGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGAGAATCCGGTTTAc                      <  2:112284/76‑1 (MQ=255)

                                                  gACGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTAccc                    <  1:985813/76‑1 (MQ=255)

                                                    cGGCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGc                  <  1:771466/76‑1 (MQ=255)

                                                      gCAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCtt                <  1:348466/76‑1 (MQ=255)

                                                       cAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTg               <  2:89301/76‑1 (MQ=255)

                                                       cAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTg               <  2:1587869/76‑1 (MQ=255)

                                                       cAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTg               <  2:1746931/76‑1 (MQ=255)

                                                       cAGAGTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTg               >  1:1695481/1‑76 (MQ=255)

                                                        agagTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGa              >  1:1364725/1‑76 (MQ=255)

                                                         gagTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAg             <  1:908562/76‑1 (MQ=255)

                                                          agTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGt            >  2:1847477/1‑76 (MQ=255)

                                                           gTCGTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTa           <  2:1163542/76‑1 (MQ=255)

                                                              gTTCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTAttt        <  2:39999/76‑1 (MQ=255)

                                                               ttCTCCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTa       <  1:1366692/76‑1 (MQ=255)

                                                                 ctcCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTAcc     >  1:548118/1‑76 (MQ=255)

                                                                 ctcCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTAcc     <  1:1661818/76‑1 (MQ=255)

                                                                  tcCAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCa    <  2:864693/76‑1 (MQ=255)

                                                                    cAGACCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGa  >  1:675435/1‑76 (MQ=255)

TTAATGCCCACGCCGTTGGCGTAAATCAAACCAAGCAGCATTTGGGCATCGACGGCAGAGTCGTTCTCCAGATCTTCCGAGGCGTTTTCTAACAGCGAAATTGCCTTTGGATAATCCGGTTTACCCGCTTGAGTATTTACCAGA  >  REL606/4275472‑4275615

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: `‑`