Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	983,925	G→A	10.2%	A99T (GCT→ACT)	ndk →	nucleoside diphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	983,925	0	G	A	10.2%	105.2 / 2.9	39	A99T (GCT→ACT)	ndk	nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (2/2);  total (23/16)
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.49e-01

TTCTACGCTGAGCACAGCGAGCGCGGCTTCTTCGGTGACCTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGCGAAAGCTTACGGGT  >  NC_002947/983778‑984072                                                                                                                                                    |                                                                                                                                                    ttCTACGCTGAGCACAGCGAGCGCGGCTTCTTCGGTGACCTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAgctgc                                                                                                                                                    >  3:125364/1‑149 (MQ=255)            tcacaGCGATCGCGGCTTATTCGGTTACCTGTTTGCCTTAATTACTTCCGGTCCGGTCGTTGTCCTGGTTCTGGAAGGCTAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACTTTAAAGAAGCTGCTGCCGGCAcc                                                                                                                                          <  7:47468/147‑1 (MQ=255)                    agcgCGGCTTCTTCGGTGACCTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCt                                                                                                                                 <  4:40405/149‑1 (MQ=255)                                  ggTGACCTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGt                                                                                                                    <  4:225006/148‑1 (MQ=255)                                        cTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGACCGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGtcgatcga                                                                                                             >  6:236149/1‑149 (MQ=255)                                                         ttCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACg                                                                                            >  5:54558/1‑149 (MQ=255)                                                               tccggtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCg                                                                                       >  6:58484/1‑148 (MQ=255)                                                                  ggtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTACTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                <  7:78005/136‑1 (MQ=255)                                                                  ggtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTACTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                <  7:82432/136‑1 (MQ=255)                                                                  ggtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTACTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                >  8:78005/1‑136 (MQ=255)                                                                  ggtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTACTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                >  8:82432/1‑136 (MQ=255)                                                                   gtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTc                                                                                  <  1:162282/149‑1 (MQ=255)                                                                   gtcGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTc                                                                                  >  5:59864/1‑149 (MQ=255)                                                                      gTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACAAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTACg                                                                                >  8:199915/1‑148 (MQ=255)                                                                       ttGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGaa                                                                              >  1:270603/1‑149 (MQ=255)                                                                          tCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                >  8:168076/1‑128 (MQ=255)                                                                          tCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGtt                                                                                                <  7:168076/128‑1 (MQ=255)                                                                            cAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAgctg                                                                          <  6:54558/148‑1 (MQ=255)                                                                            cAGGTTCTGGAAGGCGAGAACGCCATTCCCAATAAAAGTGAGCTGATGGGTGCTACCAACACTAAATAAGCTGCTGCCGGCAACATCCGCGCTTACTTCGCCGAGTCGATCGACGCCAACGCCGTTCAAGGTTCGGACTCCGAAgatgc                                                                         >  1:14829/1‑149 (MQ=255)                                                                               gTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAgctgctgc                                                                      >  8:83616/1‑149 (MQ=255)                                                                                    ggAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTc                                                                 >  8:283480/1‑149 (MQ=255)                                                                                     gAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTc                                                                 <  8:160071/148‑1 (MQ=255)                                                                                      aaGTCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGTTTCGGACTCCGAAGCTGCTGCTGCTcgc                                                               <  2:14829/149‑1 (MQ=255)                                                                                        ggCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTcgcg                                                              >  5:158867/1‑148 (MQ=255)                                                                                           gagaACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTcgcg                                                              <  3:196815/145‑1 (MQ=255)                                                                                           gagaACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTcgcg                                                              >  4:196815/1‑145 (MQ=255)                                                                                             gaACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTc                                                                                               >  4:299210/1‑110 (MQ=255)                                                                                             gaACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTc                                                                                               <  3:299210/110‑1 (MQ=255)                                                                                              aaCGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGc                                                       <  1:123565/149‑1 (MQ=255)                                                                                                   cATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTAct                                                  >  6:273929/1‑149 (MQ=255)                                                                                                   cATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTAct                                                  >  2:111724/1‑149 (MQ=255)                                                                                                       gCCAAGAACCGTGAGCTGATGGGTGCTAACAACCCTAAAGAAGCTGCTGCCCGCACCAAACGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCCTACACGCTTCGGACTCCGAAGCTGCTCCTGCTCGCGAAATCGCGAAattctt                                              >  7:264621/1‑149 (MQ=255)                                                                                                        ccAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTActtctt                                              <  5:273929/148‑1 (MQ=255)                                                                                                            gAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCgctg                                         >  5:149683/1‑149 (MQ=255)                                                                                                                        gATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTaa                             >  7:42110/1‑149 (MQ=255)                                                                                                                                    cAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAg                 >  6:258070/1‑149 (MQ=255)                                                                                                                                     aaCCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACTCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGc                <  7:199915/149‑1 (MQ=255)                                                                                                                                      aCCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAGGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGc                >  4:24846/1‑148 (MQ=255)                                                                                                                                           aaaGGAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTTAGCGagtgg          >  2:264334/1‑145 (MQ=255)                                                                                                                                               aagCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGCGAAAGCTTAc      <  6:59864/149‑1 (MQ=255)                                                                                                                                                   tgctgcCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGCGAAAGCTTACgggt  >  5:295132/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TTCTACGCTGAGCACAGCGAGCGCGGCTTCTTCGGTGACCTGGTTGCCTTCATGACTTCCGGTCCGGTCGTTGTCCAGGTTCTGGAAGGCGAGAACGCCATTGCCAAGAACCGTGAGCTGATGGGTGCTACCAACCCTAAAGAAGCTGCTGCCGGCACCATCCGCGCTGACTTCGCCGAGTCGATCGACGCCAACGCCGTTCACGGTTCGGACTCCGAAGCTGCTGCTGCTCGCGAAATCGCTTACTTCTTCGCTGCTACCGAGGTAACCACTCGCTAAGCGAAAGCTTACGGGT  >  NC_002947/983778‑984072

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