Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,013,852	A→G	39.3%	intergenic (+57/+53)	xpt → / ← PP_5266	xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,013,852	0	A	G	39.3%	18.7 / 25.2	28	intergenic (+57/+53)	xpt/PP_5266	xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base A (9/8);  new base G (7/4);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

CCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGC  >  NC_002947/6013705‑6013999                                                                                                                                                      |                                                                                                                                                    ccGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AAc                                                                                                                                                    <  1:300679/149‑1 (MQ=255)                cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGT‑AAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCtgt                                                                                                                                    <  7:70710/149‑1 (MQ=255)                 aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCtgt                                                                                                                                    <  2:24719/149‑1 (MQ=25)                  gCCATGGTTACCGCGTTTAATCGCTGGCCCGGGTGAAGTCGCTGGAATGTGGTGTGGGCAGCTTCAACGAATAAACGCGACACCTGGTTCGCGGGCCACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCtgt                                                                                                                                    <  1:248300/149‑1 (MQ=18)                   ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                                                  >  5:247999/1‑149 (MQ=25)                   ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                                                  <  6:247999/149‑1 (MQ=25)                   ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                                                  >  2:56161/1‑149 (MQ=25)                   ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                                                  <  1:56161/149‑1 (MQ=25)                    cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTggg                                                                                                                                 <  1:129941/149‑1 (MQ=255)                     aGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGa                                                                                                                                <  7:64593/149‑1 (MQ=255)                               cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGc                                                                                                                      <  6:326848/149‑1 (MQ=255)                               cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACCCCCCCTACCACAGGTGCAGCACAGGTCTCG‑CACCCTGTGCCCTACCATTGGGAGCGGGccggc                                                                                                                      >  6:114727/1‑145 (MQ=255)                                   gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCCGG‑CACGCCCGCTCCCACAGGTGCACCACAGGTCTCG‑AACGCTGTGCCCAACCTGTGGGCGCGGGCGTGccccc                                                                                                                  >  4:87548/1‑146 (MQ=39)                                                                           tCAGCTTCCTCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATgg                                                                                                            <  1:215986/115‑1 (MQ=255)                                                                                            aCGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCCACAACAAGCGCTg                                                          >  1:252148/1‑149 (MQ=255)                                                                                                        cTTCGCGGG‑CACGCCCGCTCCCACTGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttctt                                               <  2:252148/147‑1 (MQ=255)                                                                                                        attCGCGGG‑CACGCCCGCTCCCACAGGGACGGCACAGCGTTCG‑AGCCCTGTGCTGCCCCTGGGGGGGCGGGCGTGCCCGCGAATGGCCGCAACCCTTGTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGGTAAGttagtct                                             >  5:192731/2‑144 (MQ=17)                                                                                                              ggg‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGTGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGc                                       <  5:202920/149‑1 (MQ=255)                                                                                                                     gCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                 >  1:86900/1‑149 (MQ=255)                                                                                                                      cccGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                >  2:179799/1‑149 (MQ=255)                                                                                                                       ccGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg                               >  8:274443/1‑149 (MQ=255)                                                                                                                         gCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt                             >  4:302142/1‑149 (MQ=21)                                                                                                                             ccACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCACCCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                         >  7:254700/1‑149 (MQ=21)                                                                                                                             ccACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                         >  8:225436/1‑149 (MQ=21)                                                                                                                              cacaGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                         >  8:287411/1‑148 (MQ=255)                                                                                                                              cacaGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                         >  5:21460/1‑148 (MQ=21)                                                                                                                                        gcaCAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTGTTCTTTAAGACCCTGCGGCTGACCCACCTGCATcc              >  8:229736/1‑149 (MQ=39)                                                                                                                                               tctcG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCTCCAACTGCATCCGCCCCAg       >  5:205497/1‑149 (MQ=255)                                                                                                                                                     aaCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGc  <  2:86900/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGC  >  NC_002947/6013705‑6013999

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