Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,832	T→G	14.0%	D228A (GAC→GCC)	PP_3355 ←	beta‑ketothiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,798,832	0	T	G	14.0%	80.5 / 3.3	47	D228A (GAC→GCC)	PP_3355	beta‑ketothiolase
Reads supporting (aligned to +/- strand):  ref base T (16/24);  new base G (3/4);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAAGGCGGATGAGTGCAATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAG  >  NC_002947/3798698‑3798960                                                                                                                                       |                                                                                                                                 ccAAAGGCGGATGAGTGCAATCGGAGTAGCGCGTGACCCAGGCTGCCGCCGCGCCATCGACCACCGCGCCACTGTTGCCAGCCGTCGTCACGCCACCGGCGTTCCCCGGGCGCATTCGGGCCACTGCGGCCGCGTc                                                                                                                                 >  2:431327/1‑136 (MQ=255)    aaGGCGGATGAGTGCAATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGt                                                                                                                              <  1:180256/136‑1 (MQ=255)             gAGTGCAATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCgtgg                                                                                                                     >  1:173974/1‑136 (MQ=255)              tgTGCAATCGGAGTAGCGCGTGACCATGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCCGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCgtggg                                                                                                                    <  2:605472/135‑1 (MQ=255)              cgTGCAATCGGAGTAGCGCTTGACCAAGGCTGCCGCCGCGCCATCGCCCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCgtggg                                                                                                                    <  1:282968/135‑1 (MQ=255)                            aGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCt                                                                                                      <  2:334351/136‑1 (MQ=255)                              cgctTGACCACGGCGGCCGCGGGCCCAGCCCCCACCGCGCATCTGGTGCCCCCCGTCGGCACGCCCCCGGCGGGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGGGGCTGTCCTGGCTGa                                                                                                    <  2:427893/136‑1 (MQ=255)                                   gACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCAGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTTGGCCGTGGGTGGCTGGACTGGCTGCCctcc                                                                                               >  1:593288/1‑135 (MQ=255)                                     ccAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTGCTGGCTGACctcaca                                                                                             >  2:504089/1‑135 (MQ=255)                                      cAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTGGGTAGGCCGGGGGTGGCTGTGCTGGCGGACctcacta                                                                                            >  1:133038/1‑135 (MQ=255)                                      cAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGGGGGGAGGCCGGGGGGGGCTGTCCTGGCGGGcctccctc                                                                                            >  1:569270/1‑136 (MQ=255)                                         ggcggcCGCCGCGCCAGCGCCCCCCGCGCAACTGTGCCCAGCCGTCGGCCCGCCCCCGGCGGGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  1:556552/136‑1 (MQ=255)                                         ggcggcCGCCGCGCCAGCGCCCCCCCCGCAACTGTTCCCCGCCGTCTGCACGCCACCGGCGGGCACCGGGCGCAGTCTGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  2:609629/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCATCGACCCCCGCGCAACTGTTGCCAGCCGTCGGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  2:27824/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  1:154950/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  1:602894/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTACACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  2:362669/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCATCGACCACCGCGCAACTGTGGCCCGCCGTCTGCACGCCACGGGCGGGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  1:465213/136‑1 (MQ=255)                                         ggCTGCCGCCGCGCCAGCGACCACCGCGCAACTGTTGCCAGCCGGCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacg                                                                                         <  1:150494/136‑1 (MQ=255)                                           gtGCCGCCGCGCCATCGACCACCGCTCAACTGTTGCCAGCGGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacgcc                                                                                       <  2:72095/135‑1 (MQ=255)                                           cTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCCCCGGGCGCAGTCGGGCCAGTGCGGCCGCGGCGGGTGGCCGGGGGTGGCTGTGCTGGGTGCCCGCACTcacgcc                                                                                       >  1:478470/1‑136 (MQ=255)                                           cTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacgcc                                                                                       <  2:415241/136‑1 (MQ=255)                                            tGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTATGCCGTGGGGGGGTGTCCTGGGTGACCTCACTcacgccc                                                                                      >  2:128621/1‑135 (MQ=255)                                            tGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGGCCTCACTcacgccc                                                                                      >  1:485196/1‑135 (MQ=255)                                              ccgccgCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacgccacg                                                                                    <  1:348659/136‑1 (MQ=255)                                              ccgccgCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTcacgccacg                                                                                    <  1:271721/136‑1 (MQ=255)                                                gccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCCCCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCGGGCTGCCCGCACTcacgccacgcc                                                                                  >  2:348659/1‑135 (MQ=255)                                                 ccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCCCGCCCCCGGCGGGCCCCCGGCGGCGTCGGGCCAGTGCGGCCGCGGCGGGTGGGCGGGGGGGGGTGTGCTGGCGGGCCTCCCACACCCCCCCCgg                                                                                 >  1:354816/1‑136 (MQ=255)                                                 ccgcGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCGGGCCTCACTCACGCCACGCgg                                                                                 >  2:15858/1‑136 (MQ=255)                                                       cATCGCCCACCCCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgcaatgc                                                                           <  2:464895/136‑1 (MQ=255)                                                            accaccGCGCAACTGTTGCCAGCCGTCAGCACGCCCCCCGCGTGCACCGGGCGCAGTCCGGCCAGTTCGGCCGCGTGGGGAGGGCGGGGGAGGGGGTGCTGGCTGGCAACACTCACGCCCCCCGGCcatgcaaagc                                                                      >  2:424878/1‑136 (MQ=255)                                                                              cccgCCGGCGGCCCGCCACGGGGGTGCACCGGGCGCAGTCGGGCCAGGGCGGCCGCGTCGGTAGGCCGGGGGGGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCtt                                                    <  2:429552/133‑1 (MQ=255)                                                                              ccccccGTCTCCACCCCCCGGGGGTGCACGGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCtt                                                    <  1:139624/132‑1 (MQ=255)                                                                              ccAGCGGTCTCCCCGCCCCGGGGGTGCACCGGGCGCAGGCGGGCCAGTGCGGCCGGGGCGGTAGGCCGTGGGTGGCTGGCGTGGCTGACCTCACTCACGCCACGCGGAAATGCAATGCCGGGTGGCTGGCAGCCtt                                                    <  2:128956/136‑1 (MQ=255)                                                                                   cGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCCGCAAAGCAAAGCCGCGTGGCTGGCAGCCTTCTAcc                                               >  1:449158/1‑136 (MQ=255)                                                                                                    cGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGCCCTCACTCACGCCACCCGGCAATGCAAATCCGCGTGGCTGGCATCCGTCTAACTACAGTGGCTGCGCggg                              >  1:465565/1‑135 (MQ=255)                                                                                                                gcAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTTGCCGCCTTCTACCTCCAGTGCCTGCGCGGGGACGCTGCCgct                  >  1:176596/1‑134 (MQ=255)                                                                                                                gcAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGAt                  >  1:514516/1‑136 (MQ=255)                                                                                                                          ctgggcggcgGGGGCGGTAGGCGGGGGGGGGCGGGCCTGGCTCACCTCACTCCCGCCAGGCGGAAATGCATGGCCGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg        <  1:609057/132‑1 (MQ=255)                                                                                                                          cAGTGCGGCCGGGGCGGTAGGCGGTGGGTGGCTGCCCTGGCTGACCTCACTCCCGCCACGCGGCAATGCAATGCCGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg        <  2:546928/136‑1 (MQ=255)                                                                                                                          cAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg        <  1:235116/136‑1 (MQ=255)                                                                                                                          cAGTGCGGCCGCGTCGGTAGGCCGTGGGGGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg        <  1:261578/136‑1 (MQ=255)                                                                                                                          cAGGGCGGCCGCGGCGGTTGGCGGTGGGTGGCTGTCCTGGCTGACCTCACTCCGGCCACGCGCCAATCCAATGCCGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg        <  1:518371/136‑1 (MQ=255)                                                                                                                               gggCCGGGTGGGTTGGGGGGGGGTGGCTGTCGGGGCTGACCTACATCCGGCCAGGGTGAAATGAAATGCCGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCa   <  1:623667/135‑1 (MQ=255)                                                                                                                               gggCCGGGTCGGTAGGCGGGGGGGGGGGGTCGTGGCTCACCTCACTCCCCCCACGGGGAAATGCATGGCGGGGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCa   <  1:614706/135‑1 (MQ=255)                                                                                                                               cGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGCCGCTGACGATTTCCTCGTCGAAccc   >  2:9873/1‑135 (MQ=255)                                                                                                                                ggCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCAACCCGCAAAGCAAATCCCCCTGGGAGGCAGCCATAAACCTCCAAAGCCAGACCGCTGACCCAGACGATATACTCGTCCAACCac  >  2:307683/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                 CCAAAGGCGGATGAGTGCAATCGGAGTAGCGCGTGACCAAGGCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAG  >  NC_002947/3798698‑3798960

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