Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	679,996	Δ1 bp	8.5%	coding (883/1278 nt)	PP_0582 →	thiolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	679,996	0	G	.	8.5%	146.0 / 6.8	47	coding (883/1278 nt)	PP_0582	thiolase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base . (2/2);  total (26/21)
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

TTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCA  >  NC_002947/679849‑680140                                                                                                                                                    |                                                                                                                                                  tttGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGc                                                                                                                                                  >  3:44937/1‑149 (MQ=255) tttGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGc                                                                                                                                                  >  1:197118/1‑149 (MQ=255)       cGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGa                                                                                                                                            >  1:112873/1‑149 (MQ=255)               gcagggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCa                                                                                                                                    >  7:179416/1‑149 (MQ=255)                cagggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCGGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCAc                                                                                                                                   >  1:175201/1‑149 (MQ=255)                   ggcaCGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACgg                                                                                                                                 <  3:24938/148‑1 (MQ=255)                           gACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGaa                                                                                                                        <  2:197118/148‑1 (MQ=255)                            aCGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCTGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAg                                                                                                                       <  2:175201/148‑1 (MQ=255)                              ggcggGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCT‑CGGTGGATTTTGTCACGGGG‑CGCGAAGGgtt                                                                                                                   >  8:297861/1‑149 (MQ=255)                              ggcggGTAATTCCACGCCGTTGACCGATCGTGCTTCGCGCGTGCAGTAGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGCTGCTTGACTATCTGGTCGATGGCGAAACCGCT‑CGGTGGAGTTTGTCACGGGG‑CGCGAAggggt                                                                                                                   >  8:297855/1‑147 (MQ=255)                                 gggTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGgttgt                                                                                                                 <  6:78238/149‑1 (MQ=255)                                    tAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGa                                                                                                              <  4:289426/149‑1 (MQ=255)                                      aTTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATg                                                                                                            <  4:44937/149‑1 (MQ=255)                                          cACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTAACGGGG‑CGCTAAGGGTTGTTGATGgcgc                                                                                                        >  5:305763/1‑149 (MQ=255)                                                 ttGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTa                                                                                                 >  3:155254/1‑149 (MQ=255)                                                         tGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCTGAGCAGGAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTAtgccgtg                                                                                          <  6:305763/148‑1 (MQ=255)                                                         tGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGATGTAtgccgtgc                                                                                         <  2:112873/149‑1 (MQ=255)                                                           gTGCTTAGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTAtgccgtgccg                                                                                       >  2:128736/1‑149 (MQ=255)                                                              cTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACAGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCgcgc                                                                                    >  5:184585/1‑149 (MQ=255)                                                               ttCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCt                                                                                   >  3:246322/1‑149 (MQ=255)                                                                      gTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGGCCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTgg                                                                             <  1:74233/149‑1 (MQ=255)                                                                       tGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCa                                                                           <  1:136768/149‑1 (MQ=255)                                                                         cTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTTTCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑AGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACg                                                                         >  6:255395/1‑149 (MQ=255)                                                                          tGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACg                                                                         >  2:140805/1‑148 (MQ=255)                                                                           gTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGttgttg                                                                                                               >  2:116530/1‑109 (MQ=255)                                                                           gTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGttgttg                                                                                                               <  1:116530/109‑1 (MQ=255)                                                                                gCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGc                                                                  <  1:128736/149‑1 (MQ=255)                                                                                             ggCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAg                                                      <  8:179416/148‑1 (MQ=255)                                                                                                gcagcagcaGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACt                                                  >  8:64251/1‑149 (MQ=255)                                                                                                gcagcagcaGGGGTTGCAGGTGATTGCCTATCTGGTCGAAGGCGAAACCGCTGCCGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACt                                                  >  1:242833/1‑149 (MQ=255)                                                                                                  agcagcagGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCACCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACttt                                                <  6:131356/149‑1 (MQ=255)                                                                                                  agcagcagGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCACCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACttt                                                <  6:136169/149‑1 (MQ=255)                                                                                                    cagcagGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGa                                              <  4:155254/149‑1 (MQ=255)                                                                                                         ggggTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGcctac                                          >  4:24814/1‑148 (MQ=255)                                                                                                            gTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACtactac                                       >  2:262622/1‑148 (MQ=255)                                                                                                                cAGGTGATTGACTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCTGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAAt                                  <  2:242833/149‑1 (MQ=255)                                                                                                                    tGCTTGCCTATCTGGTCGATGGCGAAACCGCT‑CGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGGGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCAct                             <  7:297861/149‑2 (MQ=255)                                                                                                                    ggCTTGCCTATCTGTTCGATGGCGAAACCGCT‑CGGTGGATTTTGTCACGGGG‑CGCGACCGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCTGGACTTTGACTACTACGAAATCCACg                             <  7:297855/148‑1 (MQ=255)                                                                                                                              tCTGGTCGATGGCGAAATCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTg                    >  3:45264/1‑149 (MQ=255)                                                                                                                                tGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCg                  >  1:173325/1‑149 (MQ=255)                                                                                                                                      aTGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAg            >  3:300097/1‑149 (MQ=255)                                                                                                                                      aTGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAg            >  3:774/1‑149 (MQ=255)                                                                                                                                      aTGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTACGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAg            >  6:11862/1‑149 (MQ=255)                                                                                                                                           gAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCt       >  5:199029/1‑149 (MQ=255)                                                                                                                                           gAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGAGCt       >  1:271693/1‑149 (MQ=255)                                                                                                                                             aatcGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCtgt     <  4:45264/146‑1 (MQ=255)                                                                                                                                                cGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCTCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCa  <  4:300097/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTTGACCGCAGCGGGCAGGGCACGTTGACGGCGGGTAATTCCACGCCGTTGACCGATGGTGCTTCGCTGGTGCTGTTGGGCAGCGAGGCGTGGGCGCAGCAGCAGGGGTTGCAGGTGCTTGCCTATCTGGTCGATGGCGAAACCGCTGCGGTGGATTTTGTCACGGGG‑CGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCA  >  NC_002947/679849‑680140

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