BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,289,523	A→G	61.0%	E384E (GAA→GAG)	amn →	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,289,523	0	A	G	61.0%	112.9 / 133.4	197	E384E (GAA→GAG)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand): ref base A (51/26); new base G (19/101); total (70/127)
Fisher's exact test for biased strand distribution p-value = 1.02e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 197 total aligned reads displayed.

GAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAG  >  minE/1289461‑1289588

gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3156835/1‑69 (MQ=255)

gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3332996/1‑69 (MQ=255)

gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:640528/1‑69 (MQ=255)

gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3053933/1‑69 (MQ=255)

gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:1953187/1‑69 (MQ=255)

   gAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTaa                                                               >  1:1802926/1‑63 (MQ=255)

   gAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTAc                                                               >  1:26095/1‑64 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1012454/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:2268575/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1247531/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1321907/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1770244/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:3435301/69‑1 (MQ=255)

     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:2319304/69‑1 (MQ=255)

     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:1524640/69‑1 (MQ=255)

     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:3073125/69‑1 (MQ=255)

     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:354940/69‑1 (MQ=255)

     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:2959379/69‑1 (MQ=255)

     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:775260/69‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGTGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1737916/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCTGCTTCTGCACtt                                           <  1:3315229/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1908740/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2154727/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2017288/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2029621/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1636134/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1622166/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:365932/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2506208/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:5318/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:682749/66‑1 (MQ=255)

                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2628394/66‑1 (MQ=255)

                   tgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTTCACtt                                           <  1:1075346/66‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2981720/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2736650/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2102169/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:3222644/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:3231018/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:506752/68‑1 (MQ=255)

                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:711392/68‑1 (MQ=255)

                        actggtactgTGGTAACCACATATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:838903/69‑1 (MQ=255)

                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:1929461/69‑1 (MQ=255)

                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:2124452/69‑1 (MQ=255)

                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:497864/69‑1 (MQ=255)

                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:2871956/69‑1 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCGCTTCGTTTTAACTTAAGcc                           >  1:3146857/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACtt                                >  1:2336130/1‑64 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:943183/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2769584/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:759848/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2999080/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:302948/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2708359/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:1180815/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:1542916/1‑69 (MQ=255)

                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:157389/1‑69 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2551453/68‑1 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2520257/68‑1 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2452039/68‑1 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:3126775/68‑1 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:3250327/68‑1 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:2869034/1‑69 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:1456035/1‑69 (MQ=255)

                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:2386660/1‑69 (MQ=255)

                                         cacaGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTa                               >  1:1386825/1‑58 (MQ=255)

                                         cacaGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGt                    >  1:2816549/1‑69 (MQ=255)

                                               ttaCAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1841356/67‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCTTAGCAAt              <  1:336185/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2239282/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:10023/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:27777/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1282989/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1299632/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1375943/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:698149/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2903262/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2368235/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:438954/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1618785/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1879842/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:3127571/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:3406297/69‑1 (MQ=255)

                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2353981/69‑1 (MQ=255)

                                               tGACAGGAACTCGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2036692/69‑1 (MQ=255)

                                                gACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAtt             >  1:339913/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAtt             >  1:475643/1‑66 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1488633/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:3136985/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:903827/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:810327/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1850029/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1706085/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:521256/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:2752300/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:448008/1‑69 (MQ=255)

                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:3385041/1‑69 (MQ=255)

                                                          gggAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:2658477/1‑69 (MQ=255)

                                                          gggAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:1657117/1‑69 (MQ=255)

                                                          gggAATTACGTTACTCAGCTTCTGCACTACGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:2926844/1‑69 (MQ=255)

                                                            gAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAg  >  1:1609232/1‑68 (MQ=255)

GAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAG  >  minE/1289461‑1289588

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