Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,625,302 T→C 8.3% S334P (TCC→CCC)  lepA → Elongation factor 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,625,3020TC 8.3% 140.6 / 4.6 48S334P (TCC→CCC) lepAElongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base C (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

CAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGATCTACGAACT  >  KT2440/1625168‑1625432
                                                                                                                                      |                                                                                                                                  
cAGGCTTAAAGAAGATTCAGCCGCAGGTCTAGGCCGGCATTTTCCCGGTCAGCTCCGACGACTTCTAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCt                                                                                                                                    <  6:257952/135‑1 (MQ=255)
    cTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGa                                                                                                                                >  5:276908/1‑135 (MQ=255)
                 cAGCCTCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTCCCCCTGAACGACTCATCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCtta                                                                                                                   >  2:326141/1‑134 (MQ=255)
                            cTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTg                                                                                                        >  1:301710/1‑135 (MQ=255)
                                       tGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTgg                                                                                             >  5:6477/1‑135 (MQ=255)
                                          tCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGc                                                                                           <  8:187275/134‑1 (MQ=255)
                                                 cAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTGTGGCTTCTGTTGCGGTTTCCTGGGCATGTGGca                                                                                   >  8:48355/1‑135 (MQ=255)
                                                      ccgccgACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATgg                                                                              <  5:64089/135‑1 (MQ=255)
                                                        gaccacTGCGTGGATTTCCGCGACGCGCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATgg                                                                              <  2:20607/133‑1 (MQ=255)
                                                                 gaggagTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAg                                                                   <  4:248064/135‑1 (MQ=255)
                                                                 gAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATcc                                                                     <  4:162856/133‑1 (MQ=255)
                                                                  aGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAg                                                                   <  6:276908/134‑1 (MQ=255)
                                                                            cgacgaTCTTCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGa                                                        <  3:124352/135‑1 (MQ=255)
                                                                            cgACGCTCTGCAGAAGCGGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGa                                                        <  7:87244/135‑1 (MQ=255)
                                                                                       aGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATg                                             <  8:306050/135‑1 (MQ=255)
                                                                                            cTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATgacctg                                        <  2:91994/135‑1 (MQ=255)
                                                                                            cTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATgacct                                         >  5:205608/1‑134 (MQ=255)
                                                                                              gACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATgg                                                                              <  6:27402/95‑1 (MQ=255)
                                                                                              gACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATgg                                                                              >  5:27402/1‑95 (MQ=255)
                                                                                                 cgtGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATgacctggacct                                   <  6:315634/133‑1 (MQ=255)
                                                                                                 cgtGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATgacctggacct                                   <  1:179506/133‑1 (MQ=255)
                                                                                                 ccTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATgacctggacc                                    >  5:130844/1‑134 (MQ=255)
                                                                                                   tGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGa                                 <  3:305919/135‑1 (MQ=255)
                                                                                                   tGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGa                                 <  2:265391/135‑1 (MQ=255)
                                                                                                   tGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGa                                 <  2:182035/135‑1 (MQ=255)
                                                                                                       cGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATcac                             <  8:150256/135‑1 (MQ=255)
                                                                                                        gACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATcacc                            >  4:268805/1‑135 (MQ=255)
                                                                                                        gACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATcaca                            >  1:61402/1‑134 (MQ=255)
                                                                                                            cctcGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCg                        >  4:332769/1‑135 (MQ=255)
                                                                                                             ctcGCTGCAGTACATGCCGGAAAGCCCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGa                                 <  2:58481/125‑1 (MQ=255)
                                                                                                             ctcGCTGCAGTACATGCCGGAAAGCCCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGa                                 >  1:58481/1‑125 (MQ=255)
                                                                                                              tcGCTGCAGTACCTGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGGTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCg                        <  1:17160/133‑1 (MQ=255)
                                                                                                              tcGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCg                        <  6:414/133‑1 (MQ=255)
                                                                                                                        aCATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGa            <  1:175702/135‑1 (MQ=255)
                                                                                                                        aCATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGa            >  6:47133/1‑135 (MQ=255)
                                                                                                                        aCATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGa            <  4:26511/135‑1 (MQ=255)
                                                                                                                        aCATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGa            <  2:301710/135‑1 (MQ=255)
                                                                                                                          aTGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGAt           >  6:69908/1‑134 (MQ=255)
                                                                                                                          aTGCCGGAAAGCCCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGa                                                  <  2:4307/95‑1 (MQ=255)
                                                                                                                          aTGCCGGAAAGCCCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGa                                                  >  1:4307/1‑95 (MQ=255)
                                                                                                                               ggAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGATCTACGa     >  7:248799/1‑135 (MQ=255)
                                                                                                                               ggAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGATCTACGa     >  1:329418/1‑135 (MQ=255)
                                                                                                                               ggAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGAACACCACCGCACCCAGCGTGATCTACGa     >  1:282340/1‑135 (MQ=255)
                                                                                                                               ggAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGAGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGAGTGATCTACGa     >  8:175141/1‑135 (MQ=255)
                                                                                                                                 aaaGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATg                                                                                        >  2:221173/1‑50 (MQ=255)
                                                                                                                                 aaaGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATg                                                                                        <  1:221173/50‑1 (MQ=255)
                                                                                                                                    gCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATcc                                                                     >  4:183707/1‑66 (MQ=255)
                                                                                                                                    gCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATcc                                                                     <  3:183707/66‑1 (MQ=255)
                                                                                                                                    gCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGATCTGGACCTGATCACCACCGCACCCAGCGTGATCTACGAACt  >  1:277526/1‑133 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCTCCGACGCGCTGGGCTTTGGCTTCCGTTGCGGTTTCCTGGGCATGTTGCACATGGAGATCATCCAGGAGCGCCTGGAGCGCGAATATGACCTGGACCTGATCACCACCGCACCCAGCGTGATCTACGAACT  >  KT2440/1625168‑1625432

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