BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	REL606	70,867	T→C	100%	D92G (GAC→GGC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	70,867	0	T	C	100.0%	414.4 / NA	108	D92G (GAC→GGC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (59/49); total (59/49)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 108 total aligned reads displayed.

TTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGTCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACC  >  REL606/70819‑70914

ttGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCg                                               >  1:4229653/1‑51 (MQ=255)

 tGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGt                                              >  1:2233478/1‑51 (MQ=255)

 tGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGt                                              <  1:2016056/51‑1 (MQ=255)

 tGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGt                                              >  1:7929324/1‑51 (MQ=255)

  gAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGtt                                             >  1:5127005/1‑51 (MQ=255)

  gAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGtt                                             <  1:4359438/51‑1 (MQ=255)

    aCTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGa                                           >  1:4539739/1‑51 (MQ=255)

    aCTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGa                                           >  1:208898/1‑51 (MQ=255)

     cTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGAt                                          >  1:10059959/1‑51 (MQ=255)

      tGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATc                                         <  1:10760537/51‑1 (MQ=255)

      tGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATc                                         >  1:4341316/1‑51 (MQ=255)

       gggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATcc                                        >  1:2278424/1‑51 (MQ=255)

        ggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCca                                       >  1:10759044/1‑51 (MQ=255)

        ggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCca                                       >  1:6111081/1‑51 (MQ=255)

        ggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCca                                       >  1:6343818/1‑51 (MQ=255)

        ggTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCca                                       <  1:1598904/51‑1 (MQ=255)

         gtgtGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCcac                                      >  1:10008889/1‑51 (MQ=255)

         gtgtGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCcac                                      <  1:4787997/51‑1 (MQ=255)

          tgtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCcaca                                     <  1:2259794/51‑1 (MQ=255)

          tgtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCcaca                                     >  1:628919/1‑51 (MQ=255)

          tgtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCATTGATCcaca                                     <  1:2126858/51‑1 (MQ=38)

           gtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAt                                    <  1:9587399/51‑1 (MQ=255)

           gtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAt                                    <  1:10384867/51‑1 (MQ=255)

           gtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAt                                    <  1:950276/51‑1 (MQ=255)

           gtgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAt                                    <  1:2652697/51‑1 (MQ=255)

            tgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtt                                   >  1:3878232/1‑51 (MQ=255)

            tgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtt                                   >  1:5856937/1‑51 (MQ=255)

            tgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtt                                   <  1:2503996/51‑1 (MQ=255)

            tgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtt                                   <  1:9527728/51‑1 (MQ=255)

            tgGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtt                                   >  1:9019953/1‑51 (MQ=255)

             ggAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAttt                                  >  1:4367447/1‑51 (MQ=255)

              gAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtttt                                 >  1:8585587/1‑51 (MQ=255)

              gAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtttt                                 <  1:2178896/51‑1 (MQ=255)

               aaTTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTg                                <  1:7774075/51‑1 (MQ=255)

               aaTTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTg                                <  1:5347544/51‑1 (MQ=255)

                 ttGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGACCCACATTTTGGc                              >  1:5345591/1‑51 (MQ=38)

                  tGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGcc                             >  1:7014715/1‑51 (MQ=255)

                  tGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGcc                             <  1:6839286/51‑1 (MQ=255)

                  tGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGcc                             >  1:6773857/1‑51 (MQ=255)

                  tGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGcc                             <  1:7148563/51‑1 (MQ=255)

                   gcagcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCg                            >  1:651390/1‑51 (MQ=255)

                   gcagcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCg                            >  1:5638186/1‑51 (MQ=255)

                    cagcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgg                           <  1:8135569/51‑1 (MQ=255)

                    cagcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgg                           >  1:8225546/1‑51 (MQ=255)

                     agcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCggg                          >  1:3914877/1‑51 (MQ=255)

                     agcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCggg                          <  1:8689171/51‑1 (MQ=255)

                      gcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                         >  1:7116365/1‑51 (MQ=255)

                      gcaACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                         >  1:1656996/1‑51 (MQ=255)

                       caACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGga                        >  1:608910/1‑51 (MQ=255)

                       caACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGga                        >  1:7643138/1‑51 (MQ=255)

                       caACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGga                        <  1:6325637/51‑1 (MQ=255)

                         aCGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGgaga                      <  1:3572657/51‑1 (MQ=255)

                          cGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGaa                     <  1:4053668/51‑1 (MQ=255)

                          cGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGaa                     <  1:444577/51‑1 (MQ=255)

                           ggTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAg                    >  1:4380264/1‑51 (MQ=255)

                            gTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAgg                   <  1:4133782/51‑1 (MQ=255)

                            gTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAgg                   >  1:3679556/1‑51 (MQ=255)

                            gTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAgg                   <  1:150104/51‑1 (MQ=255)

                             tttGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgt                  >  1:4538716/1‑51 (MQ=255)

                             tttGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgt                  >  1:282459/1‑51 (MQ=255)

                             tttGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgt                  <  1:3100525/51‑1 (MQ=255)

                              ttgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAgg                   >  1:2461043/1‑49 (MQ=255)

                              ttgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtg                 >  1:1012316/1‑51 (MQ=255)

                              ttgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtg                 >  1:5505339/1‑51 (MQ=255)

                               tgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtgt                >  1:5770460/1‑51 (MQ=255)

                               tgttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGgtgt                <  1:4967918/51‑1 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               >  1:835872/1‑51 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               >  1:7535306/1‑51 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               >  1:7796120/1‑51 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               <  1:10180605/51‑1 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               <  1:8935419/51‑1 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               >  1:6569388/1‑51 (MQ=255)

                                gttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGtgtg               <  1:3746908/51‑1 (MQ=255)

                                 ttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGc              <  1:9573281/51‑1 (MQ=255)

                                 ttgAGCATGGTCAGGCCGTTGATCCACATTTTGGCC‑GGGAGAAGGTGTGCa             >  1:6309129/1‑51 (MQ=255)

                                 ttgAGCATGGTAAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGc              >  1:6599012/1‑51 (MQ=255)

                                  tgAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCa             <  1:6107577/51‑1 (MQ=255)

                                   gAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAg            >  1:9187207/1‑51 (MQ=255)

                                   gAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAg            <  1:2840329/51‑1 (MQ=255)

                                   gAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAg            <  1:1256229/51‑1 (MQ=255)

                                   gAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAg            <  1:10451780/51‑1 (MQ=255)

                                    aGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGc           >  1:5735224/1‑51 (MQ=255)

                                      cATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCca         <  1:2195608/51‑1 (MQ=255)

                                      cATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCca         <  1:10632556/51‑1 (MQ=255)

                                      cATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCca         >  1:819332/1‑51 (MQ=255)

                                      cATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCca         <  1:2527861/51‑1 (MQ=255)

                                       aTGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcac        >  1:2808366/1‑51 (MQ=255)

                                       aTGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcac        >  1:1965432/1‑51 (MQ=255)

                                        tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcaca       <  1:880307/51‑1 (MQ=255)

                                        tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcaca       >  1:7637174/1‑51 (MQ=255)

                                        tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcaca       <  1:1383442/51‑1 (MQ=255)

                                        tGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcaca       >  1:3292892/1‑51 (MQ=255)

                                         ggTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCacac      >  1:10233663/1‑51 (MQ=255)

                                           tCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCAcacca    >  1:2530166/1‑51 (MQ=255)

                                           tCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCAcacca    <  1:10431639/51‑1 (MQ=255)

                                            cAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCAcaccac   >  1:1228552/1‑51 (MQ=255)

                                            cAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCAcaccac   <  1:849138/51‑1 (MQ=255)

                                             aGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGTAGCCACaccacc  >  1:8837518/1‑51 (MQ=255)

                                             aGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACaccacc  >  1:3275470/1‑51 (MQ=255)

                                             aGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACaccacc  >  1:4513987/1‑51 (MQ=255)

TTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGTCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACC  >  REL606/70819‑70914

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