Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,060,386	G→T	12.6%	T171T (ACC→ACA)	pgm ←	phosphoglucomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,060,386	0	G	T	12.6%	107.0 / 4.6	39	T171T (ACC→ACA)	pgm	phosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base T (3/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACGTTTTCCAGGTCGGCCACATAGTTGCTGACGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATGGCGTGA  >  NC_002947/4060239‑4060527                                                                                                                                                    |                                                                                                                                              gATGACGTTTTCCAGGTCGGCCACATAGTTGCTGACGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTg                                                                                                                                               <  6:28995/148‑1 (MQ=255)   tGACGTTTTCCAGGTCGGCCACATAGTTGCTGACGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTc                                                                                                                                            <  5:232781/149‑1 (MQ=255)                                 tGACGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTtgggt                                                                                                              >  8:224390/1‑149 (MQ=255)                                    cGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTttggtggg                                                                                                           <  7:224390/149‑1 (MQ=255)                                               cGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGaa                                                                                                <  6:107314/149‑1 (MQ=255)                                                    ggtggtCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTc                                                                                                                                       <  6:137715/105‑1 (MQ=255)                                                         tCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCtt                                                                                      >  8:121477/1‑149 (MQ=255)                                                         tCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCt                                                                                       >  5:24967/1‑148 (MQ=255)                                                         tCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCt                                                                                       <  6:24967/148‑1 (MQ=255)                                                                cTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTgcggcg                                                                               >  2:251924/1‑149 (MQ=255)                                                                       gcCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTtgtg                                                                        >  6:75265/1‑149 (MQ=255)                                                                        ccTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGc                                                                       >  1:316501/1‑149 (MQ=255)                                                                                    ccATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTtgt                                                                         <  4:147589/135‑1 (MQ=255)                                                                                     cATGCGCTTGACGCCCGCCAGTTTAGCCGCCAGCATCTCGTTGGCCTTGCCCTCGATCCACTTTGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGcggcgg                                                                              <  3:271249/129‑1 (MQ=255)                                                                                     cATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTTGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGcggcgg                                                                              <  3:266506/129‑1 (MQ=255)                                                                                     cATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTTGTCACGTCGCTGTCGGCCGCGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGcggcgg                                                                              >  4:266506/1‑129 (MQ=255)                                                                                     cATGCGCTTGACGCCCGCCAGATTAGACGCCAGCAGCTCGTTGGCCTTGCCCTCCATCCACAATGTCACGTAGCTGTCGGACGGCCCGCCGTGGGGTGGCTTGTACTTGAAGCAACCGCTTTGcggcgg                                                                              >  4:271249/1‑129 (MQ=255)                                                                                      aTGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACGGGTTTGGGGGGGGTTGTGGGTGGGGGtgtttgt                                                         >  2:138555/1‑144 (MQ=255)                                                                                      aTGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGtga                                                             <  1:138555/145‑1 (MQ=255)                                                                                      aTGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGAc                                                         <  1:251924/149‑1 (MQ=255)                                                                                         cgcTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGa                                                       <  2:316501/148‑1 (MQ=255)                                                                                              gACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGtgat                                                            <  3:80348/138‑1 (MQ=255)                                                                                              gACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGtgat                                                            >  4:80348/1‑138 (MQ=255)                                                                                                    cgccAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCa                                           >  2:128806/1‑149 (MQ=255)                                                                                                             aGCCGCCAGCAGCTCGTTGGCCTTGCCCTCAATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCtgc                                   <  4:169380/148‑1 (MQ=255)                                                                                                               ccgccAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCtgcgtg                                >  1:153234/1‑149 (MQ=255)                                                                                                                 gccAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTc                                                <  3:60628/131‑1 (MQ=255)                                                                                                                        gcTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCgg                       >  7:39441/1‑149 (MQ=255)                                                                                                                          tCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCgg                       >  4:200890/1‑147 (MQ=255)                                                                                                                          tCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCgg                       <  3:200890/147‑1 (MQ=255)                                                                                                                            gTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGGCCTGCGTGCGGCCGCGGtt                     <  8:39441/147‑1 (MQ=255)                                                                                                                             ttGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTtgtg                  <  1:138063/149‑1 (MQ=255)                                                                                                                               ggCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTgg                 <  2:153234/148‑1 (MQ=255)                                                                                                                                  cTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCa               >  2:71792/1‑147 (MQ=255)                                                                                                                                     gCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCagag            <  7:121477/147‑1 (MQ=255)                                                                                                                                        cTCGATCCACTTTGTCACGTCGCTGTCGGACGCGCCGCATTTGGGTGCGTTGTACTTGAAGCCAACGCTTTGCGGCGGGTTGTGCTACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGACGCGGTtgt                   >  5:236188/1‑137 (MQ=255)                                                                                                                                        cTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATgg       >  3:306068/1‑149 (MQ=255)                                                                                                                                        cTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATgg       >  8:18591/1‑149 (MQ=255)                                                                                                                                           gATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATGGCg     <  7:18591/148‑1 (MQ=255)                                                                                                                                             tCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATGGCGTga  <  1:71792/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GATGACGTTTTCCAGGTCGGCCACATAGTTGCTGACGTAGTCGTGGCGGTGGGTGGTCGGCGCCTGCAGCGCCTGGGCATGGTCCATGCGCTTGACGCCCGCCAGATTAGCCGCCAGCAGCTCGTTGGCCTTGCCCTCGATCCACTTGGTCACGTCGCTGTCGGCCGGGCCGCCGTTGGGTGGGTTGTACTTGAAGCCACCGCTTTGCGGCGGGTTGTGCGACGGGGTGATGACGATACCGTCAGCCAGGCCCTGCGTGCGGCCGCGGTTGTGGCAGAGGATGGCGTGA  >  NC_002947/4060239‑4060527

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