Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,306,660 G→A intergenic (‑63/+210) yjcS ← / ← alsK linear primary‑alkylsulfatase/D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,6600GA97.1% 118.6 / NA 34intergenic (‑63/+210)yjcS/alsKlinear primary‑alkylsulfatase/D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/19);  minor base T (0/1);  total (14/20)
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

CAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACG  >  NC_000913/4306420‑4306905
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTAt                                                                                                                                                                                                                                             <  1:317631/251‑1 (MQ=255)
cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTAt                                                                                                                                                                                                                                             <  1:504795/251‑1 (MQ=255)
               cGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAAtt                                                                                                                                                                                                                              >  2:173855/1‑251 (MQ=255)
                   aaaCGGGAATGATTTGGAGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCCCTTTCGACTTAGCGGAATTTACCATCATGCCGCTGGCGGCCGGTAAAGCCATAACAATCCTGCTCAAAAGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATAATTTATGTAGTAAATATGTAATACATGCAAAATACCTCATTTAAAAATAGCGGTATTCTAATTGGGAAAATGCgg                                                                                                                                                                                                                          >  2:47973/1‑250 (MQ=255)
                           aTGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTATTTCGCGTTAGCGGTATTTACCATCAGGCCGCTGGCGGCAGTTAACGCAATAACAATCCTTCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTtaata                                                                                                                                                                                                                  <  2:537493/251‑1 (MQ=255)
                                       aaTTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGtgtg                                                                                                                                                                                                      <  1:443013/251‑1 (MQ=255)
                                       aaTTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGACATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGtgtg                                                                                                                                                                                                      <  2:33325/251‑1 (MQ=255)
                                                tAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTTGCGTTAGCGGTATTTTCCATCATGCCGCTGGCGGCCGTTAACGCAATAACCATCCTGCTCAAACGGAATAACAGCGAGGTATTCATATAAATCTACATTTATCATTTATGTCGTAAATATGTAATCCCTGCCAAATACGTCATTTAAAAATGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAAcg                                                                                                                                                                                             >  2:507946/1‑251 (MQ=255)
                                                                       gcggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTtgtg                                                                                                                                                                      <  2:243544/251‑1 (MQ=255)
                                                                                          tCGCGTTATCAGTATTTACCATCATGCCGCGGGCGGCAGTTAACGCAATAACAATCCTTCTCAAACGGAATAAACGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACTGCAGTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAACGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTAc                                                                                                                                                    <  2:18382/250‑1 (MQ=255)
                                                                                                 aGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATGTGTTCTGAAGAGGGAAAATTACAGTaaaa                                                                                                                                             <  1:461799/250‑1 (MQ=255)
                                                                                                     gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATc                                                                                                                                        <  1:284918/251‑1 (MQ=255)
                                                                                                       aTTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGTCGTTATCATAATTGGTAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCGATTTTATTAGGCATTTGGGCTGAAGAGGGAAAATTACAGTAAAGCTATCAt                                                                                                                                      <  2:325601/251‑1 (MQ=255)
                                                                                                           aCCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGc                                                                                                                                   >  1:432562/1‑250 (MQ=255)
                                                                                                                              cAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCGCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGtt                                                                                                               <  2:432562/251‑1 (MQ=255)
                                                                                                                                         aTAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTtta                                                                                                    <  1:173855/251‑1 (MQ=255)
                                                                                                                                             cAATCCTTCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTATGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCAGTGCAGAAAATAATGTGAGTCGGTTTATGTGTTtattat                                                                                                <  2:58576/251‑1 (MQ=255)
                                                                                                                                                    gCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATtatatata                                                                                           <  1:518592/249‑1 (MQ=255)
                                                                                                                                                       cAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTa                                                                                      <  1:410171/251‑1 (MQ=255)
                                                                                                                                                           cGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTg                                                                                  >  1:223550/1‑251 (MQ=255)
                                                                                                                                                                 aaCCGAGATGTATTCATATAAAACTCCATGTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGGCATGTAAATCTGGCGTTATCATAATTGCGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACACGCCAATTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATTATTGCAGATAAGAATATGAGTCGTATAATGTCTTTATTATATTTATGGTACCcgtcgta                                                                             <  2:205984/250‑2 (MQ=255)
                                                                                                                                                                                     aaTCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATGGCAGAAAATAATATGAGTCGGTTTATGTGTTTAGTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCt                                                        <  2:274219/251‑1 (MQ=255)
                                                                                                                                                                                       tctcCATTTATCATTTATGTCGTAAATATGAAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTAc                                                      >  1:428320/1‑251 (MQ=255)
                                                                                                                                                                                         tcCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACtg                                                    >  1:271320/1‑251 (MQ=255)
                                                                                                                                                                                          ccATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACtgt                                                   >  2:5772/1‑251 (MQ=255)
                                                                                                                                                                                                       atGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTa                                                                                      >  1:277297/1‑203 (MQ=255)
                                                                                                                                                                                                       atGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTa                                                                                      <  2:277297/203‑1 (MQ=255)
                                                                                                                                                                                                           cGTAAATATGTAATCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCAt                                  <  1:171916/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTCCATGCAAAATACGTCATTTAAAACTGGCGTTATCATAATTGCGAAATTGCGACTTTAATAAGAGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGc                       >  2:506732/1‑250 (MQ=255)
                                                                                                                                                                                                                               aaaaTACGTCATTTAAAACTGGCGTTATCATAATTGGGAACTTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACtg              >  1:303906/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gTCATTTAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCaa        >  2:125269/1‑250 (MQ=255)
                                                                                                                                                                                                                                           ttAAAACTAGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTCTTAGGCATTTGTGATGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCAAAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACg  >  1:334310/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTg                                                                                  <  2:129924/170‑1 (MQ=255)
                                                                                                                                                                                                                                            tAAAACTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTg                                                                                  >  1:129924/1‑170 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
CAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACG  >  NC_000913/4306420‑4306905

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