Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	358,031	G→A	9.3%	*459* (TGA→TAA)	tdcG‑I →	L‑serine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	358,031	0	G	A	9.3%	112.9 / 2.7	43	*459* (TGA→TAA)	tdcG‑I	L‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (2/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GATGGCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGC  >  NC_002947/357891‑358178                                                                                                                                             |                                                                                                                                                     gATGGCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCt                                                                                                                                               >  5:140717/1‑148 (MQ=255)  aTGGCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGa                                                                                                                                             <  3:236974/149‑1 (MQ=255)    ggCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTAATTTTCTGAGc                                                                                                                                           >  7:181229/1‑149 (MQ=255)    ggCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTAATTTTCTGAGc                                                                                                                                           >  7:184762/1‑149 (MQ=255)       cTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGCCCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGGTGATTTTCTGAGC‑TGa                                                                                                                                       >  3:216533/1‑149 (MQ=255)          gcgcGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGAcc                                                                                                                                     <  6:140717/148‑1 (MQ=255)                             ttCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCggg                                                                                                                  >  5:240226/1‑148 (MQ=255)                             atCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCCCGCGGCGCCCTGGCGGTCAGCGCTATCGAGTGCTGAATTTCTGAGC‑TGACCTTCTCGGCCTCTTCGCggg                                                                                                                  >  3:221853/2‑148 (MQ=255)                               cATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGAcc                                                                                                                                     >  8:128492/1‑127 (MQ=255)                               cATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGc                                                                                                                 <  6:240948/147‑1 (MQ=255)                                    cccTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCTGCCTGGCTGTCAGCGCTATCGAGTGCTAATTTTCTGAGC‑TGACCTTCAGGGCCTCTTCGCGGGCAAGccc                                                                                                           <  8:181229/148‑1 (MQ=255)                                    cccTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCTGCCTGGCGGTCAGCGCTATCGAGTGCTAATTTTCTGAGC‑TGACCTTCAGGGCCTCTTCGCGGGCAAGccc                                                                                                           <  8:184762/148‑1 (MQ=255)                                     ccTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGAGCTTCTGGGCCTCTTCGCGGGCTTGCCCGc                                                                                                         >  7:94274/1‑149 (MQ=255)                                     ccTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGc                                                                                                         >  5:242981/1‑149 (MQ=255)                                      cTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCt                                                                                                        >  2:188362/1‑149 (MQ=255)                                      cTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATAATCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCt                                                                                                        >  1:224297/1‑149 (MQ=255)                                        cGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCGcc                                                                                                      >  1:143043/1‑149 (MQ=255)                                         gACCATCTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTcc                                                                                                      <  3:163567/148‑1 (MQ=255)                                         gACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGATCTTCTGGGCCTCTTCGCGGGCTTGCCCGcgaa                                                                                                      <  8:94274/148‑4 (MQ=255)                                          aCCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTa                                                                                                    >  4:86024/1‑149 (MQ=255)                                           ccAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTg                                                                                                                                <  6:176171/120‑1 (MQ=255)                                           ccAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTg                                                                                                                                >  5:176171/1‑120 (MQ=255)                                              ggTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTTGACCTTCTGGGCCTCTTCGCGGGCAAGccc                                                                                                           <  1:24924/139‑1 (MQ=255)                                              ggTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGccc                                                                                                           >  2:24924/1‑138 (MQ=255)                                              ggTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAgg                                                                                                <  3:86024/149‑1 (MQ=255)                                               gTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGCCCTTCTGGGCCTCTTCGCGGGCACGCCCGCGCCTAGAgg                                                                                                >  8:193696/1‑148 (MQ=255)                                               gTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAgg                                                                                                >  7:251831/1‑148 (MQ=255)                                                     cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                         <  6:240226/149‑1 (MQ=255)                                                              cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGTCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTg                                                                                <  7:208728/149‑1 (MQ=255)                                                                     ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGc                                                                         <  1:188362/149‑1 (MQ=37)                                                                     ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGc                                                                         <  3:71936/149‑1 (MQ=37)                                                                      cGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCg                                                                        >  6:78688/1‑149 (MQ=255)                                                                             gATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGCACTGTGCAATCCCTGGAGGAGCGGGCTTGc                                                                 >  1:210566/1‑149 (MQ=255)                                                                                    aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCcgcg                                                            >  1:961/1‑147 (MQ=255)                                                                                     cGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAag                                                         <  3:96015/149‑1 (MQ=255)                                                                                             aCAAAGATACCTCGCGCGGCGGCCGGCCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCt                                                  <  8:251831/148‑1 (MQ=255)                                                                                               aaatccaCCTCGCGCGGCGGCCTGGGGGGCTGGGCTATCGATTGCTGATTTTCTCTGC‑TCACCTTCGGGGCCTCTTCGCGGGAAAGCCCGCTCATACAGCAACTGTGCAATCCCTGTAGTAGGGGGCTTGCCCGCGAAGAGGCCCTTGa                                               <  4:232513/143‑1 (MQ=21)                                                                                                 agagCCCTCGGGCGGCGGCCTGGCGTTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTTGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                             <  2:210566/149‑1 (MQ=255)                                                                                                 agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGa                                               >  3:232513/1‑147 (MQ=255)                                                                                                     cccgcgcgcgGCGGCCGGGCGGTCAGCGCTATCGAGTCCTGATTGCCTGAGC‑TCCCCTTCTGGGCCGCTTCGCGGGCAAGCCCGCTCCTACAGGAACGGGGCAATCCCTGTAGCAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGa                                         <  7:193696/145‑1 (MQ=25)                                                                                                          gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCCCTCCTAAACGAACTGTGCAAACCCAGTAGGAGCCGGCATGCCCGCGAAGAGGCCCTTGAGGCTGATGCac                                    >  2:220513/1‑148 (MQ=38)                                                                                                                     tGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTCCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCATGAGGCTGTTGCACCTGCAtcagcc                         <  6:242981/149‑4 (MQ=25)                                                                                                                                            tGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGc  >  3:88717/1‑149 (MQ=255)                                                                                                                                             gATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTg                                                                  >  6:236840/1‑84 (MQ=255)                                                                                                                                             |                                                                                                                                                     GATGGCCTTGCGCGGCGATGGTGAGCACTTCATTTCCCTCGACCAGGTAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGC‑TGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGC  >  NC_002947/357891‑358178

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