BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,366,599	A→G	L178S (TTA→TCA)	yfcI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,366,599	0	A	G	99.0%	429.2 / ‑6.4	100	L178S (TTA→TCA)	yfcI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (47/52); total (47/53)
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

CCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAA  >  REL606/2366527‑2366671

ccTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACt                                                                       >  2:444647/1‑76 (MQ=255)

 cTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTc                                                                      >  2:16524/1‑76 (MQ=255)

  tGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTcc                                                                     <  2:840013/76‑1 (MQ=255)

  tGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTcc                                                                     <  2:1411905/76‑1 (MQ=255)

   gTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCa                                                                    <  2:1508344/76‑1 (MQ=255)

   gTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCa                                                                    <  1:1232546/76‑1 (MQ=255)

    tAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCaa                                                                   <  1:315619/76‑1 (MQ=255)

     aaCTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAAc                                                                  <  1:1796056/76‑1 (MQ=255)

     aaCTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAAc                                                                  >  2:837075/1‑76 (MQ=255)

     aaCTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAAc                                                                  >  2:84390/1‑76 (MQ=255)

      aCTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACa                                                                 <  2:1117710/76‑1 (MQ=255)

      aCTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACa                                                                 <  1:101404/76‑1 (MQ=255)

        tAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAgc                                                               >  1:1033817/1‑76 (MQ=255)

         agcagcGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAgcg                                                              >  2:1983835/1‑76 (MQ=255)

            agcGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCa                                                           <  1:255036/76‑1 (MQ=255)

              cGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCAtt                                                         <  1:1017281/76‑1 (MQ=255)

               gAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCAttt                                                        <  1:1615459/76‑1 (MQ=255)

               gAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCAttt                                                        >  2:1628784/1‑76 (MQ=255)

                aaaCAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCAtttt                                                       <  1:707356/76‑1 (MQ=255)

                 aaCAATTTGGTCGGCTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTg                                                      <  2:334215/76‑1 (MQ=255)

                  aCAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGc                                                     >  1:1743612/1‑76 (MQ=255)

                    aaTTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCgg                                                   <  1:766657/76‑1 (MQ=255)

                    aaTTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCgg                                                   >  1:400328/1‑76 (MQ=255)

                     aTTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGgt                                                  <  1:1320933/76‑1 (MQ=255)

                     aTTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGgt                                                  <  2:1698475/76‑1 (MQ=255)

                      tttGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGgtg                                                 <  2:630544/76‑1 (MQ=255)

                        tGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGtt                                               >  2:1826266/1‑76 (MQ=255)

                        tGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGtt                                               >  1:609006/1‑76 (MQ=255)

                         ggTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTg                                              <  2:900415/76‑1 (MQ=255)

                          gTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGc                                             <  2:1824083/76‑1 (MQ=255)

                          gTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGc                                             >  2:1970885/1‑76 (MQ=255)

                           tCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCa                                            >  1:536622/1‑76 (MQ=255)

                           tCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCa                                            >  1:655512/1‑76 (MQ=255)

                           tCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCa                                            >  2:618486/1‑76 (MQ=255)

                            cGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                           >  1:101838/1‑76 (MQ=255)

                             gACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATa                                          <  2:655512/76‑1 (MQ=255)

                              aCTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATaa                                         <  1:156443/76‑1 (MQ=255)

                               cTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAt                                        <  2:446827/76‑1 (MQ=255)

                               cTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAt                                        <  1:1757313/76‑1 (MQ=255)

                                tAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtc                                       >  2:680303/1‑76 (MQ=255)

                                tAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtc                                       >  2:182574/1‑76 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      >  1:1797435/1‑76 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      >  2:1216676/1‑76 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  2:1235959/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  1:1273702/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  1:1261210/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  1:739963/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      >  2:1598489/1‑76 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  1:1143237/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  2:1956455/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      <  1:431086/76‑1 (MQ=255)

                                 aaTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtct                                      >  2:708380/1‑76 (MQ=255)

                                  aTCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtctc                                     <  2:1726098/76‑1 (MQ=255)

                                   tCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCg                                    <  1:75416/76‑1 (MQ=255)

                                   tCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCg                                    <  1:79824/76‑1 (MQ=255)

                                    cccAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCg                                                    >  2:693170/1‑59 (MQ=255)

                                    cccAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCg                                                    <  1:693170/59‑1 (MQ=255)

                                    cccAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGt                                   <  2:1102557/76‑1 (MQ=255)

                                    cccAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGt                                   >  2:1352722/1‑76 (MQ=255)

                                      cAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCa                                 <  2:1938461/76‑1 (MQ=255)

                                      cAACAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCa                                 <  1:183302/76‑1 (MQ=255)

                                      cAACAGAGCGCGGTGACGAATATGGTGCTGAATTGACTCAAACAGCGCCATTTGGCGGTGTTGCAGAAGCTCGTCa                                 <  2:1938479/76‑1 (MQ=255)

                                       aaCAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCAt                                <  2:37602/76‑1 (MQ=255)

                                       aaCAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCAt                                >  2:1281498/1‑76 (MQ=255)

                                        aCAGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATc                               <  1:876239/76‑1 (MQ=255)

                                          aGATCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCg                             <  1:1235019/76‑1 (MQ=255)

                                            aTCGCGCTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGc                           >  1:341166/1‑76 (MQ=255)

                                               gcgcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcacc                        <  1:407807/76‑1 (MQ=255)

                                               gcgcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcaca                        <  2:892188/76‑2 (MQ=255)

                                                cgcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGGGTTGCATAATCTCGTCATCCGGcacca                       >  2:1575193/1‑76 (MQ=38)

                                                 gcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcaccac                      >  2:992628/1‑76 (MQ=255)

                                                 gcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcaccac                      >  2:1243713/1‑76 (MQ=255)

                                                 gcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcaccac                      >  1:1748817/1‑76 (MQ=255)

                                                 gcTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGcaccac                      >  2:890767/1‑76 (MQ=255)

                                                  tggACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     <  2:467851/74‑1 (MQ=255)

                                                  cTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     <  2:1519137/76‑1 (MQ=255)

                                                  cTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     <  2:588015/76‑1 (MQ=255)

                                                  cTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     >  1:691365/1‑76 (MQ=255)

                                                  cTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     >  1:1683864/1‑76 (MQ=255)

                                                  cTGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACg                     >  1:1317365/1‑76 (MQ=255)

                                                   tGACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACgg                    >  1:457694/1‑76 (MQ=255)

                                                    gACGAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGt                   >  1:828495/1‑76 (MQ=255)

                                                       gAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAat                <  1:1292813/76‑1 (MQ=255)

                                                       gAATATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAat                >  1:1792256/1‑76 (MQ=255)

                                                        aaTATGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAata               >  2:894240/1‑76 (MQ=255)

                                                          tatGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATc             >  1:605760/1‑76 (MQ=255)

                                                          tatGTTTCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATc             >  1:50089/1‑76 (MQ=255)

                                                              tttCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATccacc         <  1:101573/76‑1 (MQ=255)

                                                              tttCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATccacc         >  1:157548/1‑76 (MQ=255)

                                                              tttCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATccaac         <  1:89048/76‑3 (MQ=255)

                                                                tCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCaa       <  1:1592717/76‑1 (MQ=255)

                                                                tCTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCaa       >  1:1524862/1‑76 (MQ=255)

                                                                 cTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAAc      >  1:1835328/1‑76 (MQ=255)

                                                                 cTGAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAAc      >  2:1945306/1‑76 (MQ=255)

                                                                 cTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAAc      <  1:375521/76‑1 (MQ=255)

                                                                   gAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACgg    <  1:1586278/76‑1 (MQ=255)

                                                                   gAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACgg    >  1:792537/1‑76 (MQ=255)

                                                                   gAATTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACgg    >  1:792548/1‑76 (MQ=255)

                                                                    aaTTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGa   >  2:704873/1‑76 (MQ=255)

                                                                     aTTGACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGaa  <  1:84390/76‑1 (MQ=255)

CCTGTAACTAGCAGCGAAACAATTTGGTCGACTAATCCCAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATATCCACCAACGGAA  >  REL606/2366527‑2366671

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