Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,013,838 A→G 13.1% 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_0029476,013,8380AG13.1% 93.8 / 5.2 38intergenic (+43/+67)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base G (3/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATC  >  NC_002947/6013690‑6013986
                                                                                                                                                    |                                                                                                                                                       
gTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgca                                                                                                                                                         <  1:38713/149‑1 (MQ=255)
   gTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcagca                                                                                                                                                      >  4:267363/1‑149 (MQ=255)
   gTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcagc                                                                                                                                                       <  7:146224/148‑1 (MQ=255)
    ttCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGcac                                                                                                                                                     <  7:29624/149‑1 (MQ=255)
    ttCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGcac                                                                                                                                                     <  5:189323/149‑1 (MQ=255)
       cAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgca                                                                                                                                                         <  4:227726/142‑1 (MQ=255)
       cAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgca                                                                                                                                                         >  3:227726/1‑142 (MQ=255)
       cAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAg                                                                                                                                                   <  5:360899/148‑1 (MQ=255)
          cgcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCCTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGCTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtct                                                                                                                                               >  7:263005/2‑149 (MQ=255)
           gcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                              <  3:267363/149‑1 (MQ=255)
             ggcCGTGCCGAACTGGTCAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAgcgttcga                                                                                                                                            <  1:408090/149‑8 (MQ=39)
              gcCGTGCCGAACTGGACAGCCAGGGTTACCGCTTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                           <  8:263005/149‑1 (MQ=255)
               ccGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                          >  4:372880/1‑149 (MQ=255)
                cGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑g                                                                                                                                        <  3:372880/149‑1 (MQ=255)
                      gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  >  5:295848/1‑149 (MQ=255)
                      gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  >  7:232310/1‑149 (MQ=255)
                      gAACTGGACAGCCAGGGTTACCGCGGTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTTGCGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCGCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  <  8:339059/149‑1 (MQ=255)
                               aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGACGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCtgt                                                                                                                         <  8:212360/149‑1 (MQ=255)
                               aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  >  5:60238/1‑140 (MQ=255)
                               aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  <  6:60238/140‑1 (MQ=255)
                               aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCtgt                                                                                                                         <  6:361298/149‑1 (MQ=255)
                                gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGA‑CCTGTGCTGCACCtgtg                                                                                                                        <  1:265535/149‑1 (MQ=25)
                                 ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTgg                                                                                                                       <  3:81390/149‑1 (MQ=255)
                                  cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTgg                                                                                                                       <  7:35930/148‑1 (MQ=255)
                                       ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCgg                                                                                                                 <  8:232310/149‑1 (MQ=255)
                                       ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCgg                                                                                                                 <  6:295848/149‑1 (MQ=255)
                                             cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCCTACCTGTGGGAGCGGGCGTGc                                                                                                           <  6:363869/149‑1 (MQ=255)
                                                 gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  >  7:153745/1‑122 (MQ=255)
                                                 gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                                  <  8:153745/122‑1 (MQ=255)
                                                                  gAAGTCGCTGGATGGTGGTGTGGGCAGATTCATCGAATAAACGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGACC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCtt                                                                                    <  2:320928/149‑1 (MQ=37)
                                                                                gtggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                               <  2:192299/110‑1 (MQ=255)
                                                                                gtggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                               >  1:192299/1‑110 (MQ=255)
                                                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                    >  6:240836/1‑149 (MQ=255)
                                                                                                                                    ccGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GATGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg                  >  3:156535/1‑149 (MQ=255)
                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAACGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc             >  1:17463/1‑149 (MQ=255)
                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc            >  6:289538/1‑149 (MQ=255)
                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc            >  6:148393/1‑149 (MQ=255)
                                                                                                                                          ccACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc            >  1:239309/1‑149 (MQ=255)
                                                                                                                                          ccACAGGTACGGCACAGCGTTCGAGA‑CCTGTGCTGCACCTGTGGCAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTAGCAGCCAGCCCCGCCAACAACAACCCCTGGTAAAGATCTTCTTTAAGCCCCTGCGGCTGcc             >  2:107906/1‑148 (MQ=18)
                                                                                                                                           cacaGGTACGGCACAGCGTTCGAGA‑CCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTCGTAAAGTTCTTCTTTAAGCACATGCGGCTGccc            >  8:304633/1‑148 (MQ=18)
                                                                                                                                            acaGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGAGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc            >  1:115948/1‑149 (MQ=39)
                                                                                                                                            acaGGTACGGCACAGCGTTCGAGA‑CCTGTGCTGCACCTGTGGGAGTGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc           >  5:10368/1‑148 (MQ=21)
                                                                                                                                              aGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc        >  3:168887/1‑149 (MQ=255)
                                                                                                                                                    agcaCAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc  >  8:252001/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATC  >  NC_002947/6013690‑6013986

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