Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,240,914	A→G	9.3%	F262S (TTC→TCC)	PP_3718 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,240,914	0	A	G	9.3%	115.1 / 2.5	43	F262S (TTC→TCC)	PP_3718	aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (2/2);  total (22/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 = 9.58e-01

TCGGAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTC  >  NC_002947/4240767‑4241062                                                                                                                                                     |                                                                                                                                                      tCGGAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGaa                                                                                                                                                      >  2:276715/1‑149 (MQ=255)   ggAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACa                                                                                                                                                    <  1:98039/149‑1 (MQ=255)    gAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAg                                                                                                                                                   >  5:315442/1‑149 (MQ=255)    gAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAg                                                                                                                                                   <  7:2083/149‑1 (MQ=255)             cacGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTa                                                                                                                                          <  7:128816/149‑1 (MQ=255)              acGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATg                                                                                                                                               <  3:51437/143‑1 (MQ=255)              acGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATg                                                                                                                                               >  4:51437/1‑143 (MQ=255)              acGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAg                                                                                                                                         >  7:81945/1‑149 (MQ=255)              acGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAg                                                                                                                                         >  8:16803/1‑149 (MQ=255)                    ccAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGACGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCtcg                                                                                                                                   <  8:305037/149‑3 (MQ=255)                    ccAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTgg                                                                                                                                   <  7:16803/149‑1 (MQ=255)                     cAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGAAGTCTGGGGTGACGCCAAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTgg                                                                                                                                   <  8:252894/149‑1 (MQ=255)                        cGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGc‑a                                                                                                                                >  1:243255/1‑147 (MQ=255)                              gTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGcc                                                                                                                        <  2:243255/149‑1 (MQ=255)                                aCGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCca                                                                                                                       >  3:316303/1‑148 (MQ=255)                                  gCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACAt                                                                                                                    <  8:657/149‑1 (MQ=255)                                  gCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACAt                                                                                                                    <  8:81945/149‑1 (MQ=255)                                   ccGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCt                                                                                                                                     <  8:166812/132‑1 (MQ=255)                                   ccGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCt                                                                                                                                     >  7:166812/1‑132 (MQ=255)                                               cTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGc                                                                                                                                  >  1:305153/1‑123 (MQ=255)                                               cTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGc                                                                                                                                  <  2:305153/123‑1 (MQ=255)                                                    gCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCt                                                                                                   <  2:114048/148‑1 (MQ=255)                                                              tgatgTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATg                                                                                        <  6:315442/149‑1 (MQ=255)                                                                   tCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCa                                                                                    <  4:316303/148‑1 (MQ=255)                                                                       gggTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAg                                                                                                          <  1:221310/122‑1 (MQ=255)                                                                       gggTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAg                                                                                                          >  2:221310/1‑122 (MQ=255)                                                                       gggTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACcccg                                                                               >  1:252936/1‑149 (MQ=255)                                                                              gCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                        <  1:276715/149‑2 (MQ=255)                                                                              gCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGc                                                                        >  7:24704/1‑149 (MQ=255)                                                                                     cAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGa                                                                  >  3:72798/1‑149 (MQ=255)                                                                                                tGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcggctcggc                                                        >  7:222118/1‑148 (MQ=255)                                                                                                  gCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCg                                                       >  6:23508/1‑147 (MQ=255)                                                                                                     aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGGACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCgaagaa                                                  >  6:147169/1‑149 (MQ=255)                                                                                                     aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGGACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCgaaga                                                   >  6:142267/1‑148 (MQ=255)                                                                                                            tGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCGACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccgc                                            >  1:96521/1‑149 (MQ=255)                                                                                                             gCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCa                                          >  6:36163/1‑149 (MQ=255)                                                                                                              cccAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCACGATCGGCTCGGCGAAGAATGCCGCCAc                                         >  2:41755/1‑149 (MQ=255)                                                                                                               ccAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACg                                        >  3:168759/1‑149 (MQ=255)                                                                                                                          aaCGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCACGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCAcc                             <  1:41755/149‑1 (MQ=255)                                                                                                                              aGGTCACCACTTCGTCGGCGACGGACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGc                         <  5:147169/149‑1 (MQ=255)                                                                                                                              aGGTCACCACTTCGTCGGCGACGGACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGc                         <  5:142267/149‑1 (MQ=255)                                                                                                                                   accacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGagag                    >  7:228494/1‑149 (MQ=255)                                                                                                                                     cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGa                  >  3:50208/1‑149 (MQ=255)                                                                                                                                       cTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATc                <  8:24704/149‑1 (MQ=255)                                                                                                                                                  acgaACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGaa     >  8:131810/1‑149 (MQ=255)                                                                                                                                                     aaCAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTc  >  5:231721/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCGGAAAAAATACACGCGCCCAGCGGCAGGTACGCCGAGGTCAGGCCTTTGGCGGTGGTGATGATGTCTGGGGTGACGCC‑AAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGC‑ACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTC  >  NC_002947/4240767‑4241062

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: ‑