BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,116,047	(C)_7→8	coding (345/1218 nt)	ECB_02003 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,116,040	1	.	C	97.1%	734.9 / 10.3	171	coding (352/1218 nt)	ECB_02003	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/3); new base C (81/85); total (83/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

CTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑C‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGTATACA  >  REL606/2115965‑2116111

cTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTaa                                                                            >  2:3023784/1‑76 (MQ=255)

cTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTaa                                                                            <  1:4855864/76‑1 (MQ=255)

cTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTaa                                                                            >  1:4214187/1‑76 (MQ=255)

cTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTaa                                                                            >  2:1104256/1‑76 (MQ=255)

 tAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c                                                                          >  2:7393080/1‑76 (MQ=255)

  aaCGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c                                                                        >  1:2829668/1‑76 (MQ=255)

  aaCGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c                                                                        <  2:4190728/76‑1 (MQ=255)

  aaCGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c                                                                        <  1:3785947/76‑1 (MQ=255)

   aCGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑c                                                                      <  1:4770936/76‑1 (MQ=255)

   aCGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑c                                                                      <  1:6911656/76‑1 (MQ=255)

     gAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccc                                                                    <  2:5672499/76‑1 (MQ=255)

     gAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccc                                                                    >  2:5655642/1‑76 (MQ=255)

     gAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccc                                                                    <  2:6806952/76‑1 (MQ=255)

      aaaTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccc                                                                   >  2:6460647/1‑76 (MQ=255)

       aaTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccccc                                                                  >  2:1300590/1‑76 (MQ=255)

       aaTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccccc                                                                  <  2:4357589/76‑1 (MQ=255)

       aaTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑ccccc                                                                  <  1:3098368/76‑1 (MQ=255)

        aTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccccc                                                                 <  2:6652255/76‑2 (MQ=255)

        aTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccccc                                                                 <  1:5028933/76‑2 (MQ=255)

        aTGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccccc                                                                 >  2:7901489/1‑75 (MQ=255)

         tGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccccca                                                                <  1:6118930/76‑3 (MQ=255)

         tGAATCCAATTAAATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑c‑c‑cccccca                                                                <  2:632398/76‑3 (MQ=255)

             tCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACt                                                             >  2:3350011/1‑76 (MQ=255)

             tCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACt                                                             <  1:7079404/76‑1 (MQ=255)

             tCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACt                                                             <  2:4377062/76‑1 (MQ=255)

             tCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACt                                                             <  2:7132552/76‑1 (MQ=255)

             tCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACt                                                             <  1:2978689/76‑1 (MQ=255)

               cAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTAt                                                           <  2:397057/76‑1 (MQ=255)

               cAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTAt                                                           >  1:6205983/1‑76 (MQ=255)

                aaTTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATg                                                          <  2:4027587/76‑1 (MQ=255)

                aaTTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATg                                                          <  1:6501766/76‑1 (MQ=255)

                aaTTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATg                                                          <  1:1403111/76‑1 (MQ=255)

                   tataTTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGc                                                       <  2:5994853/76‑1 (MQ=255)

                    atatTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCg                                                      >  2:3098053/1‑76 (MQ=255)

                      atTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGat                                                    >  1:6166363/1‑76 (MQ=255)

                       ttAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatt                                                   >  2:1719928/1‑76 (MQ=255)

                        tAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatta                                                  <  1:4336461/76‑1 (MQ=255)

                        tAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatta                                                  >  2:2855133/1‑76 (MQ=255)

                        tAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatta                                                  >  1:5761188/1‑76 (MQ=255)

                        tAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatta                                                  <  2:5103068/76‑1 (MQ=255)

                        tAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGatta                                                  <  1:5652291/76‑1 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               <  1:7157031/76‑1 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               <  1:6363946/76‑1 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               >  1:2900787/1‑76 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               >  1:2908708/1‑76 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               <  1:7799152/76‑1 (MQ=255)

                           aTTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGAttatta                                               >  1:7692235/1‑76 (MQ=255)

                            ttGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATtattat                                              <  2:6857319/76‑1 (MQ=255)

                             tGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATtata                                             >  1:4106126/1‑76 (MQ=255)

                             tGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATtata                                             >  2:7246630/1‑76 (MQ=255)

                             tGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATtata                                             >  1:2046124/1‑76 (MQ=255)

                             tGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATtata                                             >  1:4896920/1‑76 (MQ=255)

                               ttCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAAt                                           >  2:4811553/1‑76 (MQ=255)

                                 ccTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTc                                         <  2:1938874/76‑1 (MQ=255)

                                 ccTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTc                                         >  2:2064072/1‑76 (MQ=255)

                                  cTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCg                                        <  1:6737716/76‑1 (MQ=255)

                                    cATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCgag                                      <  2:1523400/76‑1 (MQ=255)

                                    cATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCgag                                      >  2:4739275/1‑76 (MQ=255)

                                      tGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGagag                                    >  2:4618356/1‑76 (MQ=255)

                                      tGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGagag                                    <  2:8183288/76‑1 (MQ=255)

                                         tCAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTaa                                 <  1:8267151/76‑1 (MQ=255)

                                         tCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑CTC‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTaa                                 <  2:3533570/76‑1 (MQ=255)

                                          cAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAAt                                >  1:989246/1‑76 (MQ=255)

                                          cAACTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAAt                                >  1:4351750/1‑76 (MQ=255)

                                          cAACTTTTTTTAGAATTAGAATTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAAt                                >  2:7547603/1‑76 (MQ=255)

                                            aCTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCa                              <  1:4216042/76‑1 (MQ=255)

                                            aCTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCa                              <  2:6098303/76‑1 (MQ=255)

                                             cTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCAc                             <  2:3452107/76‑1 (MQ=255)

                                             cTTTTTTTAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCAc                             <  1:3602195/76‑1 (MQ=255)

                                               ttttttAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACtt                           >  2:7276353/1‑76 (MQ=255)

                                               ttttttAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACtt                           <  1:6640551/76‑1 (MQ=255)

                                                  tttAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAgt                        >  2:5231635/1‑76 (MQ=255)

                                                   ttAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAgtg                       >  2:6971919/1‑76 (MQ=255)

                                                    tAGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAgtgt                      <  2:6136224/76‑1 (MQ=255)

                                                    tAGAATTAGATTTTTTGATATTAA‑C‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGtt                     <  1:3961998/76‑1 (MQ=255)

                                                     aGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGtt                     >  2:5921435/1‑76 (MQ=255)

                                                     aGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGtt                     >  2:8006060/1‑76 (MQ=255)

                                                     aGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGtt                     <  1:3663037/76‑1 (MQ=255)

                                                     aGAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGtt                     <  2:1643347/76‑1 (MQ=255)

                                                      gAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTa                    <  2:807799/76‑1 (MQ=255)

                                                      gAATTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTa                    <  1:4488339/76‑1 (MQ=255)

                                                       aaTTAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTAc                   >  1:3275579/1‑76 (MQ=255)

                                                         ttAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTActc                 <  1:4732729/76‑1 (MQ=255)

                                                          tAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTActct                >  2:3074632/1‑76 (MQ=255)

                                                          tAGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTActct                >  2:7473848/1‑76 (MQ=255)

                                                           aGATTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTa               >  1:246215/1‑76 (MQ=255)

                                                             aTTTTTTGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAAt             >  2:2177413/1‑76 (MQ=255)

                                                                ttttGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAAtctc          <  1:3472916/76‑1 (MQ=255)

                                                                 tttGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCa         >  1:5274307/1‑76 (MQ=255)

                                                                 tttGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCa         <  1:2497633/76‑1 (MQ=255)

                                                                 tttGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCa         >  1:7725307/1‑76 (MQ=255)

                                                                 tttGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCa         >  1:7253872/1‑76 (MQ=255)

                                                                   tGATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGt       <  2:552459/76‑1 (MQ=255)

                                                                   tGATATTAA‑C‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGta      <  1:7611502/76‑1 (MQ=255)

                                                                    gATATTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGta      <  2:2952307/76‑1 (MQ=255)

                                                                     atatTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGtat     >  2:7761406/1‑76 (MQ=255)

                                                                     atatTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGtat     >  1:1142443/1‑76 (MQ=255)

                                                                      tatTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGtata    >  2:5398814/1‑76 (MQ=255)

                                                                       atTAACC‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGTATAc   <  2:7879843/76‑1 (MQ=255)

                                                                        ttAA‑C‑CTCCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGTATAct  >  2:6924354/1‑75 (MQ=255)

CTAACGAAATGAATCCAATTATATTAGATTGTTCCTCATGCTCAACTTTTTTTAGAATTAGATTTTTTGATATTAA‑C‑C‑CCCCCATACTATGGGCGATTATTATAATTCGAGAGTAATCACTTAGTGTTACTCTAATCTCAGTATACA  >  REL606/2115965‑2116111

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