BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,591,559	A→G	61.2%	A54A (GCA→GCG)	hdeB →	acid‑resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,591,559	0	A	G	61.2%	36.4 / 45.8	72	A54A (GCA→GCG)	hdeB	acid‑resistance protein
Reads supporting (aligned to +/- strand): ref base A (16/12); new base G (32/12); total (48/24)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GAATCCGCTAAAGATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGAA  >  minE/2591491‑2591619

gAATCCGCTAAAGATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCa                                                              >  1:1875512/1‑69 (MQ=255)

 aaTCCGCTAAAGATATGACCTGCCAGGAATTTATTGATCTGAACCCAAAAGCAATGACCCCGGTTGCa                                                              >  1:401851/1‑68 (MQ=255)

       cTAAAGATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGa                                                       <  1:1954301/69‑1 (MQ=255)

            gATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGAtgctg                                                  <  1:3253997/69‑1 (MQ=255)

            gATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGAtgctg                                                  <  1:1981225/69‑1 (MQ=255)

              tatGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGAtgctgc                                                 >  1:1116392/1‑68 (MQ=255)

              tatGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGAtgctgc                                                 >  1:2970085/1‑68 (MQ=255)

               atGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCa                                                              <  1:1379866/54‑1 (MQ=255)

                tGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGGTGCATGGTGGATGCTGCATg                                              <  1:1522195/69‑1 (MQ=255)

                    cTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATgaaga                                          >  1:3219275/1‑69 (MQ=255)

                    cTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATgaaga                                          >  1:2337788/1‑69 (MQ=255)

                    cTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATgaaga                                          >  1:2511546/1‑69 (MQ=255)

                        cAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGAtgctg                                                  <  1:2737703/57‑1 (MQ=255)

                          ggAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGt                                    <  1:2514785/69‑1 (MQ=255)

                             aTTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGtat                                  >  1:309353/1‑68 (MQ=255)

                             aTTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGtat                                  >  1:1128867/1‑68 (MQ=255)

                             aTTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGta                                   >  1:3169053/1‑67 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGCCCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:998995/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGtatata                               >  1:1042107/1‑68 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:1596527/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:811720/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:2192688/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:205086/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:3262756/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaa                              >  1:895921/1‑69 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:3344762/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:1263590/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:3384322/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:1769718/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:222171/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:1484541/69‑1 (MQ=255)

                                tATTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:2817414/69‑1 (MQ=255)

                                 aTTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:529735/68‑1 (MQ=255)

                                 aTTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:703579/68‑1 (MQ=255)

                                 aTTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:841863/68‑1 (MQ=255)

                                 aTTGATCTGAATCCAAAAGCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:3187504/68‑1 (MQ=255)

                                  ttGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaaa                             <  1:134944/68‑1 (MQ=255)

                                  ttGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATaaa                             <  1:3146391/68‑1 (MQ=255)

                                            tCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGc                       <  1:2739413/64‑1 (MQ=255)

                                            tCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGc                       <  1:2387040/64‑1 (MQ=255)

                                                 gaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3144527/2‑41 (MQ=25)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3133940/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3107325/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1431768/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3239106/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:131199/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2696069/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3292566/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:365905/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2331697/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1666618/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1803127/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1835813/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:188718/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1897718/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1936990/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1615050/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2241121/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:3091388/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2368548/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:53431/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2854973/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:2928584/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCTGCATgaagaa                                         >  1:1458018/1‑41 (MQ=37)

                                                 aaGCAATGACACCAGTTGCGTGGTGGATGCT‑CATgaagaa                                         >  1:1535720/1‑40 (MQ=19)

                                                   gCAATGACACCAGTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATaa                              <  1:1765783/50‑1 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGaa  >  1:117666/4‑69 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGaa  >  1:1464052/4‑69 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGaa  >  1:2589049/4‑69 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGaa  >  1:780738/4‑69 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGaa  >  1:1016765/4‑69 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGa   >  1:1501019/4‑68 (MQ=255)

                                                            ccagTTGCGTGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGa   >  1:1655957/4‑68 (MQ=255)

                                                            ccGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGc                       <  1:3338598/48‑1 (MQ=255)

                                                            ccGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGc                       <  1:2055842/48‑1 (MQ=255)

GAATCCGCTAAAGATATGACCTGCCAGGAATTTATTGATCTGAATCCAAAAGCAATGACCCCGGTTGCATGGTGGATGCTGCATGAAGAAACAGTATATAAAGGTGGCGATACCGTTACTTTAAATGAA  >  minE/2591491‑2591619

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