BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,809,486	A→C	T455P (ACT→CCT)	katE →	catalase HPII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,809,486	0	A	C	96.3%	90.4 / ‑6.0	27	T455P (ACT→CCT)	katE	catalase HPII
Reads supporting (aligned to +/- strand): ref base A (0/0); major base C (11/15); minor base T (1/0); total (12/15)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGA  >  NZ_CP009273/1809353‑1809592

aCGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAAc                                                                                                       >  1:712994/1‑139 (MQ=255)

         ctcCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAAt                                                                                              >  2:89439/1‑139 (MQ=255)

             tataCCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTAc                                                                                           >  2:321272/1‑138 (MQ=255)

              ataCCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTAcg                                                                                          >  1:675686/1‑138 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCTCTTACCCTAATTTCCAGCGTGACCGCATGCATCGCATGGGGATCGACTCTAACCCGGCGAATTAcgtaccgaa                                                                                   >  2:554479/1‑139 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGc                                                                                                  <  2:748653/124‑1 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGc                                                                                                  >  1:748653/1‑124 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTAcgaaccgaa                                                                                   >  2:664304/1‑139 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTAcgaaccgaa                                                                                   >  1:760699/1‑139 (MQ=255)

                    aTACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTAcgaaccgaa                                                                                   >  1:750715/1‑139 (MQ=255)

                      acacAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACt                                                                                 >  1:685970/1‑139 (MQ=255)

                       cacaAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACt                                                                                 <  2:492274/138‑1 (MQ=255)

                         caAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGa                                                                              <  2:90819/139‑1 (MQ=255)

                            aTCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTACCCCGGCGAATTACGAACCGAACTCGAtt                                                                            >  2:545894/1‑138 (MQ=255)

                                             gCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCcgcgc                                                          <  2:750715/139‑1 (MQ=255)

                                                       cATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACAccgc                                                <  1:384561/139‑1 (MQ=255)

                                                                                ccTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa                       <  2:608242/139‑1 (MQ=255)

                                                                                    cccTTGCCATAATTACCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGGGAATTGCGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  1:554479/137‑1 (MQ=255)

                                                                                    cccTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  1:560197/137‑1 (MQ=255)

                                                                                    cccTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  2:47300/137‑1 (MQ=255)

                                                                                    cccTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  1:176201/137‑1 (MQ=255)

                                                                                    cccTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  1:115874/137‑1 (MQ=255)

                                                                                             tAATTTCCAGCATGACGGCATGCATCGCATGGGGATCGTCCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACCCCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAg          <  1:716662/139‑1 (MQ=255)

                                                                                                  tCCAGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGcgcg      <  1:152455/138‑1 (MQ=255)

                                                                                                     aGCGTGACGGCATGCATCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGa  <  2:605204/139‑1 (MQ=255)

                                                                                                                 tgcatCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     <  2:119532/108‑1 (MQ=255)

                                                                                                                 tgcatCGCATGGGGATCGACCCTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                     >  1:119532/1‑108 (MQ=255)

ACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGA  >  NZ_CP009273/1809353‑1809592

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