Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,669,107 T→C 7.3% K27R (AAG→AGG)  yceI ← Protein YceI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,669,1070TC 7.3% 172.2 / 3.6 55K27R (AAG→AGG) yceIProtein YceI
Reads supporting (aligned to +/- strand):  ref base T (19/32);  new base C (2/2);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CTGGCTTCAGGCTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGTTGAGAGTCAAGCTTTGCGGCCCAGC  >  KT2440/5668978‑5669237
                                                                                                                                 |                                                                                                                                  
cTGGCTTCAGGCTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCtt                                                                                                                                   <  7:331537/131‑1 (MQ=255)
     ttCAGGCTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGc                                                                                                                          >  6:61068/1‑135 (MQ=255)
           cTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAt                                                                                                                    <  3:14522/135‑1 (MQ=255)
           cTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCATTCCTTGTCGATCTTGTCCTCGGCAGCCAt                                                                                                                    <  5:214145/135‑1 (MQ=255)
                cACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcct                                                                                                               <  2:171377/135‑1 (MQ=255)
                cACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcct                                                                                                               >  5:67926/1‑135 (MQ=255)
                cACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCCTGTACTCGGCAGCCAtggcct                                                                                                               >  2:121261/1‑135 (MQ=255)
                cACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTGCCTTGTCGATCCTGTACTCGGCACCCAtggcct                                                                                                               >  2:126646/1‑135 (MQ=255)
                cACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGGGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcct                                                                                                               <  4:102427/135‑1 (MQ=255)
                 aCTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcctg                                                                                                              >  3:198825/1‑135 (MQ=255)
                     tCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          <  6:319169/135‑1 (MQ=255)
                      cccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCcag                                                                                                         <  2:10482/135‑1 (MQ=255)
                      cccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          <  5:126198/134‑1 (MQ=255)
                      cccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          <  5:61068/134‑1 (MQ=255)
                      cccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGGGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCcag                                                                                                         <  2:17606/135‑1 (MQ=255)
                      cccAGCTGAAGTTGCCGGCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          <  3:315138/134‑1 (MQ=255)
                       ccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTATCCCAGGTTGATGATTTTCCAGTCAACGAACTCGTCCTGGCCATCCTTGTCGATGTTGTACTCGGCAGCCATGGCCTGGCcagc                                                                                                        >  6:214145/1‑135 (MQ=255)
                       ccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          >  4:175350/1‑133 (MQ=255)
                       ccAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          >  4:341034/1‑133 (MQ=255)
                        cAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTAGTCGGCAGCCATGGCCTGGCca                                                                                                          >  2:345239/1‑132 (MQ=255)
                        cAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCca                                                                                                          >  4:221891/1‑132 (MQ=255)
                          gCTGAAGTTGCCGTCGAACTCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGGGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  2:206657/135‑1 (MQ=255)
                          gCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGCTTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  5:174375/135‑1 (MQ=255)
                          gCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  6:63135/135‑1 (MQ=255)
                          gCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  1:115035/135‑1 (MQ=255)
                          gCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGGCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  6:352339/135‑1 (MQ=255)
                          gCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGGGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcaga                                                                                                     <  8:318505/135‑1 (MQ=255)
                            tGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcagaga                                                                                                   <  3:221891/135‑1 (MQ=255)
                             gAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcct                                                                                                               >  4:251474/1‑122 (MQ=255)
                             gAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCAtggcct                                                                                                               <  3:251474/122‑1 (MQ=255)
                             gAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcagaga                                                                                                   >  3:312589/1‑134 (MQ=255)
                                gTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCagcagagagc                                                                                                 >  6:376783/1‑133 (MQ=255)
                                                ttGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAgcgc                                                                                <  1:160319/134‑1 (MQ=255)
                                                   aaaGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCa                                                                              <  6:194381/133‑1 (MQ=255)
                                                    aaGTACAGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGATAGCAGAGCGGTACCGAGCGCCagag                                                                           >  3:321819/1‑135 (MQ=255)
                                                       ttccGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCgg                                                                        <  3:294736/133‑1 (MQ=255)
                                                           gTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCACTCAACGAACGCGTGCTGCCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCACACAGCAGAGCGGTACCGAGCGCCAGAGCgccc                                                                      >  1:164607/1‑130 (MQ=255)
                                                              gATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGt                                                                 <  7:202772/135‑1 (MQ=255)
                                                                   aGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCtttt                                                            <  4:198825/135‑1 (MQ=255)
                                                                    gccgTAGCCCAGGTCGCGGATTTTCCAGTCAACGAACGCGGCCTGGCCGTGCTTGTCGAGCTTGTACTCGGGAGCCTTGGCCTGTCCAGCAGAGAGCAGAGCGTTACCGAGCGCCAGAGCGGGAAAAGTCttttt                                                           <  5:94766/132‑1 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCATAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAATTCGTTTTCAACATCCTGttcc                                             >  6:319753/1‑135 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcc                                             <  6:270681/135‑1 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcc                                             >  1:167507/1‑135 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcc                                             <  6:351142/135‑1 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcc                                             >  3:127571/1‑135 (MQ=255)
                                                                                  ggCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcc                                             >  5:365908/1‑135 (MQ=255)
                                                                                    cTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGttcctt                                           <  6:124709/135‑1 (MQ=255)
                                                                                       aTTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTtccttcct                                        >  1:277160/1‑135 (MQ=255)
                                                                                             cAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACt                                                                                                                              >  4:104212/1‑43 (MQ=255)
                                                                                             cAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACt                                                                                                                              <  3:104212/43‑1 (MQ=255)
                                                                                              aGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTaa                                  <  2:167507/134‑1 (MQ=255)
                                                                                                     gaacgCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGtt                          <  5:376783/135‑1 (MQ=255)
                                                                                                         gcgTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGAGTAATTGAGGTTgaga                      <  5:79735/135‑1 (MQ=255)
                                                                                                                            cGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGTTGAGAGTCAAGCTTTGCGGCCCa    <  4:127571/134‑1 (MQ=255)
                                                                                                                              aTCCTGTTCTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGTTGAGAGTCAAGCTTTGCGGCCCAGc  <  1:126646/134‑1 (MQ=255)
                                                                                                                              aTCCTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGTTGAGAGTCAAGCTTTGCGGCCCAGc  <  1:121261/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
CTGGCTTCAGGCTTGGCACTGTCCCAGCTGAAGTTGCCGTCGAAATCCTTGAAAGTACCGTGGATGAAGCTGTAGCCCAGGTGGCTGATTTTCCAGTCAACGAACGCGTGCTGGCCTTCCTTGTCGATCTTGTACTCGGCAGCCATGGCCTGGCCAGCAGAGAGCAGAGCGGTACCGAGCGCCAGAGCGGCAAAAGTCTTTTTCAACATCCTGTTCCTTCCTGTGTAATTGAGGTTGAGAGTCAAGCTTTGCGGCCCAGC  >  KT2440/5668978‑5669237

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