BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,423,342	C→A	56.4%	V177V (GTG→GTT)	ccmC ←	heme exporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,423,342	0	C	A	56.4%	7.8 / 56.8	55	V177V (GTG→GTT)	ccmC	heme exporter subunit
Reads supporting (aligned to +/- strand): ref base C (14/10); new base A (13/18); total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTG  >  minE/1423283‑1423397

gCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:407240/1‑68 (MQ=255)

gCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:2120577/1‑68 (MQ=255)

gCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:3156405/1‑68 (MQ=255)

gCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:3238636/1‑68 (MQ=255)

gCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:1793981/1‑68 (MQ=255)

gCCGGATAGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTa                                                 >  1:953171/1‑68 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:2459301/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:1989915/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:880329/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:740425/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:3128184/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:2240610/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:1714939/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:2831847/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:1473992/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:2602935/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:1330377/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaata                                          >  1:2829383/1‑69 (MQ=255)

      tCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCAACGGAGTAATGaat                                           >  1:1456870/1‑68 (MQ=255)

                   gcATCCGCGTTGATCCCTGATGCAGGGTGTTCGACCACTCCACGGAGTAATGAATAATCGGCAGATTca                             <  1:407283/69‑1 (MQ=255)

                   gcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTca                             <  1:1156629/69‑1 (MQ=255)

                   gcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTca                             <  1:2010297/69‑1 (MQ=255)

                        cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                         >  1:3270782/1‑68 (MQ=255)

                        cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                         >  1:1289059/1‑68 (MQ=255)

                        cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                         >  1:1263701/1‑68 (MQ=255)

                        cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                         >  1:1077718/1‑68 (MQ=255)

                        cgcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTcaccac                         >  1:2700073/1‑68 (MQ=255)

                           gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCa                     <  1:1610262/69‑1 (MQ=255)

                           gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCa                     <  1:518087/69‑1 (MQ=255)

                           gTTGAT‑CCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCa                     <  1:1515854/68‑1 (MQ=255)

                                        gCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCa                     <  1:3010142/56‑1 (MQ=255)

                                              tgtTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTg  >  1:190478/1‑69 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:3141003/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:3143318/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:3299741/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:700562/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:279785/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:2688211/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:2473472/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1004177/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1914529/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1892922/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1870847/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1850523/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1814376/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1805495/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1668577/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:1238531/53‑1 (MQ=255)

                                                      cacTCAACGGAGTAATGAATAATCGCCAGATTCACCACGCCAATCAGCACCAg          <  1:482962/53‑1 (MQ=39)

                                                       acTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                >  1:282111/1‑46 (MQ=255)

                                                       acTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                <  1:1117461/46‑1 (MQ=255)

                                                       acTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                <  1:1410219/46‑1 (MQ=255)

                                                       acTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                <  1:2858684/46‑1 (MQ=255)

                                                       acTCCACGGAGTAATGAATAATCGGCAGATTCACCACG‑CAATCAg                >  1:1021373/1‑45 (MQ=255)

                                                       acTCAACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg          <  1:2022908/52‑1 (MQ=255)

                                                           cACGGAGTAATGAATAATCGGCAGATTCCCCACGCCAATCAGCACCAgg         >  1:1755259/1‑49 (MQ=255)

                                                           cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGa        >  1:193415/1‑50 (MQ=255)

GCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTG  >  minE/1423283‑1423397

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