Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,085 A→G 11.0% intergenic (+53/‑77) tdcG‑I → / → PP_5441 L‑serine dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947358,0850AG11.0% 83.9 / 2.2 36intergenic (+53/‑77)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/13);  new base G (2/2);  total (21/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 = 1.69e-01

TAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTG  >  NC_002947/357938‑358233
                                                                                                                                                   |                                                                                                                                                      
tAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGaa                                                                                                                                                       <  8:289749/149‑1 (MQ=255)
  aTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACt                                                                                                                                                     <  6:356806/149‑1 (MQ=255)
  aTCCGAAGCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGAGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                      <  6:283216/148‑1 (MQ=255)
          cATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATc                                                                                                                                             >  5:355120/1‑149 (MQ=255)
                     ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAg                                                                                                                                   <  8:84108/148‑1 (MQ=37)
                      cGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTATTCGCGGGCAAGCCCGCTCCTAAACGAACAGTGGAATCCCTGTAGGAGCg                                                                                                                                 >  2:279441/1‑149 (MQ=255)
                       ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCgg                                                                                                                                <  4:331585/149‑1 (MQ=37)
                                         aaGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCCCGGGCAAGCCCCCACCTACAGGAACAGTGCCAAACCTGTAGGAGCCGGCTTGCCCGCGAACAGG‑c                                                                                                             >  2:42723/1‑149 (MQ=255)
                                           gTACAAAGAGACCTCGCGCGGCGCCCTGGCCGTCACCGCTATCGAGAGCATATTTTCTGACCTGACCTTCTGCGCCTCTTCGCGGGCAAGCCCCATACTTCACGAACAGTACATTACCAGTACGAGCCGGCTTTCCCCCCAAGCGG‑ccc                                                                                                           >  6:358143/1‑149 (MQ=255)
                                            tACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑ccca                                                                                                          >  1:144488/1‑148 (MQ=255)
                                               aaaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGACAGGGACTGTGGAATGCCGGTAGGCGCGGGCGTGGCCGCGGAGAGG‑CCCGTgg                                                                                                       >  1:196953/1‑148 (MQ=37)
                                               aaaGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCGACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCtt                                                                                                         >  6:168882/1‑147 (MQ=255)
                                                       cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATg                                                                                               >  7:113998/1‑149 (MQ=255)
                                                        tCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTTCAGGAACGGTGCAATACCGGTAGGAGCGGGCTTGCCCGCGCAGCGG‑CCCGTGAGGCTGATg                                                                                               >  7:19069/1‑148 (MQ=38)
                                                             cggcggCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTg                                                                                         <  6:355120/149‑1 (MQ=255)
                                                                 ggcggggcggcCCGCGCTATCCTGTGGTTATTTTCTGCCCTCCCTTGCTGGGCCTCTTCGCGGCAACCCCCGCCCCTCCCGAAACTGTCCAACCCGTGTCGGCGGGGGCTGCCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAac                                                                                     <  2:196953/144‑1 (MQ=255)
                                                                         gggccgcgcTATCGAGTGCTGTTTTTCTTAGCGGACCTTCTGGGCCTCTTCGCGGCCAAGCCGGCTCCTACCGCAACTGTGCAATCCCTGTAGGAGGGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGc                                                                             <  8:19069/144‑1 (MQ=37)
                                                                             aGCGCTATCGTGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTa                                                                          <  1:42723/149‑1 (MQ=255)
                                                                             aGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTAt                                                                         <  2:116255/149‑1 (MQ=255)
                                                                                           cTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTg                                                          >  5:170378/1‑149 (MQ=255)
                                                                                            tGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGc                                                         >  5:338101/1‑149 (MQ=255)
                                                                                            tGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCAACAGGAACAGTGCAATCCCTGGAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGAAACACAACCGCCTATGCCGC‑TTTAAGGTGc                                                         >  5:326239/1‑149 (MQ=38)
                                                                                              tttttttGAGCTCACCTGCTGGGCTGTTTCGCGGCAATGCCCGGTCATACTGCAATGTTGCAAACCATGTTGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCg                                                        <  1:279441/147‑1 (MQ=255)
                                                                                               ttttCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCgg                                                       <  6:338101/148‑1 (MQ=255)
                                                                                                   cTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACTGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATg                                                  >  1:254741/1‑149 (MQ=255)
                                                                                                   cTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGAt                                                   <  6:170378/148‑1 (MQ=255)
                                                                                                                cTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCg                                         <  7:36522/145‑1 (MQ=255)
                                                                                                                cTGGGCCTCTTCGCGGGCAAAACCGCTCCCACAGGAACAGTGGAATCCCGGTTAGAACCCGCATGGCCCCCAAGACG‑CCCGTGGCGCGGATGCACCTGCAACCAAACCGCCAAAGGCTCGTTTAT‑GTTAGGGGATGGGCGgcccc                                         >  8:36522/1‑144 (MQ=255)
                                                                                                                       tctTCGCGGGCAAGCCCGCAACGACAGGAAAAGTGCAATACCAGGAGGAGCGGGGTTGGCCGCGAAGAGG‑CCCTTGAGGCTGATGCAAGGGAAACACAACCGCCAATTCCTA‑TTTATGGTGGGGTGATGGGCGCCCCGCCAACCGCtcg                              >  2:278785/1‑149 (MQ=255)
                                                                                                                         ttCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCGGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCGTGAGGCTGCTGCAGGTGCAACACACCCGCCTATGCCTC‑GTTATGGTGCGGGGCTGGGCGCCCCGCCTACCGAacgtc                            >  3:296624/1‑149 (MQ=37)
                                                                                                                           cgcgGGCAAGCCCGCTCCTACAGGACCTGTGCAAGCCCTGTCGTAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgt                             <  1:232059/146‑1 (MQ=255)
                                                                                                                            gcgGGCAAGCCCGCTACTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGCTCGTCGGa                         >  4:140029/1‑149 (MQ=255)
                                                                                                                                gCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTc                     >  7:273856/1‑149 (MQ=255)
                                                                                                                                gCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTc                     >  5:28433/1‑149 (MQ=255)
                                                                                                                                        gCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAg             >  7:345105/1‑149 (MQ=37)
                                                                                                                                          tCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg           >  1:95734/1‑149 (MQ=255)
                                                                                                                                            cTACTGGAACTGGGCAATCCCGGTCGGCGGGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg         <  6:326239/149‑1 (MQ=255)
                                                                                                                                                   aaCTGTGCAATCCCTGTTGGTGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTg  <  4:296624/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGG‑CCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTC‑TTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTG  >  NC_002947/357938‑358233

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