Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,194,935 T→G 26.2% intergenic (+29/+24) PP_4574 → / ← PP_4575 hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,194,9350TG26.2% 24.3 / 12.3 23intergenic (+29/+24)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base G (3/3);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

CTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCG  >  NC_002947/5194798‑5195074
                                                                                                                                         |                                                                                                                                           
cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGAAAg                                                                                                                                 <  2:510135/150‑1 (MQ=255)
         nacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTNTCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCACCCGCGAAAGGGCcagccg                                                                                                                        >  1:85778/2‑148 (MQ=255)
                 cGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGa                                                                                                                <  2:541257/150‑1 (MQ=255)
                          tGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGccc                                                                                                       >  1:79174/1‑150 (MQ=255)
                                gcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                  <  1:138565/117‑1 (MQ=255)
                                gcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                  >  2:138565/1‑117 (MQ=255)
                                         tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctg                                                                                        <  2:79174/150‑1 (MQ=255)
                                                     tGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGc                                                                            <  2:299086/150‑1 (MQ=255)
                                                            gCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGcc                                                                                <  2:363985/139‑1 (MQ=255)
                                                            gCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGcc                                                                                >  1:363985/1‑139 (MQ=255)
                                                                           gTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGGCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTg                                                      >  2:485272/1‑150 (MQ=255)
                                                                                                      cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCa                                <  1:549665/145‑1 (MQ=255)
                                                                                                      cAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCa                                >  2:549665/1‑145 (MQ=255)
                                                                                                                    tGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCATGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgc                        <  1:231025/139‑1 (MQ=255)
                                                                                                                    tGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgc                        >  2:231025/1‑139 (MQ=255)
                                                                                                                    tGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgc             >  1:431990/1‑150 (MQ=255)
                                                                                                                        ccGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagc         <  1:532900/150‑1 (MQ=255)
                                                                                                                        ccGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagc         >  1:4861/1‑150 (MQ=255)
                                                                                                                               tttCGCGGGCGAGCCCGTGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                      >  2:229222/1‑82 (MQ=255)
                                                                                                                               tttCGCGGGCGAGCCCGTGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                      <  1:229222/82‑1 (MQ=255)
                                                                                                                               tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGccagcgccagcg  >  2:176822/1‑150 (MQ=255)
                                                                                                                                    cgGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCAc                           <  1:606687/120‑1 (MQ=255)
                                                                                                                                    cgGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCAc                           >  2:606687/1‑120 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCG  >  NC_002947/5194798‑5195074

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