BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,760,795	A→G	52.1%	A62A (GCA→GCG)	yrdC →	dsRNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,760,795	0	A	G	52.1%	9.5 / 108.8	123	A62A (GCA→GCG)	yrdC	dsRNA‑binding protein
Reads supporting (aligned to +/- strand): ref base A (28/31); new base G (39/25); total (67/56)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 124 total aligned reads displayed.

ACAGCAGTGATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGCAGCGTG  >  minE/2760730‑2760856

aCAGCAGTGTTGCGACTGTTGGAGTTCAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCgcagca                                                            <  1:3354897/69‑1 (MQ=255)

        gATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACg                                                     >  1:3195780/1‑68 (MQ=255)

           gCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1013786/68‑1 (MQ=255)

           gCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1910215/68‑1 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGGTGATTTTAATCGCGGCAAATTa                                                       >  1:1320051/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:1934005/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:1030398/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2275766/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:326893/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2479856/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:169633/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:80157/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:3150712/1‑61 (MQ=255)

             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2813039/1‑61 (MQ=255)

                tGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3314133/69‑1 (MQ=255)

                tGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:304916/69‑1 (MQ=255)

                    ggggTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTaaa                                          <  1:1571319/64‑1 (MQ=255)

                    ggAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1453672/59‑1 (MQ=255)

                      aGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:2136160/68‑1 (MQ=255)

                      aGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:3212645/68‑1 (MQ=255)

                       gTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCta                                     <  1:3208709/69‑1 (MQ=255)

                          aaaaCAGCGTCCGGTTGAGAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCtata                                   >  1:390986/1‑68 (MQ=255)

                            aaCAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATAtt                                 <  1:2610808/68‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2718633/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2530608/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2354809/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:978926/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:330632/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2159048/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:578063/69‑1 (MQ=255)

                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:1052111/69‑1 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:2783391/2‑68 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:206033/2‑68 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:1944204/2‑68 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:126409/2‑68 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:2990435/2‑68 (MQ=255)

                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:945391/2‑68 (MQ=255)

                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:1814937/69‑1 (MQ=255)

                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:455096/69‑1 (MQ=255)

                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:1631637/69‑1 (MQ=255)

                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:3133682/68‑1 (MQ=255)

                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:775690/68‑1 (MQ=255)

                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:133840/68‑1 (MQ=255)

                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:1192586/68‑1 (MQ=255)

                                          gATATGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtg                  >  1:1489191/1‑69 (MQ=255)

                                          gATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:44542/48‑1 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgatga                            >  1:1921638/1‑55 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgtctg              >  1:18554/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1175092/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1603630/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1599749/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2884431/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:36580/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:3369985/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2701093/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:909443/1‑69 (MQ=255)

                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2407560/1‑69 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:783531/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:1028413/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:398459/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:3408756/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3397746/37‑1 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:12298/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:2955407/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:1889122/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:969439/1‑37 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3079634/37‑1 (MQ=255)

                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:2427949/1‑37 (MQ=255)

                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:192822/47‑1 (MQ=255)

                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGc        >  1:2718308/1‑69 (MQ=255)

                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGc        >  1:411601/1‑69 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2363331/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2458767/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:1368092/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3215605/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:218035/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3042222/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2492263/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2087600/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3386591/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2880236/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2513726/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2573513/1‑68 (MQ=255)

                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:363324/1‑68 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:459331/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1022107/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:203824/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:3414437/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:3079589/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2851150/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2850958/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2838599/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1649086/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2793584/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1766700/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2663355/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2558602/68‑1 (MQ=255)

                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2054368/68‑1 (MQ=255)

                                                           aaTCGCGGCAAAATACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2397528/68‑1 (MQ=255)

                                                            aTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1892762/67‑1 (MQ=255)

ACAGCAGTGATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGCAGCGTG  >  minE/2760730‑2760856

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: `‑`