Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,013,836 G→A 13.1% 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_0029476,013,8360GA13.1% 98.2 / 5.3 38intergenic (+41/+69)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (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 = 9.92e-01

AGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAC  >  NC_002947/6013689‑6013980
                                                                                                                                                   |                                                                                                                                                   
aGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgc                                                                                                                                                    <  8:178817/149‑1 (MQ=255)
 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                                                                                                                            <  6:60238/140‑1 (MQ=255)
                                aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGc                                                                                                                            >  5:60238/1‑140 (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      >  1:239309/1‑149 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc      >  6:289538/1‑149 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc      >  6:148393/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)
                                                                                                                                                   |                                                                                                                                                   
AGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAC‑GCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAA‑‑TGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAC  >  NC_002947/6013689‑6013980

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