Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,013,847 C→T 14.9% intergenic (+52/+58) 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,8470CT14.9% 79.6 / 5.8 37intergenic (+52/+58)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (8/21);  new base T (2/3);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC  >  NC_002947/6013704‑6013977
                                                                                                                                                |                                                                                                                                  
gcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                               <  4:138203/149‑1 (MQ=255)
       cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                        <  2:204644/149‑1 (MQ=255)
        gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                       <  3:149449/149‑1 (MQ=255)
              gACAGCCAGGGTTACCGCGTTGAATCGCCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgt                                                                                                                         >  1:175136/1‑142 (MQ=255)
              gACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgt                                                                                                                         <  2:175136/141‑1 (MQ=255)
                   ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                            <  8:3398/149‑1 (MQ=255)
                     aGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGa                                                                                                          <  5:259517/149‑1 (MQ=25)
                         ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCgg                                                                                                      <  1:296415/149‑1 (MQ=21)
                            ccGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                    <  8:138690/148‑1 (MQ=255)
                             cgcgTTTAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGt                                                                                                  <  5:108302/149‑1 (MQ=255)
                             cgcgTTGAATCGCTGGCCCGGGTGAAGTCGCGGGAAGGTGGTGTGGTCAGCTTCATCGAATAACCGCGACACCTGCTGCGCGGGCACGCCCGTTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGt                                                                                                  <  1:239058/149‑1 (MQ=255)
                               cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCACCCCCAGCTTCAGCACACGGCTCCAAACCAGTGCCGTACCTGAGCGACCCGGCGTGc                                                                                                >  3:272273/1‑149 (MQ=255)
                                      tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCTACCACCGGTGCAGCACAGGTGTGGAACGCGGTGGCGTAACGGTGGGAGCGGGCGTGGCcgcgca                                                                                         >  7:42416/1‑147 (MQ=255)
                                             gCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                    <  2:183399/148‑1 (MQ=255)
                                                    tGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCt                                                                            <  3:98128/149‑1 (MQ=255)
                                                        gTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGACCGCCCGCTCCCACAGGGTCAGCCCAGGGATCGCAAGCGGTGGCGGACCTGTGTGTGCGGGCGTGGCCGCGAAAGGCCGCAAACCgtttt                                                                        >  8:136333/1‑148 (MQ=255)
                                                           gCTCGAATGTGGTGTGGTCAGCTTCCTCGAAAAAACCCGACACCTGCTACGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCt                                                                     <  7:136333/149‑1 (MQ=37)
                                                            cgggACGGGGGTGTGGCCAGCTTCATCCAAAAAACGCGACACCTGCTTCGCGGGCCCGCCCGCTCCCCCCGGGGCAGCACAGGCCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                    <  7:249297/147‑1 (MQ=25)
                                                                     ggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGc                                                                                                    >  8:278437/1‑108 (MQ=14)
                                                                     ggtgtggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGc                                                                                                    <  7:278437/108‑1 (MQ=14)
                                                                         tgttCAGCTTCATCTAATAAACGCTACACCTGCTTCGCGGGCACGCCCGCTCCCACGGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGcc                                                       <  8:42416/146‑1 (MQ=38)
                                                                          ggtCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGGCGCCCCGGGCTGGCACGCGGTGGCGGTCCTGTGGGGGCGGGCGGGCCCGCGAATGGGCGCCACCCGGTTATATGGTTGGCGGCAGccc                                                      >  2:83753/1‑149 (MQ=19)
                                                                                cTTCATCGAATAACCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGc                                                   <  1:97220/146‑1 (MQ=255)
                                                                                cTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTACCACAGGTGCAGCAAAGGTATCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCATTTATCTGGTTGCCGCCAGCCCCGCcaa                                                >  7:1146/1‑149 (MQ=38)
                                                                                cTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCCCAACCACACGTGCCGCACAGGTCTAGCACCCGGTGGCGGACCGGTGGGAGCGGGGGGGCCCGCGCAAGGCCGCAACCCGTTTATCTGGTTTGCGCccccccccccca                                                >  8:171760/1‑147 (MQ=255)
                                                                                    aTCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaa                                             >  3:194494/1‑148 (MQ=255)
                                                                                    aTCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCAGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaa                                             >  5:143646/1‑148 (MQ=255)
                                                                                                    acCTGCTTCGCGGGCACGGGCGCTCCCTCTGGTGCAGCACAGGTCTCGATCGCTGTGCCGTACCTGTGTGAGCGGGCGTGCCCGCGAATGGCCTCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGtt                            <  4:272273/149‑1 (MQ=38)
                                                                                                          ttCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTcttctt                      <  4:194494/149‑1 (MQ=255)
                                                                                                                     gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc           <  7:40251/149‑1 (MQ=255)
                                                                                                                     gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc           <  7:75076/149‑1 (MQ=255)
                                                                                                                     gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc           >  8:75076/1‑149 (MQ=255)
                                                                                                                           gcccACCGGTGCAGCCCAGGTCTCGACCGCTGTGCCGTACCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc     <  1:128476/148‑1 (MQ=255)
                                                                                                                            cccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc    >  5:94799/1‑149 (MQ=21)
                                                                                                                             ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc   <  3:141745/149‑1 (MQ=255)
                                                                                                                             ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc    <  7:171760/148‑1 (MQ=255)
                                                                                                                               acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc  >  4:171637/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                  
GCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC  >  NC_002947/6013704‑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: