Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 411,781 C→A 8.7% L172I (CTA→ATA)  kdkA_1 → 3‑deoxy‑D‑manno‑octulosonic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440411,7810CA 8.7% 140.8 / 4.2 46L172I (CTA→ATA) kdkA_13‑deoxy‑D‑manno‑octulosonic acid kinase
Reads supporting (aligned to +/- strand):  ref base C (15/27);  new base A (2/2);  total (17/29)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAGCGCTTGCGCTTCT  >  KT2440/411650‑411913
                                                                                                                                   |                                                                                                                                    
gcTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCACGGCACGCTAt                                                                                                                                   >  5:49022/1‑135 (MQ=255)
gcTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTAt                                                                                                                                   >  3:119755/1‑135 (MQ=255)
            gcgcCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGATATTCTTCATTGAca                                                                                                                       >  6:328433/1‑134 (MQ=255)
            gcgcCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGATATTCTTCATTGAca                                                                                                                       >  6:334410/1‑134 (MQ=255)
                    gTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGTGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCa                                                                                                                 <  7:370204/133‑1 (MQ=255)
                      cGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCCCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCgg                                                                                                             >  3:135765/1‑135 (MQ=255)
                          cATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGAc                                                                                                         <  3:209699/135‑1 (MQ=255)
                           aTGTCAGTCGTCAGCTGGCGCGCCATACACGGGTCATGCACGAACACCGCTTTGCCCATAACGACATGAAGTGGCTCAACCTGCTGGTCGATGACCAACGCACGCTATTCTTCATTGACTGCCCCACCGGCGACt                                                                                                        >  5:179914/1‑135 (MQ=255)
                               cAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTg                                                                                                    <  3:67318/135‑1 (MQ=255)
                               cAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTg                                                                                                    <  7:87733/135‑1 (MQ=255)
                               cAGTCGGCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCt                                                                                                     <  5:45315/134‑1 (MQ=255)
                                           ggCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCAt                                                                                        >  5:80257/1‑135 (MQ=255)
                                            gcgcgcCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCAtg                                                                                       <  2:134166/135‑1 (MQ=255)
                                                      cccGGGTCATGCTCGAACACCGCTTTGCCCATCACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATc                                                                             <  8:371988/135‑1 (MQ=255)
                                                       ccGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATc                                                                             <  4:87314/134‑1 (MQ=255)
                                                        cGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCgg                                                                           <  6:280425/135‑1 (MQ=255)
                                                                 cACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCCTCGGATGATCAAg                                                                  >  4:469/1‑135 (MQ=255)
                                                                     aaCACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGAcc                                                              <  2:30078/135‑1 (MQ=255)
                                                                     aaCACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGAcc                                                              <  7:155105/135‑1 (MQ=255)
                                                                            cTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATg                                                       >  2:295560/1‑135 (MQ=255)
                                                                             tttGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCAAGc                                                      >  6:35692/1‑135 (MQ=255)
                                                                                    aTAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACa                                               <  6:80257/135‑1 (MQ=255)
                                                                                    aTAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTGCATTGACTGCCCCACCTGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACa                                               <  6:30911/135‑1 (MQ=255)
                                                                                       aCGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTATGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGCCCTGGCATGCCTTGACAAAg                                            >  5:328158/1‑135 (MQ=255)
                                                                                        cGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGt                                           <  4:229795/135‑1 (MQ=255)
                                                                                                aGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCaaa                                    <  7:53069/134‑1 (MQ=255)
                                                                                                aGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCaaa                                    <  1:361233/134‑1 (MQ=255)
                                                                                                    gcgcAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATacca                               <  2:64889/135‑1 (MQ=255)
                                                                                                     cgcAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATacca                               <  4:197245/134‑1 (MQ=255)
                                                                                                     cgcAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATacca                               <  4:348665/134‑1 (MQ=255)
                                                                                                       cAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGCTGATCAAGGCCCTGGCATGCCTGGACAAACTGGCCAAATACCACCt                            >  8:172404/1‑135 (MQ=255)
                                                                                                       cAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCCCCt                            >  6:174520/1‑135 (MQ=255)
                                                                                                         aCCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGt                          <  5:179293/135‑1 (MQ=255)
                                                                                                                 gTCGATGCCCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTATGGCGCGGGTTGATGTGGCGCCATCGGCTGATGAAGGACCGGGCCTGCCCGGGCCAAGTGGGCAAATACCACCGGTCCGCTccc                  >  1:391406/1‑135 (MQ=255)
                                                                                                                      tGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcg             <  3:85400/135‑1 (MQ=255)
                                                                                                                              gCACGCTATTCTTCCTTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  5:166217/135‑1 (MQ=255)
                                                                                                                              gCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  6:113320/135‑1 (MQ=255)
                                                                                                                              gCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  5:35692/135‑1 (MQ=255)
                                                                                                                              gCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  2:104680/135‑1 (MQ=255)
                                                                                                                              gCACGATATTCTTCATTGACCGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  5:334410/135‑1 (MQ=255)
                                                                                                                              gCACGATATTCTTCATTGACCGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAgcgcttgcgct     <  5:328433/135‑1 (MQ=255)
                                                                                                                               cACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGt                                           <  5:78173/96‑1 (MQ=255)
                                                                                                                               cACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGt                                           >  6:78173/1‑96 (MQ=255)
                                                                                                                                aCGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCCCCTGTCAGCTACCCAgcgcttgcgctt    >  3:174064/1‑134 (MQ=255)
                                                                                                                                aCGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAGCGCTTGCGCTTc   <  6:62411/135‑1 (MQ=255)
                                                                                                                                 cGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAGCGCTTGCGCTTCt  >  8:27028/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                    
GCTTGGCCGACCGCGCCTGGGTCGACCATGTCAGTCGTCAGCTGGCGCGCCATACCCGGGTCATGCACGAACACCGCTTTGCCCATAACGACCTGAAGTGGCGCAACCTGCTGGTCGATGACCAAGGCACGCTATTCTTCATTGACTGCCCCACCGGCGACTTCTGGCGCGGTTTCATGTGGCGCCATCGGATGATCAAGGACCTGGCATGCCTGGACAAAGTGGCCAAATACCACCTGTCAGCTACCCAGCGCTTGCGCTTCT  >  KT2440/411650‑411913

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: