Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 305,368 C→G 25.2% K19N (AAG→AAC ydeP_1 ← Protein YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440305,3680CG25.2% 75.0 / 6.7 52K19N (AAG→AACydeP_1Protein YdeP
Reads supporting (aligned to +/- strand):  ref base C (14/23);  new base G (5/8);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGATAAGGTTGAGGTT  >  KT2440/305235‑305498
                                                                                                                                     |                                                                                                                                  
cGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGcg                                                                                                                                   >  5:150220/1‑134 (MQ=255)
   ggCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTa                                                                                                                               <  7:251279/136‑1 (MQ=255)
     cAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAg                                                                                                                              <  1:254579/135‑1 (MQ=255)
          gAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTAAGCGCATTAACCCTTCCGACCCCAGCCGTGGCCCCGCTACGCAGTGCGCCCCCACCGCCGGCGGGGCCCTGGTAGGGGTGGTTGGGAggg                                                                                                                        >  7:182017/1‑136 (MQ=25)
          gAAACCGCCGTTCTGGTTGGGCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAgg                                                                                                                         <  6:65169/135‑1 (MQ=255)
            acccgccgTTCTGGGTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAggggg                                                                                                                      <  7:394530/134‑1 (MQ=255)
            aaCCGCCGTTCTGGTTGGTCTTGAGCTGGGCGCGAATGTTCTTCAGCGCATTATCGCGGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCGGTGGTAGGGCTTGTAGCGAgggg                                                                                                                       <  4:407071/135‑1 (MQ=255)
                gccgTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGGAGGGCGTGTAGCGAgggggggg                                                                                                                   >  5:33995/1‑134 (MQ=255)
                     tCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCATCGCATTATCGCTGCCGACCCAGGCCTTGGCCACCCTACGCAGTGCGCCCCAACCCCCGGCCGGGCCGTGGTAGGGGATGTAGCGAGGGGGGGAGGCTg                                                                                                              >  6:36542/1‑135 (MQ=255)
                          ttggtCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTggg                                                                                                            <  4:89116/132‑1 (MQ=255)
                                ttGAGCCGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTCCGCAGTGCGCCCCAACCGCGGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTcc                                                                                                   <  6:359732/135‑1 (MQ=255)
                                           gcgAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGttgttg                                                                                        >  6:137683/1‑135 (MQ=255)
                                           gcgAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGttgttg                                                                                        >  5:172639/1‑135 (MQ=255)
                                           gccatggTTCTTCAGCTCTTTTTGGCTGCGCACCCTGGCCTGGGCCTCTCTACGCCTGGCGCCCCACGCGCCGGCGGGGGGGTGGTAGTGATTGTCGCTTGGGGGGGAGGATGGGGTGTTGTCCGGAAGTTGTTgg                                                                                       <  8:182017/130‑1 (MQ=255)
                                                  ttcttcAGCGCATTCTCGCGGCCGCCCCTGGCGTGGGCCACGCTCCGCAGTGCGCCCCACCCGCCGGGGGGGCGGGGGTTGGGCGGGTTGCGAGGGGGGGAGGCTGGGGTGGTGTCCGGAAGTTGTTGGTACGa                                                                                  <  6:371904/134‑1 (MQ=255)
                                                     ttcAGCGCATTAACGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCCCCCCAACCCCCCGCGGGGCCGGGGGAGGGCGGGTTGAGCGGGGGGGAGGCGGGGGGGTTGTGGGGTAGTTTTTTGTACGAGGTc                                                                              >  2:191106/1‑135 (MQ=255)
                                                       cAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCAc                                                                            >  3:118154/1‑135 (MQ=255)
                                                             aTTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTg                                                                       >  1:88502/1‑134 (MQ=255)
                                                             aTTATCCCTGCGCACCACGGCGTTCGCCACGCACCGCCGGGCCCCCCTCCCGGCGGGGTGGCTGTGGAAGGCGTTGTAGCGAGGGGGGGAGGCTTGGGTGTTGCCCGGATGTTGTTGGAACGAGGCCACGACTGt                                                                      <  4:14039/135‑1 (MQ=255)
                                                               tATCGCTGCCGACCCAGGCCTGGGCCACGCTACGCAGGGCGCCCCACCCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTa                                                                    <  3:435825/135‑1 (MQ=255)
                                                                   gcggCCGACCCCGGCCTGGGCCCCGCTCCGCGGTGCCCCCCCCCCGGGGGGGGGGGGGGGTTTGGGGTGGTAGCGCGGGGGGGAGGCGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGa                                                                 <  1:212222/131‑1 (MQ=255)
                                                                    cTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGGTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACt                                                               <  4:377780/135‑1 (MQ=255)
                                                                       ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTgg                                                                                       >  3:378590/1‑108 (MQ=255)
                                                                       ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTgg                                                                                       <  4:378590/108‑1 (MQ=255)
                                                                       ccGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgc                                                           >  7:326159/1‑136 (MQ=255)
                                                                         cccccGGGCGTGGGCCCCGCTACGCCGGGCGCCCCGCCGGCCGGGGGGGCGGTGGTTGGGCTTGTAGGGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                                         <  7:392391/134‑1 (MQ=255)
                                                                          cccccGGCCTGGGCCCCGCCCCGCGGGGCGCCCCCACCGCGGGGGGGGCCGGGGGAGGGCGGGCAGCGCGGGGGGGCGGGGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGGCACGACTGTTACGACTccgccg                                                         <  6:172639/134‑1 (MQ=255)
                                                                          aCCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCGGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                                         <  4:184483/135‑1 (MQ=255)
                                                                          aCCCAGGCCTGGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTccgccg                                                         <  2:295709/135‑1 (MQ=255)
                                                                                  cTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTGCGGAAGTTGGTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTg                                                 >  5:189546/1‑135 (MQ=255)
                                                                                    tGGCCACGCACCGCTGGGCGCCCCCCCCGCCGGGGGGGCGGGGGGAGGGCTGGGCGCGCGGGGGGGAGGGGGGGGGGGTGTCCGGAAGGTGTTGTTACGAGGTCACGACTGGTACGACTCCGCCGCTGGGCTGTa                                               <  5:248085/135‑1 (MQ=255)
                                                                                        cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGGGTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTTGTACGACGTCACCCCTGTTACGACACCGCCGCGGGGCGGTAGccc                                           >  5:301121/1‑135 (MQ=255)
                                                                                        cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGAccc                                          >  8:120416/1‑136 (MQ=255)
                                                                                        cACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGAcc                                           >  5:171833/1‑135 (MQ=255)
                                                                                             gaCGCCGTGCGCCCCAACCCGCGGGGGGGGCGTGGTTGGGCGTGTAGCGCGGGGGGGCGGCGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGGCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                                     <  8:103597/135‑1 (MQ=255)
                                                                                               cgctgggcgcCCCCCCCGCCGGGGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTg                                    <  1:381967/129‑1 (MQ=255)
                                                                                               cGCAGGGCGCCCCAACCGCCGGCGGGGCCGTGGGTGGGCTTGTCGCGTGGGGGGGTGGCGGGGGTGTTGTCCGGACGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                                     <  8:86514/134‑1 (MQ=255)
                                                                                                gccgggcgcCCCACCCGCCGGGGGGGCGGGGGTCGGGCTGGCCGCGAGGGGGGGAGGGGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                                     <  3:247135/128‑1 (MQ=255)
                                                                                                gccgTGCGCCCCCCCCGCCGGGGGGGCGGTGGTAGGGGTGGTCGCGAGGGGGGGAGGGTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGt                                     <  5:256325/130‑1 (MQ=255)
                                                                                                  aGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTgcgc                                 >  3:3330/1‑135 (MQ=255)
                                                                                                    ggcgcCCCCCCCCCGGGCGGGGCGGGGGGGGGGGTGGTCGGGCGGGGGGGCGGGGGGGGGGGTGGCCGGAAGTTGTTGGTCCGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCtg                               <  6:422311/134‑1 (MQ=255)
                                                                                                     ccgcCCCACCCGCGGGGGGGGCGGGGGTTGGGGTTTTAGGGTGGGGGGGAGGCTGGGGTGTTGTCCGGATGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCtgt                              <  5:85349/134‑1 (MQ=255)
                                                                                                           cAACCGCCGGGGGGGCGGTGGTAGGGCTTGTCGCGAGGGGGGGTGGCTGGGGTGTTGGCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTggg                       <  8:275298/136‑1 (MQ=255)
                                                                                                            aaCCGCCGGGGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTggg                       <  3:147473/135‑1 (MQ=255)
                                                                                                            aaCCGCCGGCGGGGCCGTGGTAGGGGTTGTAGCGTGGGGGGGAGGCTGGGGTGGTGTCCGGGAGTTGGGGGTAGGAGGTCCCGGCTGTTACGCCTCCGCCGCTGGGCTGTTGCCCCGCCGGGCGCTGTGCGTggg                       >  2:43511/1‑135 (MQ=255)
                                                                                                                gccgGCGGGGCCGGGGTAGGGCTGGTAGCGAGGGGGGGAGGCGGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGt                  <  8:291854/136‑1 (MQ=255)
                                                                                                                 cggggggggGCGGGGTAGGGGTTGTCGGGTGGGGGGGAGGGTGGGGTGTTGTCCGGATGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGt                  <  6:150220/130‑1 (MQ=255)
                                                                                                                 cgggcggggcggGGGTAGGGCTTGTAGGGAGGGGGGGAGGCGGGGGTGTTGGCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGt                  <  3:421330/133‑1 (MQ=255)
                                                                                                                  gggcggGGGGGTGGTTGGGTTTGTAGGTTGGGGGGGTGGCTGGGGTGTTGTCCGTAAGTTGTTGTTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCTGTGCGCTGTGCTTGGGCAAGTg                 <  5:188567/134‑1 (MQ=255)
                                                                                                                      ggggCCGTGGTAGGGGTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGa               >  1:375537/1‑133 (MQ=255)
                                                                                                                               gTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGATAaggttgaggt   >  7:30598/1‑136 (MQ=255)
                                                                                                                               gTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGATAaggttgagg    >  1:49401/1‑135 (MQ=255)
                                                                                                                                 aGGGCGTGTCGCGCGGGGGGGCGGCTGGGGGGTTGTCCGGAAGTTGTTGGTACGAGGGCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGATAaggttgaggtt  <  3:67500/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
CGGGGCAGTCGAAACCGCCGTTCTGGTTGGTCTTGAGCAGGGCGCGAATGTTCTTCAGCGCATTATCGCTGCCGACCCAGGCCTTGGCCACGCTACGCAGTGCGCCCCAACCGCCGGCGGGGCCGTGGTAGGGCTTGTAGCGAGGGGGGGAGGCTGGGGTGTTGTCCGGAAGTTGTTGGTACGAGGTCACGACTGTTACGACTCCGCCGCTGGGCTGTAGACCCGCGGTGCGCTGTGCTTGGGCAAGTGGATAAGGTTGAGGTT  >  KT2440/305235‑305498

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: