Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,975,241 T→G 43.3% L684R (CTG→CGG)  sasA_19 → Adaptive‑response sensory‑kinase SasA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,975,2410TG43.3% 12.3 / 2.4 28L684R (CTG→CGG) sasA_19Adaptive‑response sensory‑kinase SasA
Reads supporting (aligned to +/- strand):  ref base T (2/9);  new base G (5/4);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTC  >  KT2440/5975109‑5975353
                                                                                                                                    |                                                                                                                
aaGCCTTGATTGGCGCGGGCTTCCTGCAACCGGGCCAGAACGGCTGGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTgg                                                                                                                <  8:347200/135‑1 (MQ=255)
                         gCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTg                                                                                                                 >  4:316511/1‑109 (MQ=255)
                         gCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGGGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTg                                                                                                                 <  3:316511/109‑1 (MQ=255)
                          cAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATg                                                                                                           >  4:37575/1‑114 (MQ=255)
                          cAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATg                                                                                                           <  3:37575/114‑1 (MQ=255)
                               gggCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGGCCAGCGCCCCGCATGGCCCCGGCGacacca                                                                                 >  6:292043/1‑130 (MQ=255)
                                         ggCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGCCAATCGCCCGGCTTGGCCACGGCGAGCCCACCAGGggttg                                                                       >  5:62461/1‑135 (MQ=255)
                                                      cTCATGCCGGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCGGCATGCACGTATCGCCCGGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCgcc                                                           <  2:2866/134‑1 (MQ=255)
                                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCGTGCACGTATCGCCCCGGCCCAGCCCCAGGCGTGGGCCGGGGGGACCCCCGCGGGGTTGGTCGAGCTAGCCCCCCAGCGTCAAcc                                             >  6:291663/1‑134 (MQ=255)
                                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGGCCAGCGCCCCGCGTGGCCAGGGCGCACGCCCCCGGGGTTGTTGGAGCGCGCCCGCCCCCaacaacc                                             >  2:81180/1‑134 (MQ=255)
                                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGCCAGTCGCCAGGCTTGGCCAGGGCGCACGCACGAGGTGTTGTTGGCGCTCGCCAGCCCGCGACACcg                                             >  7:371872/1‑135 (MQ=255)
                                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAAGCCCAGCGCCCCGCGTGGGCCCGGGGCACGCCCGACGTGGTGTTGGCGCTCCCCACGCAGCAAAAAcc                                             >  6:355767/1‑134 (MQ=255)
                                                                    aTGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGCCCAGCGCCACGCGTGGGCCCGGCGCCCCCCCGCGGGGTTGTTGGAGCGCGCCCCCCCGCGACCAcgc                                            >  5:163419/1‑135 (MQ=255)
                                                                    aTGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGGCAAGCGCCACGCGTGGCCCGGGGGCACGCCAGAGGGGTTGTTCGCGCGCGCCCCCCCGCTACAAccc                                            >  2:280998/1‑133 (MQ=255)
                                                                           ccGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGCCCTGCGCCAGGCGTGGGCCGGGCGCCCGCCCGCGGTGTGGGTGGCGCTCGCCCGCCCGCGACCACCCGcccccc                                     >  7:303775/1‑131 (MQ=255)
                                                                              cTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCCTGCGCCCGGCTTGGCCAGGGCGCTCGCCCCCGGGGTTGTTGGCGCGCGCCCCCCCGCCACCACGCGCCAGccccgc                                 >  4:70361/1‑134 (MQ=255)
                                                                                         aCGCCCGTGCACCACCGCGCCTCGCGGCATGCACGTATCGCCCTGGCCTTGCGCCAGGCTGGGCCTGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCa                       <  7:88736/135‑1 (MQ=255)
                                                                                          cGCCGGTGCACCCGCCCGCCTCGCTGAATGCACGTAGCGCGCGGGCCATGCGCGTGGGTTGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGc                     <  3:352799/136‑1 (MQ=255)
                                                                                          cGCCCGGGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGc                     <  3:282420/136‑1 (MQ=255)
                                                                                            cccGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGGCCAGCGCCCCGCGTGGGCCGGGCGCACGCCCGCGGGGTTGGTGGCGCTCGCCCCCCCGCGCCACCGCGCCCCCCCCGGCGGGGTTggggcg                    >  6:135516/1‑130 (MQ=255)
                                                                                            cccGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGCCCTGCGCCCGGCGTGGGCCGGGGGGACGCCCGAGGTGGTGTTGGCGCGCGCCAGCCCGCGACCACCCGCCCGCCCCCGCGGGTTTTGAGCg                    >  6:107878/1‑135 (MQ=255)
                                                                                             ccGTGCCCCAGCGCGCCGCGCTGCAGGCACGTAGCGCCCGGCCCATGCGCCCGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCgg                   <  8:228724/135‑1 (MQ=255)
                                                                                                      cgcgcgcCTCGATGCAGGCACGTATCGCCCGGGCGATGAGGTTGGGTGGGCGTGGGCTAACGCACGAGTTGTTGTTCGAGCTCTCGAGCCAGCTACAACTCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTgag          <  8:253211/134‑3 (MQ=255)
                                                                                                        cgcgcCTCGCTGCAGGCCCGTAGCGCCCTGCCCAGGCGCCAGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGaa        <  6:258514/135‑1 (MQ=255)
                                                                                                        cgcgcCTCGCTGCAGGCACGACTCCCCCTGCCCAGGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGGTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGaa        <  6:253785/135‑1 (MQ=255)
                                                                                                          cgcCTCGCTGCATGCACGTATCGCCCGGGCCCTGCGCCCGGGGTTGCGCCGGCGCACGCCCCCGGGGTTTTTGGCGCGGGCCCGGCCGCGACCACGCGCCACCCCCGGCGGGTGTGGAGCGGGTGTCGGCGCAGc      >  8:263282/1‑135 (MQ=255)
                                                                                                             cgcGCTGCATGCAGGTTTCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTc  <  3:179978/134‑1 (MQ=255)
                                                                                                              tCGCTGCTGGCAGGTATCGCCCTGGCCATGCGCCTGGCTTGGCCTGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTc  <  8:180772/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                
AAGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTC  >  KT2440/5975109‑5975353

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