Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,205,006 G→T 8.9% P257H (CCC→CAC)  eptA_2 ← Phosphoethanolamine transferase EptA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,205,0060GT 8.9% 126.6 / 4.2 45P257H (CCC→CAC) eptA_2Phosphoethanolamine transferase EptA
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base T (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

GGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGCCATAA  >  KT2440/5204872‑5205137
                                                                                                                                       |                                                                                                                                   
ggTTCTTCGCTACGCGGGCTTCGTAGTCCTTCCGCTTAATCCCC‑GATAAAATGCAAGGCACCTAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTgg                                                                                                                                     <  5:263000/135‑1 (MQ=255)
ggTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTgg                                                                                                                                     >  5:21307/1‑135 (MQ=255)
        ggtACGCGGGCAGCGTAGTCCTTGCGCTGC‑TGCCCGGAGAACTTGAAACGCACCGAACCCGCCGTTTCGGTCCCGCACGTGTGCAGGGCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTa                                                                                                                             <  7:183217/133‑1 (MQ=255)
         cTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCACGGCACCGAAACCGCCGTTTCCGTACCGCACGAGGGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTa                                                                                                                             <  8:194447/134‑1 (MQ=255)
          tACGCGGGCATCGTAGTCCTTGCGCGTAATGACC‑GCGAACATGCAAGGCACCCAAACCCCCCTTTCCCTTCCCCACGAGTGCACCTAGGAAAACGCCCGCCGGCCCTGCAACCAGGCCAGAGTGGGCGTGGGAta                                                                                                                           >  2:358238/1‑134 (MQ=255)
            cgcgGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGGCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATcgc                                                                                                                         <  6:113955/135‑1 (MQ=255)
             gcgGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGGACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATcgc                                                                                                                         <  7:264354/134‑1 (MQ=255)
              ggggCATCGTAGTCCTTGCGCTTCATGCAG‑GTGAACATGCATGGCACCGAAAGCGCGGTTTCCGTACCGCACGAGTGCACGTCGTAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCgg                                                                                                                       <  6:264760/134‑1 (MQ=255)
              cgGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCgg                                                                                                                       <  3:52285/135‑1 (MQ=255)
                    tCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCCCGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGGGGTATCGCGGTCATAg                                                                                                                 >  6:224550/1‑135 (MQ=255)
                      gTAGTCATTGCGATTCATCCCC‑TAGAACAGGCAAGGACCCTAACCCGCCGTTCCCTTACCGCACGAGGGCACGTCGCAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGcc                                                                                                               <  1:252841/135‑1 (MQ=255)
                       tAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTGCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGcc                                                                                                               <  5:202056/134‑1 (MQ=255)
                       tAGGCCTGGCGCTTCAGGCCC‑GAGAACATCCAAGCCACCGAACCCGCCGTTTCCGTCCCGCACGAGTGAACGTCGCAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAgac                                                                                                               <  6:21307/134‑3 (MQ=255)
                        aGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCa                                                                                                             <  4:100210/135‑1 (MQ=255)
                                                aaCATGCAAGGCACCGAAACCGCCGTTTCCGAACCGCACGAGGGCACGTCGGAAAACGCCAGCCGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGGCTGCCCgg                                                                                      <  6:312265/135‑1 (MQ=255)
                                                 aCATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCggg                                                                                     >  4:36596/1‑135 (MQ=255)
                                                 aCATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCggg                                                                                     >  8:153131/1‑135 (MQ=255)
                                                      cAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACtt                                                                                >  8:216671/1‑135 (MQ=255)
                                                        aGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTc                                                                              >  7:221684/1‑135 (MQ=255)
                                                                 aaCCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGaccac                                                                     >  8:216395/1‑135 (MQ=255)
                                                                 aaCCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGaccac                                                                     >  3:40965/1‑135 (MQ=255)
                                                                   cagccgTTTCCGTACCGCACGAGTGCACGTCGTAAAACGCCAGCAGGCCCTGCTCCCTGGCTAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGAccacca                                                                   <  5:320318/133‑1 (MQ=255)
                                                                        gTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACg                                                              >  6:314056/1‑135 (MQ=255)
                                                                                      cGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTgcg                                                <  8:359763/135‑1 (MQ=255)
                                                                                       gAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTgcgc                                               >  4:380748/1‑135 (MQ=255)
                                                                                           gCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCAt                                           <  1:101587/135‑1 (MQ=255)
                                                                                            cACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGCCCACCAGCACGGTCAGCGACTTGCGCTCATg                                          >  4:181501/1‑135 (MQ=255)
                                                                                            cACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATg                                          >  4:74150/1‑135 (MQ=255)
                                                                                                tCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGtttt                                      <  7:9004/135‑1 (MQ=255)
                                                                                                tCGCAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGtttt                                      <  1:358238/135‑1 (MQ=255)
                                                                                                    aaatCGCCAGCAGGCGCTGCTCCGTGGCGAGTTTGGTCGTGGTATCTCGGTCATTGCCCAACACGCCGATATGGTCTGCCCTGGCACTTTCCCAGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGcc                                  <  2:269475/135‑1 (MQ=255)
                                                                                                    aaaaCGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGcc                                  <  8:247186/135‑1 (MQ=255)
                                                                                                         gCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGcc                             >  8:132570/1‑135 (MQ=255)
                                                                                                          ccAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGGCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCg                            <  4:294438/135‑1 (MQ=255)
                                                                                                                      gCTCCCTGGCGAGATTGTGCTTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTcccc                                                                           <  4:332915/76‑1 (MQ=255)
                                                                                                                      gCTCCCTGGCGAGATTGTGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTcccc                                                                           <  4:332913/76‑1 (MQ=255)
                                                                                                                      gCTCCCTGGCGAGATTGTGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTcccc                                                                           >  3:332913/1‑76 (MQ=255)
                                                                                                                      gCTCCCTGGCGAGATTGTGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTcccc                                                                           >  3:332915/1‑76 (MQ=255)
                                                                                                                        tCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCa              >  7:188831/1‑135 (MQ=255)
                                                                                                                               cgtgaTTGGGGGTGGTATCGCGGTCATTGCCCAACACGCGGAAATGGGCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGc       <  6:343485/132‑1 (MQ=255)
                                                                                                                               cGAGATTGGGCGTGGTATCGCGGTCATAGCCCACCACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGc       >  3:248528/1‑135 (MQ=255)
                                                                                                                               cGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCTACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCg        >  3:114532/1‑134 (MQ=255)
                                                                                                                               cGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCg        >  1:33472/1‑134 (MQ=255)
                                                                                                                                gagaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGcc      <  5:224550/135‑1 (MQ=255)
                                                                                                                                gagaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGcc      >  5:115447/1‑135 (MQ=255)
                                                                                                                                  gaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGCCAt    <  4:40965/135‑1 (MQ=255)
                                                                                                                                    ttGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGCAATGGTCTCCCCGGGCACTTTCCCCGACCACCAGCACGCTCAGCGACTAGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGCCATaa  >  5:343485/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GGTTCTTCGCTACGCGGGCATCGTAGTCCTTGCGCTTCATGCCC‑GAGAACATGCAAGGCACCGAAACCGCCGTTTCCGTACCGCACGAGTGCACGTCGGAAAACGCCAGCAGGCCCTGCTCCCTGGCGAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAGGCCGCATCACGCTTGGCATCTTCGCCATAA  >  KT2440/5204872‑5205137

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