Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,013,861 C→T 16.4% intergenic (+66/+44) 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,8610CT16.4% 51.6 / 2.7 25intergenic (+66/+44)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (11/9);  new base T (2/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAC  >  NC_002947/6013722‑6014007
                                                                                                                                           |                                                                                                                                                        
gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGCCCTGTGCTGCACCtgtg                                                                                                                                                 >  8:328339/1‑149 (MQ=25)
      gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCg                                                                                                                                           <  7:328339/149‑1 (MQ=21)
          ccGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGCTGTGGGCAGCTTCATCGAATAAACGCGACACCTGCGTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                                        <  7:174463/148‑1 (MQ=255)
                  aaTCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTggg                                                                                                                                               <  8:33166/133‑1 (MQ=21)
                    tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                      >  5:72011/1‑124 (MQ=255)
                    tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                      <  6:72011/124‑1 (MQ=255)
                                     gTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                                        <  4:134592/121‑1 (MQ=255)
                                     gTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                                        >  3:134592/1‑121 (MQ=255)
                                                                         aaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTg                                                                  <  1:316397/149‑1 (MQ=255)
                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                                 <  1:99400/51‑1 (MQ=255)
                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                                 >  2:99400/1‑51 (MQ=255)
                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                         >  2:100837/1‑149 (MQ=255)
                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                         >  2:282444/1‑149 (MQ=255)
                                                                                                  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                         >  2:265939/1‑149 (MQ=255)
                                                                                                          ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                 >  5:21880/1‑149 (MQ=255)
                                                                                                          ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                 >  1:172019/1‑149 (MQ=21)
                                                                                                            acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                               >  1:123819/1‑149 (MQ=255)
                                                                                                               ggTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                            <  6:21880/149‑1 (MQ=255)
                                                                                                                gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTg                           >  6:115000/1‑149 (MQ=255)
                                                                                                                           gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGccc                  >  1:132735/1‑147 (MQ=255)
                                                                                                                                gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTg            <  1:265939/148‑1 (MQ=255)
                                                                                                                                gAACGCTGTGCCGTAACTGTGGGAGCGGGCGGGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACACTTCAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg              <  1:100837/148‑1 (MQ=255)
                                                                                                                                 aaCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGc          <  2:123819/149‑1 (MQ=255)
                                                                                                                                  aCGCTGTGCGGTACCTGTGGTTGGGGGCGTGCCCGCTAATGGCCGAAACCCTTTTTTCATTGGTTGCGGCATCTCGATCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTACCCCAACTGCATCCGCCCCAgg              <  1:282444/148‑1 (MQ=255)
                                                                                                                                      tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa     >  7:327714/1‑149 (MQ=255)
                                                                                                                                          ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                        >  8:126178/1‑110 (MQ=255)
                                                                                                                                          ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                            >  4:21674/1‑122 (MQ=255)
                                                                                                                                          ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                            <  3:21674/122‑1 (MQ=255)
                                                                                                                                          ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAc  <  8:327714/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                        
GCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATC‑‑TGGTTGCCGC‑‑CAGCCCCGCCAACA‑‑ACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAC  >  NC_002947/6013722‑6014007

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