BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,815,057	G→A	V397V (GTG→GTA)	katE →	catalase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,815,057	0	G	A	97.2%	112.9 / ‑4.0	36	V397V (GTG→GTA)	katE	catalase II
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (18/17); total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

CTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTT  >  NC_000913/1814926‑1815184

ctcGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGc                                                                                                                                <  1:543963/133‑1 (MQ=255)

        aaccaaACTTATCCCGGATGAAGTGGTGCCGGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGggactg                                                                                                                        <  1:203272/133‑1 (MQ=255)

          ccaaACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGACCAGGCGGCGTTCCATCGTGGGCATATCGTACCGggactgga                                                                                                                      >  2:203272/1‑133 (MQ=255)

              aCTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTc                                                                                                                  >  1:109404/1‑133 (MQ=255)

                   tCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGggactg                                                                                                                        <  2:616634/122‑1 (MQ=255)

                   tCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGggactg                                                                                                                        >  1:616634/1‑122 (MQ=255)

                   tCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCaa                                                                                                             >  1:568196/1‑133 (MQ=255)

                   tCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCaa                                                                                                             >  1:433385/1‑133 (MQ=255)

                    cccGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGggactggac                                                                                                                     >  2:418802/1‑124 (MQ=255)

                    cccGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGggactggac                                                                                                                     <  1:418802/124‑1 (MQ=255)

                        gaagaaCTGGTGCCCGTTCAGCGTATCGTCAATATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATc                                                                                                        <  2:275889/133‑1 (MQ=255)

                          agaaCTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCg                                                                                                      <  1:331624/133‑1 (MQ=255)

                                   gCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAg                                                                                             >  2:407153/1‑133 (MQ=255)

                                     ccGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAggg                                                                                           >  1:483983/1‑133 (MQ=255)

                                       gTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGAc                                                                                         <  2:68565/133‑1 (MQ=255)

                                        ttCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGAc                                                                                         <  1:150489/132‑1 (MQ=255)

                                        ttCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGAc                                                                                         >  2:150489/1‑132 (MQ=255)

                                            gCGTGTCGGCAAAATGGTGCTCAATCGCAATCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTg                                                                                    <  1:407153/133‑1 (MQ=255)

                                            gCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTg                                                                                    <  1:88272/133‑1 (MQ=255)

                                            gCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTg                                                                                    <  1:373880/133‑1 (MQ=255)

                                             cGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGt                                                                                   <  2:544659/133‑1 (MQ=255)

                                                 tCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTtctc                                                                               <  1:417893/133‑1 (MQ=255)

                                                         aTGGTGCTCAATCGCTACACGGATAACTTCCTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCg                                                                       <  1:202834/133‑1 (MQ=255)

                                                               cTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCg                                                                       <  2:349144/127‑1 (MQ=255)

                                                               cTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCg                                                                       >  1:349144/1‑127 (MQ=255)

                                                                       cAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAgt                                                         >  2:356555/1‑133 (MQ=255)

                                                                                tAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTtggt                                                <  1:356555/133‑1 (MQ=255)

                                                                                         tGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAAt                                       >  1:520974/1‑133 (MQ=255)

                                                                                                        ggcggcTTTCCATCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCg                        >  2:620116/1‑133 (MQ=255)

                                                                                                                   aTCCTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGGCTTGGTGGGCCGAACTACCATGAGAATCCGAGTAACCGTcc             >  2:283651/1‑133 (MQ=255)

                                                                                                                      cTGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGAc          <  2:483983/133‑1 (MQ=255)

                                                                                                                       tGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGAcc         >  2:478012/1‑133 (MQ=255)

                                                                                                                       tGGGCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGAcc         >  1:468147/1‑133 (MQ=255)

                                                                                                                         ggCATATCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTg       >  1:346305/1‑133 (MQ=255)

                                                                                                                            atatCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGccc    <  2:109404/133‑1 (MQ=255)

                                                                                                                              atCGTACCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCtt  >  2:13372/1‑133 (MQ=255)

CTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTATACCGATACACAAATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTT  >  NC_000913/1814926‑1815184

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