Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,013,846 T→G 14.4% intergenic (+51/+59) 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,8460TG14.4% 87.8 / 6.8 34intergenic (+51/+59)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base G (3/2);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCC  >  NC_002947/6013700‑6013991
                                                                                                                                                  |                                                                                                                                                    
cgcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCCTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGCTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtct                                                                                                                                                    >  7:263005/2‑149 (MQ=255)
 gcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                                   <  3:267363/149‑1 (MQ=255)
    gcCGTGCCGAACTGGACAGCCAGGGTTACCGCTTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                                <  8:263005/149‑1 (MQ=255)
     ccGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                               >  4:372880/1‑149 (MQ=255)
      cGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACg                                                                                                                                             <  3:372880/149‑1 (MQ=255)
            gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       >  7:232310/1‑149 (MQ=255)
            gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       >  5:295848/1‑149 (MQ=255)
            gAACTGGACAGCCAGGGTTACCGCGGTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTTGCGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCGCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       <  8:339059/149‑1 (MQ=255)
                     aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGACGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                              <  8:212360/149‑1 (MQ=255)
                     aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       <  6:60238/140‑1 (MQ=255)
                     aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       >  5:60238/1‑140 (MQ=255)
                     aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                              <  6:361298/149‑1 (MQ=255)
                      gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                             <  1:265535/149‑1 (MQ=25)
                       ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                            <  3:81390/149‑1 (MQ=255)
                        cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                            <  7:35930/148‑1 (MQ=255)
                             ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                      <  6:295848/149‑1 (MQ=255)
                             ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                      <  8:232310/149‑1 (MQ=255)
                                   cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCCTACCTGTGGGAGCGGGCGTGc                                                                                                                <  6:363869/149‑1 (MQ=255)
                                       gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       <  8:153745/122‑1 (MQ=255)
                                       gAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                       >  7:153745/1‑122 (MQ=255)
                                                        gAAGTCGCTGGATGGTGGTGTGGGCAGATTCATCGAATAAACGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCtt                                                                                         <  2:320928/149‑1 (MQ=37)
                                                                      gtggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                    <  2:192299/110‑1 (MQ=255)
                                                                      gtggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                    >  1:192299/1‑110 (MQ=255)
                                                                                                                        gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                         >  6:240836/1‑149 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGATGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg                       >  3:156535/1‑149 (MQ=255)
                                                                                                                                ccACAGGTGCAGCACAGGTCTCGAACGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                  >  1:17463/1‑149 (MQ=255)
                                                                                                                                ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                 >  6:289538/1‑149 (MQ=255)
                                                                                                                                ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                 >  6:148393/1‑149 (MQ=255)
                                                                                                                                ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                 >  1:239309/1‑149 (MQ=255)
                                                                                                                                ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGCAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTAGCAGCCAGCCCCGCCAACAACAACCCCTGGTAAAGATCTTCTTTAAGCCCCTGCGGCTGcc                  >  2:107906/1‑148 (MQ=18)
                                                                                                                                 cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTCGTAAAGTTCTTCTTTAAGCACATGCGGCTGccc                 >  8:304633/1‑148 (MQ=18)
                                                                                                                                  acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGAGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                 >  1:115948/1‑149 (MQ=39)
                                                                                                                                  acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGTGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                >  5:10368/1‑148 (MQ=21)
                                                                                                                                    aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc             >  3:168887/1‑149 (MQ=255)
                                                                                                                                          agcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc       >  8:252001/1‑149 (MQ=255)
                                                                                                                                                 gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGccc  >  3:269992/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCC  >  NC_002947/6013700‑6013991

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: