Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,013,836	G→A	40.8%	intergenic (+41/+69)	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,836	0	G	A	40.8%	16.9 / 24.0	27	intergenic (+41/+69)	xpt/PP_5266	xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base A (7/4);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

TCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCC  >  NC_002947/6013691‑6013976                                                                                                                                                   |                                                                                                                                              tCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTgcag                                                                                                                                            <  4:267211/149‑1 (MQ=255)   gTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGT‑AAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTgcagca                                                                                                                                          >  8:70710/1‑148 (MQ=255)      ccAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGcaca                                                                                                                                        <  2:196484/148‑1 (MQ=255)               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                                                                                                           <  1:56161/149‑1 (MQ=25)                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                           <  6:247999/149‑1 (MQ=25)                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                           >  2:56161/1‑149 (MQ=25)                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTACGGCACAGCGTTCG‑AGACCTGTGCTGCACCTGTgg                                                                                                           >  5:247999/1‑149 (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)                                                                                                                                                   |                                                                                                                                              TCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGG‑CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCG‑AACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCC  >  NC_002947/6013691‑6013976

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