Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,496,381	T→G	8.9%	C270W (TGT→TGG)	mdeA →	methionine gamma‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,496,381	0	T	G	8.9%	162.6 / 5.4	54	C270W (TGT→TGG)	mdeA	methionine gamma‑lyase
Reads supporting (aligned to +/- strand):  ref base T (28/21);  new base G (2/3);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGGTTGGTGGTGGGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCTTTGA  >  NC_002947/1496234‑1496521                                                                                                                                                    |                                                                                                                                             cGGGTTGGTGGTGGGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATtgtg                                                                                                                                             <  1:367135/149‑1 (MQ=255)  gggTTGGTGGTGGGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGc                                                                                                                                            <  5:110577/149‑1 (MQ=255)    gTTGGTGGTGGGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGcc                                                                                                                                           >  6:191922/1‑148 (MQ=255)         tggtggGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGcc                                                                                                                                     <  2:48489/149‑1 (MQ=255)          ggtggGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGGCACAATGccg                                                                                                                                    <  8:94607/149‑2 (MQ=255)          gggggggCGCAAGCGGATGGGCGATCCCATTCGCCGCGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGGCACAATGccg                                                                                                                                    <  3:199229/146‑2 (MQ=255)           gtggGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGcc                                                                                                                                     >  6:137064/1‑147 (MQ=255)             ggggCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGGCACAATGccg                                                                                                                                    >  6:85053/1‑145 (MQ=255)             ggggCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGGCACAATGccg                                                                                                                                    <  5:85053/146‑2 (MQ=255)                    aaGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCcagca                                                                                                                                <  8:92040/143‑1 (MQ=255)                    aaGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCcagca                                                                                                                                >  7:92040/1‑143 (MQ=255)                          ttGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGtt                                                                                                                    >  5:137313/1‑149 (MQ=255)                            ggTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTtc                                                                                                                  >  7:247710/1‑149 (MQ=255)                             gTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTtc                                                                                                                  <  3:211469/148‑1 (MQ=255)                              tCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCTCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTtctc                                                                                                                <  7:252044/149‑1 (MQ=255)                              tCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTtctc                                                                                                                <  1:150517/149‑1 (MQ=255)                                gATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGt                                                                                                              >  1:19295/1‑149 (MQ=255)                                     cATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCgg                                                                                                          >  2:109478/1‑148 (MQ=255)                                       ttCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAg                                                                                                       >  5:153220/1‑149 (MQ=255)                                        tCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGc                                                                                                      >  2:173811/1‑149 (MQ=255)                                            cTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCg                                                                                                    >  8:283307/1‑147 (MQ=255)                                                aaGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTg                                                                                              <  7:283307/149‑1 (MQ=255)                                                    gCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAAt                                                                                          <  1:173811/149‑1 (MQ=255)                                                        aaGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGAt                                                                                      <  3:91773/149‑1 (MQ=255)                                                        aaGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCAGTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGAt                                                                                      >  7:160478/1‑149 (MQ=255)                                                          ggATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGAtt                                                                                     >  7:39814/1‑148 (MQ=255)                                                             tatGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCAtt                                                                                 <  6:137313/149‑1 (MQ=255)                                                                gACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGGCGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGTGTTTATCGTCCGGCAGCCGCAGGTGGAATTGATTCATTAcc                                                                              >  1:164675/1‑149 (MQ=255)                                                                   cGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCgg                                                                           >  7:24728/1‑149 (MQ=255)                                                                    ggggCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCggg                                                                          <  6:348380/149‑1 (MQ=255)                                                                             cTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTc                                                                 <  1:109478/149‑1 (MQ=255)                                                                               gTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGAAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCggatt                                                               >  4:29871/1‑145 (MQ=255)                                                                                 cGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCtttt                                                             <  5:191922/149‑1 (MQ=255)                                                                                  gCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTg                                                            >  4:76679/1‑149 (MQ=255)                                                                                    cccATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGcc                                                          >  4:276189/1‑149 (MQ=255)                                                                                    cccATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGATACAGTCTTTTGcc                                                          >  6:291563/1‑149 (MQ=255)                                                                                       aTGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGccc                                                         <  1:114318/147‑1 (MQ=255)                                                                                       aTGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGccc                                                         >  2:114318/1‑147 (MQ=255)                                                                                         gACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTc                                                                                    >  8:220105/1‑118 (MQ=255)                                                                                         gACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTc                                                                                    <  7:220105/118‑1 (MQ=255)                                                                                                     ctgATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCa                                         <  3:359699/149‑1 (MQ=255)                                                                                                      tgATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAg                                        >  3:204451/1‑149 (MQ=255)                                                                                                                gATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGc                              >  1:124570/1‑149 (MQ=255)                                                                                                                    aaGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGttt                          <  5:137064/149‑1 (MQ=255)                                                                                                                    aaGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGGCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGgtt                          >  8:275224/1‑149 (MQ=255)                                                                                                                     aGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTg                         >  4:196335/1‑149 (MQ=255)                                                                                                                                 tgcGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCggg                >  3:309811/1‑146 (MQ=255)                                                                                                                                 tgcGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCggg                <  4:309811/146‑1 (MQ=255)                                                                                                                                      aTGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGc        >  8:290001/1‑149 (MQ=255)                                                                                                                                        ggACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCt      >  3:269786/1‑149 (MQ=255)                                                                                                                                          aCCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCttt    <  2:111720/149‑1 (MQ=255)                                                                                                                                          aCCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTGTTGCCCAGTACGCGTTGGCGCCGCGACAGATGCGGTTGCCTGGCGGGCTGATTGCCttg    >  4:208718/1‑148 (MQ=255)                                                                                                                                           tgggCATTGGCACAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGACTGGCGGGATGATTGCCttt    >  8:288060/3‑148 (MQ=255)                                                                                                                                            cGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCTTTGa  <  3:76679/149‑1 (MQ=255)                                                                                                                                            cGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCTTTGa  <  8:39814/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             CGGGTTGGTGGTGGGGCGCAAGGTGTTGGTCGATCGCATTCGCCTGGAAGGGCTGAAGGATATGACCGGGGCAGTGCTGTCGCCCCATGACGCGGCGCTGCTGATGCGCGGGATCAAGACGCTTGCCTTGCGCATGGACCGGCATTGTGCCAATGCCCAGCAGGTGGCGGAGTTTCTCGTCCGGCAGCCGCAGGTGGAATTGATTCATTACCCGGGGCTACCGTCTTTTGCCCAGTACGCGTTGGCCCAGCGACAGATGCGTTTGCCTGGCGGGATGATTGCCTTTGA  >  NC_002947/1496234‑1496521

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