Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,091 G→A 12.5% intergenic (+59/‑71) 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,0910GA12.5% 75.2 / 3.0 33intergenic (+59/‑71)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/17);  new base A (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGT  >  NC_002947/357943‑358238
                                                                                                                                                      |                                                                                                                                                   
cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                                     <  1:98769/149‑1 (MQ=255)
 gCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                                     >  4:330755/1‑148 (MQ=255)
    ccATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAAt                                                                                                                                                 <  4:388281/149‑1 (MQ=255)
      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                             <  2:168873/128‑1 (MQ=255)
                                                                                                                                               gAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                             >  1:168873/1‑128 (MQ=255)
                                                                                                                                                   tgtgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGGACAGTGTGGGCGTTGTCg    >  1:135799/1‑149 (MQ=255)
                                                                                                                                                     tgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGt  >  6:315130/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTG‑‑ACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGT  >  NC_002947/357943‑358238

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: