Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,013,854	G→C	16.5%	intergenic (+59/+51)	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,854	0	G	C	16.5%	111.2 / 14.8	48	intergenic (+59/+51)	xpt/PP_5266	xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base C (4/4);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

CGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTT  >  NC_002947/6013706‑6014001                                                                                                                                                     |                                                                                                                                                    cGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACg                                                                                                                                                     <  6:108055/149‑1 (MQ=255)             gACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                          <  5:140763/148‑1 (MQ=255)                aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                      <  1:101432/149‑1 (MQ=255)                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                                     <  6:358053/149‑1 (MQ=25)                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                                     <  6:1618/149‑1 (MQ=25)                           ccGCGTTGAATCGCTGGCCCGGGGGAAGTCGCTGGAAGGGGGTGGGGTCAGCTTCATCGAATAACCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                            <  7:245665/148‑1 (MQ=255)                                 tGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGcccg                                                                                                                     >  1:64623/1‑147 (MQ=255)                                     tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                      >  1:105714/1‑128 (MQ=255)                                     tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                      <  2:105714/128‑1 (MQ=255)                                         tGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCGCACCCACAGGTGCAGCACAGGTCTAGAACGCGGTCCCGTACCTGTGGGCGCGGGCGTGACCGCGAATGGc                                                                                                             >  8:92608/1‑149 (MQ=255)                                              cGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                        >  5:111691/1‑149 (MQ=255)                                                ggTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAAcc                                                                                                      <  8:241501/149‑1 (MQ=255)                                                   gAAGGCGCTGGAAGGGGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtt                                                                                                   <  7:92608/149‑1 (MQ=255)                                                        cGCTGGAAGGTGGTGTGGTCAGCTTCATCTAATAACCGCTACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGTTCTCGAACGCTGTGCGTTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTAt                                                                                               <  8:139563/148‑1 (MQ=37)                                                         gCTGGATGTGGGTGTGTTCTGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                              <  2:162027/148‑1 (MQ=37)                                                          cTGGACGGGTGTGTGTGCAGCTTCATCTATTAAACTCGACACCTGCTGCGCGGGCACGCCCGCCCCCACAGGTGCAGCACAGGTCTCGACCGCGGTGCCGTACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                                            <  6:34855/149‑1 (MQ=25)                                                           tGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCCCGCCCGCTCCCACACGTGCACCACAGGGCTCGAACCCGGTGCCGTACCAGTGGGAGCGGGCGTGCCCGCGAATGGCCGCCACCCGTTTTTCTgg                                                                                           >  2:77774/1‑149 (MQ=37)                                                           tGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTgg                                                                                           <  2:64623/149‑1 (MQ=255)                                                             gacggtggtGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGt                                                                                          <  7:228452/145‑1 (MQ=255)                                                             gAAGGGGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCTCAGGTCTAGATCGCTTTGAGGTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGtt                                                                                         <  5:266842/149‑1 (MQ=38)                                                                  tggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccgc                                                                                    <  6:111691/149‑1 (MQ=255)                                                                     tgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCGCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAg                                                                                 <  7:39664/149‑1 (MQ=255)                                                                                                     tGCTTCGCGGGCACGCCCGCTACCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGGGCCCGCGAATGGCCCCAACCCGTTTATATGGTTGCCGCCAGCCCCGCCAACCACAAGCGCTGGTAAAGttgt                                                  >  4:340566/1‑146 (MQ=37)                                                                                                             gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGc                                         >  4:36024/1‑149 (MQ=255)                                                                                                                 aCGCCCGCTACCACAGGTGCCGGACAGGTCTCGAACGGTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCt                                     <  5:134622/149‑1 (MQ=39)                                                                                                                   gCCCGCTCCCCCCGGTGCAGCACAGGGCTCGAACGCTGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                   <  1:353563/149‑1 (MQ=39)                                                                                                                   gCCCGCTCCCACAGTGGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                   <  1:347901/149‑1 (MQ=255)                                                                                                                   gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                   >  4:327293/1‑149 (MQ=255)                                                                                                                    cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                  <  3:36024/149‑1 (MQ=255)                                                                                                                          cccACAGGTGCAGCACAGGTCTCGAACGCGGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                            <  3:224334/149‑1 (MQ=255)                                                                                                                          cccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaac                                                                    >  6:332565/1‑109 (MQ=17)                                                                                                                          cccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaac                                                                    <  5:332565/109‑1 (MQ=17)                                                                                                                           ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                           >  6:38264/1‑149 (MQ=21)                                                                                                                            cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                  >  6:305974/1‑141 (MQ=21)                                                                                                                            cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                  <  5:305974/141‑1 (MQ=21)                                                                                                                             acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                          >  5:96108/1‑148 (MQ=255)                                                                                                                              caGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                        <  2:112824/149‑1 (MQ=255)                                                                                                                                ggTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACTAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                      <  1:77774/149‑1 (MQ=255)                                                                                                                                 gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTg                     >  7:173325/1‑149 (MQ=255)                                                                                                                                      gcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                             <  7:242196/136‑1 (MQ=255)                                                                                                                                      gcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                             >  8:242196/1‑136 (MQ=255)                                                                                                                                        acaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                   <  5:296569/128‑1 (MQ=255)                                                                                                                                        acaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                   >  6:296569/1‑128 (MQ=255)                                                                                                                                          aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGccc            >  7:253010/1‑149 (MQ=255)                                                                                                                                              ctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCACAACTGCATCCGCCCCAgg        >  6:318993/1‑149 (MQ=255)                                                                                                                                               tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTCGTAAAGTTCTACTTTAAGCCCCTGCGGCTGCCCCAACTGCAACCGCCCCATGt       >  6:82405/1‑149 (MQ=39)                                                                                                                                                 gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg     >  6:353057/1‑149 (MQ=255)                                                                                                                                                  aaCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGc    >  8:42353/1‑149 (MQ=255)                                                                                                                                                    cGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCtt  <  6:96108/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTT  >  NC_002947/6013706‑6014001

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