Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,922,351	G→T	13.3%	P712P (CCC→CCA)	gcvP‑II ←	glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,922,351	0	G	T	13.3%	80.4 / 3.8	30	P712P (CCC→CCA)	gcvP‑II	glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base T (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

GATGCTTGCGCTGCCCCAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGC  >  NC_002947/5922204‑5922496                                                                                                                                                    |                                                                                                                                                  gATGCTTGCGCTGCCCCAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTgggg                                                                                                                                                 >  2:471097/1‑150 (MQ=255)          gcTGCCCCAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGaa                                                                                                                                        <  1:274119/150‑1 (MQ=255)            tGCCCCAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACg                                                                                                                                      >  1:159160/1‑150 (MQ=255)              ccccAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGTGGAATGCAGAACGTc                                                                                                                                    >  1:185738/1‑150 (MQ=255)                ccAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCtt                                                                                                                                  >  1:139594/1‑150 (MQ=255)                ccAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTGGGCGCGCGAGCGGATACCGAGCGGCCCCATACCCGGCCAGCCGCCGCCGTTGGGAATGCAGAACGTCtt                                                                                                                                  >  1:434017/1‑150 (MQ=255)                            gcTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCAt                                                                                                                      <  2:190388/150‑1 (MQ=255)                            gcTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCAt                                                                                                                      <  1:43541/150‑1 (MQ=255)                              tGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGTGGAATGCAGAACGTCTTGTGCAGGTTCAtgt                                                                                                                    <  2:185738/150‑1 (MQ=255)                                aCCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCAtgtgt                                                                                                                  <  2:550322/150‑1 (MQ=255)                                aCCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGACCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCAtgtgt                                                                                                                  <  1:323933/150‑1 (MQ=255)                                 ccGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCACCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCAtgtgtg                                                                                                                 <  2:477557/150‑1 (MQ=255)                                  cgcgCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGATGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGa                                                                                                                >  2:368618/1‑150 (MQ=255)                                                ggTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCaa                                                                                                  >  1:223367/1‑150 (MQ=255)                                                        ccAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCg                                                                                          <  1:358750/150‑1 (MQ=255)                                                        ccAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCg                                                                                          <  2:223367/150‑1 (MQ=255)                                                            gcaccggcaccACCGGATGGCTGGCGACGAACGGCTTGATGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGgc                                                                                      <  1:368618/150‑1 (MQ=255)                                                            gcaccggcaccACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGgc                                                                                      <  2:123220/150‑1 (MQ=255)                                                                cggcaccACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGTGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCaa                                                                                                  >  2:521850/1‑134 (MQ=255)                                                                cggcaccACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGTGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCaa                                                                                                  <  1:521850/134‑1 (MQ=255)                                                                              ggctggcGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTggg                                                                    >  2:434089/1‑150 (MQ=255)                                                                                    cgacgaACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGa                                                              >  1:400212/1‑150 (MQ=255)                                                                                        gaACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGc                                                                   >  2:147596/1‑141 (MQ=255)                                                                                        gaACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGc                                                                   <  1:147596/141‑1 (MQ=255)                                                                                               ttGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCg                                                   >  2:227895/1‑150 (MQ=255)                                                                                               ttGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCg                                                   >  2:3551/1‑150 (MQ=255)                                                                                                          gcgGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCTGCCATGTAAAc                                        <  2:139594/150‑1 (MQ=255)                                                                                                                          cccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTgg                        <  1:471097/150‑1 (MQ=255)                                                                                                                           ccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGt                       >  1:392678/1‑150 (MQ=255)                                                                                                                              tACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGca                    <  2:203868/150‑1 (MQ=255)                                                                                                                                                cgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTcgc  <  1:227895/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GATGCTTGCGCTGCCCCAGGGCGCGGCGCTGACCGCGCTGTTGTTCGGGTCCAACCCAGGCACCGGCACCACCGGATGGCTGGCGACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGC  >  NC_002947/5922204‑5922496

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