BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	709,001	A→G	9.4%	Y102H (TAT→CAT)	yceI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	709,001	0	A	G	9.4%	283.6 / 16.8	96	Y102H (TAT→CAT)	yceI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (31/56); new base G (7/2); total (38/58)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCTG  >  Exported/708931‑709071

aGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGata                                                                        >  1:4210891/1‑71 (MQ=255)

  aTCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtt                                                                      <  1:856362/71‑1 (MQ=255)

  aTCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtt                                                                      <  1:1424643/71‑1 (MQ=255)

   tCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     >  1:4747887/1‑71 (MQ=255)

   tCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     >  1:1386783/1‑71 (MQ=255)

            atatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTAtt                                                            <  1:812136/71‑1 (MQ=255)

            atatCCAGTTCGGCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:3623531/71‑1 (MQ=255)

             tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     <  1:1325690/61‑1 (MQ=255)

             tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:2690065/63‑1 (MQ=255)

             tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:3753033/63‑1 (MQ=255)

             tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:3214466/63‑1 (MQ=255)

             tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTg                                                           <  1:1668441/71‑1 (MQ=255)

               tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCg                                                                 >  1:1753245/1‑63 (MQ=255)

               tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:1820984/65‑1 (MQ=255)

               tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:3858419/65‑1 (MQ=255)

               tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:1024008/70‑1 (MQ=255)

               tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         <  1:3721345/71‑1 (MQ=255)

                ccAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:3140962/1‑70 (MQ=255)

                 cAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAgg                                                        <  1:2629709/70‑1 (MQ=255)

                  aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:2385814/71‑1 (MQ=255)

                  aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:4196178/71‑1 (MQ=255)

                  aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:4195083/71‑1 (MQ=255)

                  aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:3119284/71‑1 (MQ=255)

                   gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:2741610/64‑1 (MQ=255)

                   gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:323022/64‑1 (MQ=255)

                   gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:2058575/64‑1 (MQ=255)

                   gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaaa                                                     <  1:71023/71‑1 (MQ=255)

                      cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtttt                                                                    >  1:2710125/1‑53 (MQ=255)

                      cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:2919258/1‑64 (MQ=255)

                      cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:3136263/1‑64 (MQ=255)

                      cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:190927/1‑64 (MQ=255)

                      cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCt                                                  >  1:269371/1‑71 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 >  1:2969289/1‑71 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:537186/71‑1 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2665645/71‑1 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:4188940/71‑1 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:1695755/71‑1 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 >  1:2274470/1‑71 (MQ=255)

                       gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2154732/71‑1 (MQ=255)

                        tCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGc                                                <  1:594416/71‑1 (MQ=255)

                         cACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCa                                               >  1:3960237/1‑71 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:2257355/59‑1 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:1304137/59‑1 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:2070267/59‑1 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:1701364/1‑71 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:4747249/1‑71 (MQ=255)

                          aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:3566381/1‑71 (MQ=255)

                            cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:602545/52‑1 (MQ=255)

                              tctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3381919/71‑1 (MQ=255)

                              tctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3599000/71‑1 (MQ=255)

                               ctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3544634/70‑1 (MQ=255)

                                tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2768733/62‑1 (MQ=255)

                                tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:3852228/62‑1 (MQ=255)

                                tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:3863720/62‑1 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:650069/1‑49 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:4561198/1‑49 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:766149/1‑49 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTg                                            >  1:1344090/1‑65 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:877286/1‑67 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:2842289/1‑67 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:4643470/1‑67 (MQ=255)

                                  tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGaa                                        <  1:1540179/69‑1 (MQ=255)

                                    ttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAAt                                                    >  1:643488/1‑55 (MQ=255)

                                    ttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAAt                                                    >  1:3804013/1‑55 (MQ=255)

                                      cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATg                                    <  1:3591480/69‑1 (MQ=255)

                                       aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:963352/71‑1 (MQ=255)

                                       aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:2468774/71‑1 (MQ=255)

                                       aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:11483/71‑1 (MQ=255)

                                       aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:1793961/71‑1 (MQ=255)

                                          cTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               <  1:873697/70‑1 (MQ=255)

                                                 aGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:3627036/71‑1 (MQ=255)

                                                 aGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:1389041/71‑1 (MQ=255)

                                                  tgTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:881780/69‑1 (MQ=255)

                                                  ggTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:3021797/1‑71 (MQ=255)

                                                     gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:3441264/1‑59 (MQ=255)

                                                     gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:55913/1‑59 (MQ=255)

                                                     gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:2713551/1‑59 (MQ=255)

                                                     gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:2464745/1‑59 (MQ=255)

                                                     gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:1692551/1‑68 (MQ=255)

                                                     gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:4597172/1‑68 (MQ=255)

                                                     gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:2630899/1‑68 (MQ=255)

                                                      aaTGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:2010665/71‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:1122831/70‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:3878876/70‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:2640879/70‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:4458898/70‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   >  1:1448243/1‑70 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:999372/71‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:3731636/71‑1 (MQ=255)

                                                      aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:1366492/71‑1 (MQ=255)

                                                       aTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGa                 >  1:867164/1‑71 (MQ=255)

                                                        tGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGt                                   <  1:2044292/52‑1 (MQ=255)

                                                          ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGa                 >  1:2420128/1‑68 (MQ=255)

                                                          ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:2293915/1‑69 (MQ=255)

                                                          ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:951578/1‑69 (MQ=255)

                                                          ttGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:2280829/66‑1 (MQ=255)

                                                          ttGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:769369/1‑69 (MQ=255)

                                                                    atatTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGc    >  1:2360724/1‑71 (MQ=255)

                                                                      atTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCtg  >  1:3213359/1‑71 (MQ=255)

                                                                      atTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCtg  >  1:3914724/1‑71 (MQ=255)

AGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCTG  >  Exported/708931‑709071

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