Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,690,862 T→G 46.0% L281R (CTG→CGG)  yhaJ_2 → HTH‑type transcriptional regulator YhaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,690,8620TG46.0% 3.0 / 2.7 27L281R (CTG→CGG) yhaJ_2HTH‑type transcriptional regulator YhaJ
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base G (2/11);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTG  >  KT2440/4690736‑4690979
                                                                                                                               |                                                                                                                     
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGGCTGGgcgcgc                                                                                                               >  5:322589/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGGCTGGgcgcgc                                                                                                               >  3:27076/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGgcgcgc                                                                                                               >  6:100986/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGgcgcgc                                                                                                               >  4:235152/1‑135 (MQ=255)
 gCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGGCTGGGcgcgcg                                                                                                              >  1:236468/1‑135 (MQ=255)
 gCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGGTTGGCCTTGgcgccgg                                                                                                              >  3:353448/1‑132 (MQ=255)
   cGCAGCCATTGGGCGAACCGCTGATGGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGCCCAACCGCCGGCGCTTGTGGGTAGACATCATCTGGATCAAGGCGCGCCCGTTGGGCCTGGGcgcgcg                                                                                                              <  7:253352/133‑1 (MQ=255)
     cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGGTGGGGCGGGGCGCGCGgcgc                                                                                                          >  3:359260/1‑135 (MQ=255)
     cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGGTGGGCCGGGGCGCGcgccgc                                                                                                          >  6:10150/1‑135 (MQ=255)
                                      gctgtGGCGGCGGTGCGGTTGGAAAACTGGGCAACCCCCGGGCGTTTGGGGGTACAATTAATCTGGAACACGCCGCGCCCGTGGGGCGGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAg                                                                         <  6:301120/131‑1 (MQ=255)
                                            gCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCGACCTGGAAGTCATGCGCCAGCGGTTCGCGAAGGCGccc                                                                   >  8:96617/1‑135 (MQ=255)
                                              cgtggTGCGCTTCGAAAACATGGCCAACCGCCTGCGCTGGTGGGTAGACATCATCTGGAGCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGccccc                                                                 <  3:106705/132‑1 (MQ=255)
                                              cgcggggcgcTCCGACAACATGGCCAACCGCCTGCGCTGGGGGTTCTACATCATCTGTATCAAGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGccccc                                                                 <  7:96617/129‑1 (MQ=255)
                                               gtggggcgcTTCTAAAACATGCCCAACCGCGGGCGGTTGGGGTTCCACATAAGGTGGATCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcccccc                                                                <  5:10150/130‑1 (MQ=255)
                                                      gcTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCt                                                         <  8:263442/135‑1 (MQ=255)
                                                                caTGCCCAACGGCCTGCGGTTGGGGGTTCAAATCATCTGCATCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCtt                                               <  8:367456/135‑1 (MQ=255)
                                                                   ggCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGCCAGGGCGCGCGGCGCCACCCGGAACGCCCGGGCCCGCCGGTAGACCAGGGGCCCCCCACGGGCAGCTATACCTTGcc                                            >  6:336592/1‑134 (MQ=255)
                                                                          cGCCTGCGCTGGGGGTCCAACACCTTCTGGATCACGGCCCGCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTcc                                     <  2:275642/135‑1 (MQ=255)
                                                                                         ggTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGCCGCGACCGGGAAGTAATGCGCCAGCGGTTGGAGAAGGAGCCCCCCCCGGCCAGAGATGCCATGCACGCCTCCAGGCCCAGCAGGcc                       >  1:281831/1‑134 (MQ=255)
                                                                                              acctcatcGGGATCACGGCCCGCCGGTTGGGCCGGGGCGCGCGACGCTACCTGGAACTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGGCCAGGCCCTGCAGGCGGTCGc                  <  7:20482/132‑1 (MQ=255)
                                                                                              aCATCATCTGGAGCACGCCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGc                  <  2:281831/135‑1 (MQ=255)
                                                                                                atcatcTGCATCACGGCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGc                  <  6:269285/133‑1 (MQ=255)
                                                                                                        gATCATGGCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCttt        <  6:322589/135‑1 (MQ=255)
                                                                                                         aTCCAGTCGCGCCCGGGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCCCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTg       <  5:299234/135‑1 (MQ=255)
                                                                                                          ccccggcgcgCCCGTGGGGCGGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTgg      <  1:330995/129‑1 (MQ=255)
                                                                                                              ggcgggcCGGTGGGGCGTGGGCGCGCGACGCTTCCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCtg  <  1:80033/130‑1 (MQ=255)
                                                                                                              ggcgcgCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCtg  <  6:41970/133‑1 (MQ=255)
                                                                                                                               |                                                                                                                     
TGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTG  >  KT2440/4690736‑4690979

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