BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,179,127	C→A	N4K (AAC→AAA)	yjaG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,179,127	0	C	A	99.2%	482.1 / ‑2.8	122	N4K (AAC→AAA)	yjaG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (78/43); total (79/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 128 total aligned reads displayed.

ATACAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAACCCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGCGAACGC  >  REL606/4179054‑4179199

ataCAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAacc                                                                        <  1:3392170/76‑1 (MQ=255)

 taCAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCa                                                                       >  1:3081161/1‑76 (MQ=255)

  aCAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCaa                                                                      >  1:4715229/1‑76 (MQ=255)

   cAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAAt                                                                     >  2:1795569/1‑76 (MQ=255)

   cAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAAt                                                                     >  1:1777099/1‑76 (MQ=255)

    aGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAAtt                                                                    <  2:40494/76‑1 (MQ=255)

    aGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAAtt                                                                    <  2:1570478/76‑1 (MQ=255)

     gCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTc                                                                   <  1:3463917/76‑1 (MQ=255)

          tCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTg                                                              <  1:2793400/76‑1 (MQ=255)

          tCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTg                                                              >  1:3754076/1‑76 (MQ=255)

          tCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTg                                                              >  1:719461/1‑76 (MQ=255)

            aGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCg                                                            <  1:4542486/76‑1 (MQ=255)

               ccTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCt                                                         >  2:1159052/1‑76 (MQ=255)

               ccTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCt                                                         >  2:382011/1‑76 (MQ=255)

                cTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAACCCAATTCATCTGCGTCTg                                                        >  1:4094254/1‑76 (MQ=255)

                cTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTg                                                        >  2:3494918/1‑76 (MQ=255)

                cTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATATTACAAAAACCAATTCATCTGCGTCTg                                                        >  2:856252/1‑76 (MQ=255)

                 tAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTgg                                                       >  2:3842138/1‑76 (MQ=255)

                   aTTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAg                                                     <  2:1853162/76‑1 (MQ=255)

                   aTTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAg                                                     >  2:3839952/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  <  1:2005488/76‑1 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  1:1715807/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  2:2575461/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  2:3247803/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  1:2466994/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  1:3384931/1‑76 (MQ=255)

                      cAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAgcgc                                                  >  2:729885/1‑76 (MQ=255)

                       aGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGcc                                                 >  2:4955595/1‑76 (MQ=255)

                       aGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGcc                                                 >  1:3086399/1‑76 (MQ=255)

                        ggTAAACTGCGGATAATTTCCTTCTTTGAGAGCTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCt                                                <  2:744585/76‑1 (MQ=255)

                        ggTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCt                                                <  1:3266273/76‑1 (MQ=255)

                        ggTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCt                                                <  1:3765936/76‑1 (MQ=255)

                         gTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTg                                               <  2:2664182/76‑1 (MQ=255)

                         gTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTg                                               <  1:2034150/76‑1 (MQ=255)

                          tAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTgg                                              >  2:1242316/1‑76 (MQ=255)

                          tAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTgg                                              >  2:3892192/1‑76 (MQ=255)

                          tAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTgg                                              <  1:2636722/76‑1 (MQ=255)

                           aaaCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGa                                             >  1:2888549/1‑76 (MQ=255)

                           aaaCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGa                                             >  2:2144245/1‑76 (MQ=255)

                           aaaCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGa                                             >  2:4108154/1‑76 (MQ=255)

                             aCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGaaa                                           <  1:1010546/76‑1 (MQ=255)

                             aCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGaaa                                           <  1:862008/76‑1 (MQ=255)

                             aCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGaaa                                           <  1:3725715/76‑1 (MQ=255)

                             aCTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGaaa                                           <  1:4758284/76‑1 (MQ=255)

                              cTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAg                                          >  1:2524197/1‑76 (MQ=255)

                              cTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAg                                          >  1:815721/1‑76 (MQ=255)

                              cTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAg                                          <  2:2235708/76‑1 (MQ=255)

                              cTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAg                                          >  1:4394010/1‑76 (MQ=255)

                               tGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGc                                         >  1:1583769/1‑76 (MQ=255)

                                gCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCt                                        <  1:2449146/76‑1 (MQ=255)

                                 cGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTg                                       <  2:3816508/76‑1 (MQ=255)

                                 cGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTg                                       >  2:4355779/1‑76 (MQ=255)

                                  ggATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTgg                                      <  1:3433479/76‑1 (MQ=255)

                                  ggATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTgg                                      >  2:487982/1‑76 (MQ=255)

                                  ggATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTgg                                      <  2:1561902/76‑1 (MQ=255)

                                   gATAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgc                                     >  2:3417585/1‑76 (MQ=255)

                                    aTAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgca                                    >  1:1960729/1‑76 (MQ=255)

                                    aTAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgca                                    >  2:3092454/1‑76 (MQ=255)

                                     tAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcag                                   >  2:3623028/1‑76 (MQ=255)

                                     tAATTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcag                                   >  2:2401835/1‑76 (MQ=255)

                                      aaTTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcagc                                  >  2:3446204/1‑76 (MQ=255)

                                      aaTTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcagc                                  >  1:929774/1‑76 (MQ=255)

                                      aaTTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcagc                                  >  1:120038/1‑76 (MQ=255)

                                       aTTTCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGgcagca                                 <  2:3567019/76‑1 (MQ=255)

                                         ttCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATg                               >  1:3188474/1‑76 (MQ=255)

                                         ttCCGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATg                               <  2:4279267/76‑1 (MQ=255)

                                           ccGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTc                             >  2:4947182/1‑76 (MQ=255)

                                            cGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCa                            >  1:2748043/1‑76 (MQ=255)

                                            cGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCa                            >  2:4882085/1‑76 (MQ=255)

                                            cGTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCa                            >  2:4102819/1‑76 (MQ=255)

                                             gTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCAc                           >  2:3239106/1‑76 (MQ=255)

                                             gTATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCAc                           >  1:2512949/1‑76 (MQ=255)

                                              tATTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACt                          >  1:4677034/1‑76 (MQ=255)

                                               aTTTGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACtt                         >  2:1522416/1‑76 (MQ=255)

                                                tttGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACttt                        >  2:1179736/1‑76 (MQ=255)

                                                 ttGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTc                       <  2:3913669/76‑1 (MQ=255)

                                                 ttGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTc                       <  2:1386512/76‑1 (MQ=255)

                                                 ttGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTc                       >  1:4671757/1‑76 (MQ=255)

                                                  tGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCa                      >  1:3300490/1‑76 (MQ=255)

                                                  tGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCa                      <  2:4419835/76‑1 (MQ=255)

                                                  tGAGAACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCa                      >  1:4619654/1‑76 (MQ=255)

                                                   gagaACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCAt                     >  2:4629914/1‑76 (MQ=255)

                                                    agaACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATg                    >  1:3338199/1‑76 (MQ=255)

                                                    agaACTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATg                    <  1:3790688/76‑1 (MQ=255)

                                                      aaCTCAACATGTTACAAAAACCAATTCATCTGCGTATGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGgc                  <  2:2574817/76‑1 (MQ=255)

                                                        cTCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGgctt                <  2:1784725/76‑1 (MQ=255)

                                                         tCAACATGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGgcttt               >  1:879286/1‑75 (MQ=255)

                                                              aTGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTAt          <  2:2700190/76‑1 (MQ=255)

                                                               tGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATg         <  1:4853063/76‑1 (MQ=255)

                                                               tGTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATg         <  1:605939/76‑1 (MQ=255)

                                                                gTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTTCTTTTgt        >  1:4281537/1‑75 (MQ=255)

                                                                gTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGc        <  2:2179475/76‑1 (MQ=255)

                                                                gTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGc        <  2:1218656/76‑1 (MQ=255)

                                                                gTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGc        <  1:1800261/76‑1 (MQ=255)

                                                                gTTACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGc        <  1:657269/76‑1 (MQ=255)

                                                                  tACAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTTTTTTTGCGa      >  1:1506565/1‑76 (MQ=255)

                                                                   aCAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTTCTTTTTCGaa     >  2:3624816/1‑76 (MQ=255)

                                                                    cAAAAACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGCGAAc    <  1:1275553/76‑1 (MQ=255)

                                                                     aaaaACCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTTCTTTTTCGAACg   >  2:4663107/1‑76 (MQ=255)

                                                                      aaaaCCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTTTTTGAACGc  >  2:3456491/1‑76 (MQ=255)

ATACAGCAAATCAGCCCTAATTCAGGTAAACTGCGGATAATTTCCGTATTTGAGAACTCAACATGTTACAAAACCCAATTCATCTGCGTCTGGAGCGCCTGGAAAGCTGGCAGCATGTCACTTTCATGGCTTGCTTATGCGAACGC  >  REL606/4179054‑4179199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`