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

TGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTG  >  KT2440/4690736‑4690979
                                                                                                                               |                                                                                                                     
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGGCTGGgcgcgc                                                                                                               >  5:322589/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGGCTGGgcgcgc                                                                                                               >  3:27076/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGgcgcgc                                                                                                               >  6:100986/1‑135 (MQ=255)
tGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGgcgcgc                                                                                                               >  4:235152/1‑135 (MQ=255)
 gCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGGCTGGGcgcgcg                                                                                                              >  1:236468/1‑135 (MQ=255)
 gCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGGTTGGCCTTGgcgccgg                                                                                                              >  3:353448/1‑132 (MQ=255)
   cGCAGCCATTGGGCGAACCGCTGATGGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGCCCAACCGCC‑GGCGCTTGTGGGTAGACATCATCTGGATCAAGGCGCGCCCGTTGGGCCTGGGcgcgcg                                                                                                              <  7:253352/133‑1 (MQ=255)
     cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGGTGGGGCGGGGCGCGCGgcgc                                                                                                          >  3:359260/1‑135 (MQ=255)
     cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGGTGGGCCGGGGCGCGcgccgc                                                                                                          >  6:10150/1‑135 (MQ=255)
                                      gctgtGGCGGCGGTGCGGTTGGAAAAC‑TGGGCAACCCCCGGGCGTTTGGGGGTACAATTAATCTGGAACACGCCGCGCCCGTGGGGCGGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAg                                                                         <  6:301120/131‑1 (MQ=255)
                                            gCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCGACCTGGAAGTCATGCGCCAGCGGTTCGCGAAGGCGccc                                                                   >  8:96617/1‑135 (MQ=255)
                                              cgtggTGCGCTTCGAAAACATGGCCAACCGCC‑TGCGCTGGTGGGTAGACATCATCTGGAGCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGccccc                                                                 <  3:106705/132‑1 (MQ=255)
                                              cgcggggcgcTCCGACAACATGGCCAACCGCC‑TGCGCTGGGGGTTCTACATCATCTGTATCAAGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGccccc                                                                 <  7:96617/129‑1 (MQ=255)
                                               gtggggcgcTTCTAAAACATGCCCAACCGCG‑GGCGGTTGGGGTTCCACATAAGGTGGATCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcccccc                                                                <  5:10150/130‑1 (MQ=255)
                                                      gcTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCt                                                         <  8:263442/135‑1 (MQ=255)
                                                                caTGCCCAACGGCC‑TGCGGTTGGGGGTTCAAATCATCTGCATCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCtt                                               <  8:367456/135‑1 (MQ=255)
                                                                   ggCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGCCAGGGCGCGCGGCGCCACCCGGAACGCCCGGGCCCGCCGGTAGACCAGGGGCCCCCCACGGGCAGCTATACCTTGcc                                            >  6:336592/1‑134 (MQ=255)
                                                                          cGCC‑TGCGCTGGGGGTCCAACACCTTCTGGATCACGGCCCGCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTcc                                     <  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)
                                                                                                                               |                                                                                                                     
TGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCC‑TGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTG  >  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: ‑