Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,690,862 T→G 40.6% 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,8620TG40.6% 7.7 / 7.9 40L281R (CTG→CGG) yhaJ_2HTH‑type transcriptional regulator YhaJ
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base G (4/13);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CTACCTGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTGTGTGCGGGTTATGCTGT  >  KT2440/4690731‑4690996
                                                                                                                                   |                                                                                                                                      
cTACCTGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCCTTGGGCCTggg                                                                                                                                     >  5:72321/1‑135 (MQ=255)
cTACCTGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCGTCGACAACATGGCCAACCGCCTGCGCGTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGgc                                                                                                                                    >  3:172104/1‑136 (MQ=255)
 tACCTGCCGCAGCCATTGGTCGCACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCGGGgc                                                                                                                                    >  7:32416/1‑135 (MQ=255)
          cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGGTGGGGCGGGGcgcgcgcccc                                                                                                                           >  2:151326/1‑131 (MQ=255)
          cAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTGCGACAACATGGCCTACCGCCTGCGCTTGTGGGTAGACATCATCTGGAGCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGc                                                                                                                           >  2:293509/1‑135 (MQ=255)
           aGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCt                                                                                                                          >  7:288293/1‑135 (MQ=255)
              cttgggTCGACCGGTTGATTGCCGCGGGCGCACTGGCCGAGGTGCGCTCCGACAACTTGCCAACCCGCCTGCTGTGGTGGTTAGACATAATGTGGAGCACGCCGCGCCCGTTGGGCCTGGGCGCGCGACGCTAcc                                                                                                                       <  2:108196/131‑1 (MQ=255)
                ttGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGGTGGGCCGGGGCGCGCGACGcgcccgg                                                                                                                     >  2:96688/1‑129 (MQ=255)
                ttGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCCTTGGGGCGGGGCGCGCGCCGCTCCCggg                                                                                                                    >  4:207204/1‑136 (MQ=255)
                  ggTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGCCTGGGCGCGcgacccgaccgg                                                                                                                     >  8:196436/1‑131 (MQ=255)
                              gATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGGCGCGACCGGGAAGTCATTCGccc                                                                                                       >  2:5480/1‑134 (MQ=255)
                                        ggCTCACTGGCCGTGTTGGGCTTCTACAACTTGGCCACCCGCTTGCGCTTGTTGGTAGACATCATGTGGATCAAGTGGCGCCGGTGGGGCCGGGGGGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCg                                                                                               <  6:279886/133‑1 (MQ=255)
                                         gcgcAGGGGCCGTGGGGCGCTTCCAAAACATGCCCAACCGCCGGCGGTGGGGGGTAGACATCATCTGCAGCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCgaga                                                                                            <  8:20449/135‑1 (MQ=255)
                                           gcACTGGCCGTGGTGCGCTTCGACAACATGGCCAACGGCCTGCGCTTGGGGGTAGACATCATCTGTAGCCAGGCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAgg                                                                                         <  3:47729/136‑1 (MQ=255)
                                                ggcggtggTGGGCTTCGACAACCTGCCCACCCGCCTGGGGTGGGGGGTATCCATCATGTGGATCACGGCGCGCCCGGGGGGCGGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcc                                                                                     <  1:96688/129‑1 (MQ=255)
                                                ggcggcggTGCGCGTCGACAACAGGGCCAACGGCCGGCGCTTGGGGGTCGACATCAGCGGGATCACGGCGCGCCCGTTGGGCGTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGccc                                                                                    <  4:35447/130‑1 (MQ=255)
                                                ggCCGAGGGGCGCTGCGACAACATGGCCAACGGCCTGCGCTTGGGGGTAGACATCATCGGCAACAAGGCGCGCCCGTGGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcc                                                                                     <  8:361691/135‑1 (MQ=255)
                                                    gtggggcgcTTCCCCAACCTGCCCCACCGCCGGCGCTGGTGGGACGACAGCTGCGGCACCACGCCGCGCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcccccc                                                                                 <  2:109867/130‑1 (MQ=255)
                                                    gtggTGGGCTTCCAAAACAGGCCCACCCGCGGGCGTTTGTGGGTACACATCATCTGCACCACGCCCCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcccccc                                                                                 <  5:351911/133‑1 (MQ=255)
                                                    gAGGTGCGCTTCGACAACATGCCCACCCCCCTGCGCTGGTGGGTAGACATCAGCTGGATCACGGCGCGCCCGTGGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGcccccc                                                                                 <  7:114965/135‑1 (MQ=255)
                                                            ctccGACAACAGGGCCAACCGCCTGCGCTTGGGGGTAGACAACACCGGCAGCACGCCGCGCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTg                                                                         <  1:14421/132‑1 (MQ=255)
                                                            cTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGGCCGTTGGGCCTGGGCGCGCGGCCCGACCAGGCAGTCCTGCGCCAGCTGGTAGAGCTCGAGCCCCCCCGGGGCTgg                                                                        >  3:251156/1‑135 (MQ=255)
                                                            cTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGCCGCGACCGGGAAGTCATGCGCCAGCGGTTCGAGAAGGAGCCCCCCCGGGccgg                                                                         >  7:9580/1‑133 (MQ=255)
                                                                    acaGGGCCAACCGCCGGCGGTTGTGGGTAGACATAATCGTGATCACGGCGCGCCCGTGGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCtt                                                                <  4:319751/136‑1 (MQ=255)
                                                                     cctGGCCACCCGCGTGCGCTTGGGGGACCACAACTTCTGGAGCCAGGCGCGCCCTTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCtt                                                                <  5:300415/133‑1 (MQ=255)
                                                                             aCCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGTGGCAGGGCGCGCGCCCCTAAAAGTAAGTCATTCGCCAGCGGTTAGAAAAAGATCCCCCCAAGGCCTCATATACCGTGCAcgcga                                                        >  6:293056/1‑134 (MQ=255)
                                                                               cgcgTGCGCTGGTGGGCAGAAATCATGTGGAACACGCCGCGCCCGTGGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTcc                                                      <  2:130694/135‑1 (MQ=255)
                                                                               cGCCTGCGCTTGGGGGTACACAACATCGGGACCACGGCGCGCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTcc                                                      <  6:87809/135‑1 (MQ=255)
                                                                               cGCCGGCGCTGGTGGGTACACATCATGGGCAGCACGGCGCGCCCGTTGGCCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTcc                                                      <  5:197956/135‑1 (MQ=255)
                                                                                    gctgttGTGGGTAGACTTCATCTGTATCTAGGGGCGCCGGTTGGGCGTGGGCGCGCGACGCTACCTGGAAGTCATTCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCTAGGccc                                                <  4:91969/132‑1 (MQ=255)
                                                                                                acaCATCATCTGCATCATGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTc                                     <  6:187489/133‑1 (MQ=255)
                                                                                                  aCATCATGTGGACCTCGTCGCGCCCGTGGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCTCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGc                                   <  6:369882/135‑1 (MQ=255)
                                                                                                          tGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGCAGGAGCCCCCCACGGCCTGATATCCCTTGCACGCCTCCCGGCCCTGCCGGCGTTAGCCCCGccca                           >  2:118773/1‑134 (MQ=255)
                                                                                                            cagcacgTCGCGCCCGTTGGGCCTGGCCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCATGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCttt                         <  7:196436/129‑1 (MQ=255)
                                                                                                              ccatgccgcgCCGGTTGGGCCGGGCCCCGCGACGCTACCTGGAAGTCATGCGCCAGCTGGTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTgg                       <  1:76497/129‑1 (MQ=255)
                                                                                                              ccAAGCCGCGCCCGGTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTgg                       <  8:99329/134‑1 (MQ=255)
                                                                                                                  ggcgcgCCCGTTGGGCCGGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCtg                   <  8:160337/133‑1 (MQ=255)
                                                                                                                  gTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCtg                   <  8:288293/135‑1 (MQ=255)
                                                                                                                    cgcgCCCGTTTGGCCTGGGCGCGCGCCGCGACCGGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCGGATATCCCTTGCACGCGTCCAGGCCCAGCCGGCGTTCGCCCCGCCCGGTGGGCtggg                 >  2:227497/1‑133 (MQ=255)
                                                                                                                    cgcgCCCGTTGGGCCTGGGCGCGCGACCCTACCGGGAAGTCATTCGCCAGCAGTTAGAGAAGGCGCCCCCCCCGGCCGGATATCCCATGCGCGCGTCCAGGCCCTGCCGGCGTTCGCCCCGCCCGTTggggggtg                 >  6:153763/1‑130 (MQ=255)
                                                                                                                                   tGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTGTGTGCGGGTTATGCTGt  <  7:354396/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CTACCTGCCGCAGCCATTGGTCGAACCGCTGATTGCCTCGGGCGCACTGGCCGAGGTGCGCTTCGACAACATGGCCAACCGCCTGCGCTTGTGGGTAGACATCATCTGGATCAAGTCGCGCCCGTTGGGCCTGGGCGCGCGACGCTACCTGGAAGTCATGCGCCAGCTGTTCGAGAAGGAGCCCCCCAGGGCCTGATATCCCTTGCACGCGTCCAGGCCCTGCAGGCGTTCGCCCAGCCCTTTGGGCTGTGTGCGGGTTATGCTGT  >  KT2440/4690731‑4690996

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