Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,096 C→T 10.5% intergenic (+64/‑66) 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,0960CT10.5% 101.5 / 2.5 38intergenic (+64/‑66)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base T (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

ATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357949‑358242
                                                                                                                                                     |                                                                                                                                                  
aTGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATcc                                                                                                                                                   <  2:266330/149‑1 (MQ=255)
   cgcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                 <  3:330755/148‑1 (MQ=255)
    gcgATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGt                                                                                                                                               <  5:284623/149‑1 (MQ=255)
     cgATACCGGCGCCTTTATGCACGACATGTACTAAGAGAGATCGCGAGGCGGGCTGGCTGTCAGCGCTATAGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTa                                                                                                                                              <  6:377172/149‑1 (MQ=255)
        tACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGa                                                                                                                                           <  5:345476/149‑1 (MQ=37)
            ggCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGTCTCCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGc                                                                                                                                         <  1:117227/149‑1 (MQ=25)
              cgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCgggt                                                                                                                                     <  4:4644/149‑2 (MQ=255)
              cgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGc                                                                                                                                     <  4:376321/149‑1 (MQ=255)
                         aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGccc                                                                                                                               <  7:211655/144‑1 (MQ=255)
                         aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                          <  7:308389/149‑1 (MQ=255)
                         aCGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCCCTCCGACAGGCACTGTTCAATCCCTGGAGGATCGGGCTTGCCCgagat                                                                                                                          >  8:211655/1‑148 (MQ=38)
                          ctactaGTACAATTAGACCTGGGGCGGGGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTTTGTGATG‑‑ACCTTGTGTGCCTGTTCGCGGGCAATCCCGCTCCTACAGTAAGTTTGAAATCCCTGTAGTTGCGGGTTTGCCCGCGAag                                                                                                                         <  6:217872/144‑1 (MQ=255)
                                            cTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTTCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                                       >  8:19429/1‑149 (MQ=255)
                                               gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAg                                                                                                                                             >  2:250469/1‑108 (MQ=255)
                                               gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAg                                                                                                                                             <  1:250469/108‑1 (MQ=255)
                                                  cggcggCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTg                                                                                                 <  7:19429/149‑1 (MQ=255)
                                                                                                           gCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGTAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATg                                                           <  5:140597/132‑1 (MQ=255)
                                                                                                           gCCTCTTCGCGGGCAAGCCCCCAACTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCGATGCCACTTTATTGTGCGGTGATg                                                           >  6:140597/1‑132 (MQ=37)
                                                                                                              tctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCCATCCCAGTTGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGAAGCCGGTGCCACACAACCGCCGATGCCTCGTTATGGTGCGGTGAATGGCCACCCGCCTAACGAtcg                                       >  6:163709/1‑149 (MQ=25)
                                                                                                                ctCGCGGGCACGCCCGCTCATACAGGAACTGTGCAATCCCTTTAGGCGCGGGGTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgt                                      <  4:279868/147‑1 (MQ=255)
                                                                                                                   gcgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgg                                      >  3:279868/1‑144 (MQ=255)
                                                                                                                     gggCAAGCCCGCTCCTACAGGATCAGTGCAATCCCTGTAGGAGCGGGCTTGACCGCGAACAGGCCCTTGAGGATGATGAAGGAGAAACACAACCGCCAATGCCTCTTTATGGTGCGGTGATGGGCGCCCCCCATACCGAAAGTCGGAt                                 >  7:59014/1‑148 (MQ=37)
                                                                                                                             ccGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGATGATGCAGGTGCAACACAACCGCCTATGAAAATTTATGGTGCGGTGATGGGCGCCCCGCAAACCGAACGACGGATGACGCCAt                         >  8:244558/1‑148 (MQ=18)
                                                                                                                               gCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGACCTTGAGTCTGATGCAGGTGCAACACAACCGACTATGCATCTTTATGGTGCGGTGATGGGAGCCCCTCCTACCGATCGTCGGATTTCGCGATAAg                      >  4:209851/1‑149 (MQ=38)
                                                                                                                                cgccTCCAGCACCTGTGCAATCCCTGTCGGCGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgt                     <  1:150384/147‑1 (MQ=255)
                                                                                                                                 tCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAgtg                    >  4:30000/1‑149 (MQ=255)
                                                                                                                                   cTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg                  >  8:319172/1‑149 (MQ=37)
                                                                                                                                   cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGCTCGTCGGATGTCGCGATCAGtgtg                  >  8:359053/1‑149 (MQ=255)
                                                                                                                                    tACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                 >  4:173793/1‑149 (MQ=255)
                                                                                                                                      cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc               <  3:173793/149‑1 (MQ=255)
                                                                                                                                         gAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                                 >  1:168873/1‑128 (MQ=255)
                                                                                                                                         gAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                                 <  2:168873/128‑1 (MQ=255)
                                                                                                                                             tgtgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGGACAGTGTGGGCGTTGTCg        >  1:135799/1‑149 (MQ=255)
                                                                                                                                               tgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  6:315130/1‑149 (MQ=255)
                                                                                                                                                 cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      >  6:363944/1‑147 (MQ=255)
                                                                                                                                                 cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt      <  5:363944/147‑1 (MQ=255)
                                                                                                                                                  aaTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                         <  1:31591/127‑1 (MQ=255)
                                                                                                                                                  aaTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTCATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                         >  2:31591/1‑127 (MQ=255)
                                                                                                                                                   accccGGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAGCAGTGTGGGCGTTGTCGGTGaca  <  8:59014/147‑1 (MQ=255)
                                                                                                                                                   aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGaca  >  2:133709/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACA  >  NC_002947/357949‑358242

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