Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,013,838	A→G	9.9%	intergenic (+43/+67)	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,838	0	A	G	9.9%	106.5 / 3.3	42	intergenic (+43/+67)	xpt/PP_5266	xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base G (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

TCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC  >  NC_002947/6013691‑6013977                                                                                                                                                    |                                                                                                                                            tCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcag                                                                                                                                            <  7:220706/149‑1 (MQ=255) tCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACACGTgcag                                                                                                                                            >  4:317302/1‑149 (MQ=255)         gggCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtc                                                                                                                                    <  7:302947/149‑1 (MQ=255)                  tGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCt                                                                                                                           <  3:245814/149‑1 (MQ=255)                     cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                        <  7:310983/149‑1 (MQ=255)                     cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGCACGCtgtg                                                                                                                        >  7:62173/1‑149 (MQ=255)                               aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                              <  1:193308/149‑1 (MQ=25)                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTTCCGTACCTGTgg                                                                                                            >  6:108235/1‑149 (MQ=255)                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCACACACGTGCAGCACACGTCTCGAACGCTGTGCCGTAACAGTgg                                                                                                            >  7:287972/1‑149 (MQ=255)                                 ccAGGGGTACCGCGTAGAAGCGCTGGCCCGGGGGAAGTCGCTGGAAGGTGGTGGGGTCAGCGTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                            <  6:111966/149‑1 (MQ=18)                                   aGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                             >  6:243498/1‑146 (MQ=255)                                   aGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                             <  5:243498/146‑1 (MQ=255)                                      gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCCCGCCCCCCCGGGCACCACAGGGGGGGAACGCGGTGGCGGAGCGGTGGGTGCg                                                                                                       >  4:147533/1‑149 (MQ=255)                                            gcgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCGCTACCACGGGTGCAGCAAAGGTCTCGAACGCTGTGCCGTACCGGTGGGAGCGGGCGTg                                                                                                 >  7:195921/1‑149 (MQ=255)                                                           cccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGc                                                                                  <  7:36418/149‑1 (MQ=255)                                                            ccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCAAGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                        >  4:117518/1‑110 (MQ=255)                                                            ccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCAAGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                        <  3:117518/110‑1 (MQ=255)                                                                      tCGCTGGAAGGTGGTGTGATCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCCCAGGGGCCGCAGAGGGCTGGGACCCGGTGCCCTACCTGTGGGGGCGGGCGTGCCCGCGATAGGGCGCAACCCGTTTAt                                                                       >  4:81417/1‑149 (MQ=25)                                                                        gCTGGATGGTGGTGTGATCAGCTTCATCGAATAACCGCGCCACCTGCGTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCt                                                                     <  3:81417/149‑1 (MQ=37)                                                                        gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTAt                                                                       <  3:317302/147‑1 (MQ=255)                                                                            gAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCACCACACGTCTCCAAAGCTGTGCCGTACCTGTGGGAGACGGCCTGCCCGCGAATGGCCGCAACCCAATAATCTGGtt                                                                 >  5:110391/1‑149 (MQ=37)                                                                                     gtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                    <  1:108121/137‑1 (MQ=255)                                                                                      gggtCAGCTTCATCGAATACACGCGACACCTGCTGCGCGGGCCCGCCCGCTCCCACAGGTGCAGCACAGGGCGCGAACGCTGGGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGcc                                                       <  5:306376/148‑1 (MQ=38)                                                                                           aGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGcc                                                  <  3:147533/149‑1 (MQ=255)                                                                                                 aTCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaac                                            <  5:108235/149‑1 (MQ=255)                                                                                                 aTCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCCCACCCCCACGGGCAGCCCACGGAGGGGACGCGTTGGCGGACCGTGGGGGGCGGGGGGGGCCGCCGATGGCCGCACACCTTTTATGTGGTTGCCGCCCACCCCGCcaacaac                                            >  6:90292/1‑149 (MQ=255)                                                                                                         acggcgACACCGGCTGCGGGGGCAGGCCCGCTCCCAATGGTGCAGCACCGGTCTCGACCGCTGTGCCGTACCTGTGGGAGGGGGGGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                    <  1:345767/146‑1 (MQ=255)                                                                                                              gACACCGGCTCCGCGGGCACGCCCGCTCCCACCGGGCCCGCACAGGTCTCGACCGCGGTGCCGTACCTGGGGGTGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                               <  6:21221/149‑1 (MQ=37)                                                                                                                    tGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCGGTGGGAGCGGGCGGGCCCGCGAATGGCCGCAACCCGTTTATCTGGTTGCCGCCAGCCCCGCCAACCACCAGCGCTGGTAAAgttgtt                         >  8:190291/1‑149 (MQ=37)                                                                                                                      cTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttcttc                        <  7:190291/148‑1 (MQ=255)                                                                                                                         cgcgGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                     <  8:287972/148‑1 (MQ=255)                                                                                                                          gcgGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCTAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                   <  8:62173/149‑1 (MQ=255)                                                                                                                           cgGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAg                  >  1:41887/1‑149 (MQ=255)                                                                                                                            gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGc                 >  4:143639/1‑149 (MQ=255)                                                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc           <  1:317220/149‑1 (MQ=255)                                                                                                                                    ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg         >  7:123264/1‑149 (MQ=255)                                                                                                                                    ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg         >  6:269290/1‑149 (MQ=255)                                                                                                                                    ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTCCTTTAAGCCCCTGCgg         >  7:276765/1‑149 (MQ=255)                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc   >  2:21976/1‑149 (MQ=255)                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc   >  3:2936/1‑149 (MQ=255)                                                                                                                                          ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc   >  2:105440/1‑149 (MQ=21)                                                                                                                                          caaaagttGCAGCAAAGGTCTCGAACGCTGTGCCGTAAATGTGGGAGCGGGAGTGCACGCGAATGGACGAAACAATTTTATCTGGTTGCAGCCAGACCCGCCAACACAAAGCGCTTGTAAAGTTCTTCTTTAAGCCACTGCGGATGcc    >  2:27821/8‑148 (MQ=255)                                                                                                                                           cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCACCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc   >  7:107403/1‑148 (MQ=21)                                                                                                                                            acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc  >  6:152988/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                            TCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC  >  NC_002947/6013691‑6013977

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