BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,339,403	T→C	M81T (ATG→ACG)	fis →	DNA‑binding protein Fis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,339,403	0	T	C	98.1%	673.1 / 1.1	161	M81T (ATG→ACG)	fis	DNA‑binding protein Fis
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (69/89); total (70/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 166 total aligned reads displayed.

GAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTT  >  REL606/3339330‑3339478

gAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACgg                                                                           <  2:2959090/76‑1 (MQ=255)

gAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACgg                                                                           <  2:4292086/76‑1 (MQ=255)

gAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACgg                                                                           >  2:3164779/1‑76 (MQ=255)

  aGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGc                                                                         <  1:1131600/76‑1 (MQ=255)

    tAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCAt                                                                       >  2:2207787/1‑76 (MQ=255)

      gAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCa                                                                     <  2:2291591/76‑1 (MQ=255)

      gAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCa                                                                     <  1:1420468/76‑1 (MQ=255)

      gAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCa                                                                     <  2:3205230/76‑1 (MQ=255)

        aCAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAAc                                                                   <  1:1534757/76‑1 (MQ=255)

         cAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAAcc                                                                  <  2:2272693/76‑1 (MQ=255)

         cAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAAcc                                                                  <  2:851525/76‑1 (MQ=255)

            ccccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTg                                                               >  1:2117489/1‑76 (MQ=255)

             cccTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTgg                                                              >  2:83549/1‑76 (MQ=255)

               cTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTa                                                            >  1:1996426/1‑76 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           <  2:1501318/76‑1 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           <  2:2085510/76‑1 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           <  2:1853808/76‑1 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           >  2:1591559/1‑76 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           >  1:1040084/1‑76 (MQ=255)

                tGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTAc                                                           >  1:3305180/1‑76 (MQ=255)

                 gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACg                                                          <  1:2051612/76‑1 (MQ=255)

                 gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACg                                                          <  1:822123/76‑1 (MQ=255)

                 gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACg                                                          <  1:2213759/76‑1 (MQ=255)

                 gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACg                                                          >  1:3285515/1‑76 (MQ=255)

                 gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACg                                                          >  2:4357958/1‑76 (MQ=255)

                  ttGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGc                                                         <  2:4166814/76‑1 (MQ=255)

                   tGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCt                                                        <  1:4563001/76‑1 (MQ=255)

                   tGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCt                                                        >  2:1523295/1‑76 (MQ=255)

                    ggACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTg                                                       >  2:2561744/1‑76 (MQ=255)

                    ggACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTg                                                       <  2:4959912/76‑1 (MQ=255)

                     gACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGc                                                      >  1:27053/1‑76 (MQ=255)

                        aTGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTa                                                   <  1:864485/76‑1 (MQ=255)

                         tGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaa                                                  <  2:4414377/76‑1 (MQ=255)

                         tGGTGATGCAAGACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaa                                                  <  1:1932727/76‑1 (MQ=255)

                          ggTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaaa                                                 <  2:1241670/76‑1 (MQ=255)

                          ggTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaaa                                                 <  2:2856461/76‑1 (MQ=255)

                          ggTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaaa                                                 <  2:4614937/76‑1 (MQ=255)

                           gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTaaaa                                                <  2:2404082/76‑1 (MQ=255)

                           gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATTAACCGTGGTACGCTGCGTaaaa                                                <  2:3462194/76‑1 (MQ=255)

                           gTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTaaaa                                                <  1:432400/76‑1 (MQ=255)

                              aTGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAAt                                             <  1:187127/76‑1 (MQ=255)

                                gCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTg                                           <  1:4366581/76‑1 (MQ=255)

                                 cAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGa                                          >  2:1595463/1‑76 (MQ=255)

                                 cAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGa                                          >  2:947131/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCCGACCCGTGATGCGCAGCTGACGCACATCACCAGTGGTACGCTGCGCACAACATCGaa                                         >  2:2700194/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         >  2:2700193/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         >  2:2170539/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         >  2:3245917/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         >  2:3899824/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         >  2:406320/1‑76 (MQ=255)

                                  aaTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaa                                         <  2:417563/76‑1 (MQ=255)

                                   aTACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaaa                                        <  1:1816937/76‑1 (MQ=255)

                                    tACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaaaa                                       <  2:2760878/76‑1 (MQ=255)

                                    tACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaaaa                                       <  2:923002/76‑1 (MQ=255)

                                     acacCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaaaaa                                      <  1:778230/76‑1 (MQ=255)

                                     acacCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGaaaaa                                      <  1:2619865/76‑1 (MQ=255)

                                       acCCGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAAt                                    <  2:549102/76‑1 (MQ=255)

                                         ccGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATAc                                  <  1:2880292/76‑1 (MQ=255)

                                          cGTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACg                                 <  2:2993074/76‑1 (MQ=255)

                                           gTGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACgg                                >  2:458421/1‑76 (MQ=255)

                                            tGGTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGc                               <  1:139265/76‑1 (MQ=255)

                                             ggTAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCa                              <  1:2296428/76‑1 (MQ=255)

                                               tAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATg                            <  2:4557912/76‑1 (MQ=255)

                                               tAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATg                            >  1:2954521/1‑76 (MQ=255)

                                               tAACCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATg                            <  1:2115206/76‑1 (MQ=255)

                                                aaCCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGa                           >  1:1666872/1‑76 (MQ=255)

                                                aaCCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGa                           >  1:1143292/1‑76 (MQ=255)

                                                aaCCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGa                           >  1:1539205/1‑76 (MQ=255)

                                                 aCCAGATCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGaa                          <  2:3752734/76‑1 (MQ=255)

                                                 aCCAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGaa                          <  2:3372913/76‑1 (MQ=255)

                                                  ccAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAAc                         <  1:4171303/76‑1 (MQ=255)

                                                   cAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACt                        >  2:2359145/1‑76 (MQ=255)

                                                   cAGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACt                        >  1:4999909/1‑76 (MQ=255)

                                                    aGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTa                       <  1:4614956/76‑1 (MQ=255)

                                                    aGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTa                       >  2:3942463/1‑76 (MQ=255)

                                                    aGACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTa                       <  2:4460063/76‑1 (MQ=255)

                                                     gACCCGTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTaa                      <  2:4966756/76‑1 (MQ=255)

                                                          gTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAg                 <  2:3489419/76‑1 (MQ=255)

                                                          gTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAg                 >  2:679466/1‑76 (MQ=255)

                                                          gTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAg                 >  1:3013397/1‑76 (MQ=255)

                                                          gTGCTGCGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAg                 <  1:1318764/76‑1 (MQ=255)

                                                           tgctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAgg                <  1:863337/76‑1 (MQ=255)

                                                           tgctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAgg                >  1:4436003/1‑76 (MQ=255)

                                                           tgctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAgg                >  2:695244/1‑76 (MQ=255)

                                                            gctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGt               <  1:5015383/76‑1 (MQ=255)

                                                             ctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGtt              <  1:662355/76‑1 (MQ=255)

                                                             ctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGtt              <  1:2661872/76‑1 (MQ=255)

                                                             ctgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGtt              <  1:1489411/76‑1 (MQ=255)

                                                              tgcGCTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTa             <  1:1276600/76‑1 (MQ=255)

                                                               gcgcTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAg            <  2:2996170/76‑1 (MQ=255)

                                                               gcgcTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAg            <  1:3155722/76‑1 (MQ=255)

                                                                 gcTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCt          <  1:4923775/76‑1 (MQ=255)

                                                                 gcTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCACGAACTAATTCAGGTTAGCt          <  1:588421/76‑1 (MQ=255)

                                                                  cTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTa         >  1:4003345/1‑76 (MQ=255)

                                                                  cTGATGACGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTa         <  2:2881520/76‑1 (MQ=255)

                                                                   tgatgaCGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTaa        >  1:4145765/1‑76 (MQ=255)

                                                                   tgatgaCGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTaa        >  1:1961968/1‑76 (MQ=255)

                                                                   tgatgaCGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTaa        <  1:4254156/76‑1 (MQ=255)

                                                                    gatgacGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTaaa       >  2:2333027/1‑76 (MQ=255)

                                                                      tgacgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATg     >  1:2471072/1‑76 (MQ=255)

GAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTT  >  REL606/3339330‑3339478

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