Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,831,753 T→G 27.9% T158P (ACG→CCG)  CJ019_01748 ← NAD(P)H dehydrogenase (quinone)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,831,7530TG27.9% 45.4 / 5.3 52T158P (ACG→CCG) CJ019_01748NAD(P)H dehydrogenase (quinone)
Reads supporting (aligned to +/- strand):  ref base T (16/20);  new base G (2/14);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCACGCGCCGCCTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCAGGCTGGCGGTG  >  KT2440/1831621‑1831883
                                                                                                                                    |                                                                                                                                   
aGCCACGCGCCGCCTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGGGGCGGGCGCCCTAAGCCGTgg                                                                                                                                   >  5:176178/1‑134 (MQ=255)
     cgcgCCGCCTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGccac                                                                                                                             <  7:309015/136‑1 (MQ=255)
      gcgcCGCCTCCAGGGCCTTGGCCGGGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTGGCCATCGGCGCCGGCGTGGGGGCTGGCGCCATAAGGCGTGCCGccac                                                                                                                             <  4:367096/135‑1 (MQ=255)
        gccgccTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGGGCCGccccccc                                                                                                                          >  8:76121/1‑134 (MQ=255)
          cgccTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCGGGCGGCCTAAGGGCTGccccccccccc                                                                                                                         >  6:393500/1‑132 (MQ=255)
               ccAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGtctc                                                                                                                   <  8:425301/136‑1 (MQ=255)
                cAGGGCGGTGGCCGGGGGCGCCCGGCGTGGGCCCCGGGCGGGGCACAGGGCGATTTCGGGGGGGGCGAGTGCGCGCTTGCAATCGGCGCCGGGGTGGTGGCTGGCGCCATACGGCGTGCCGCCACCACGGGtctc                                                                                                                   <  2:85464/135‑1 (MQ=255)
                cAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGtctc                                                                                                                   <  3:137112/135‑1 (MQ=255)
                 aGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCTATTTCGTGTGGGTCTAGTTCGCGCTTGCCATCGGCGCGTGGGTGGTGGCTGGCGCAATAAGGCGTGCCGCCACCACGGGtctc                                                                                                                   <  6:176136/134‑1 (MQ=255)
                 aGGGCCGTGGCCGTGGTCGCCAGGCGTGGGCCCCGGGCGCGTCACAGGGCGATTTCGTGTGGGGGGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATACGGCGTGCCGCCACCACGGGtctc                                                                                                                   <  6:129238/134‑1 (MQ=255)
                          gCCGTGGTCGCCAGGCGTGGGCCCAGGGCGCGGCACAGGGGGATTTCGTGTGGGGCGAGTGCGCGCTTGCCAGCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCg                                                                                                        <  8:362171/136‑1 (MQ=255)
                          gCCGGGGGCGCCTGGCGTGGGCCCTGGGCGGGGCACAGGGCTATTGCGGGTGGGGCGAGTTCGCGCTTGCCATCGGCGCCGGGGGGGGGGCGGGCGCCATACGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCg                                                                                                        <  8:82388/136‑1 (MQ=255)
                           cggTGGGCGCCAGGCGTGGGCCCTGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCg                                                                                                        <  3:78876/133‑1 (MQ=255)
                           ccGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCCCCCCGGGTCTCGCGCAGTGGCg                                                                                                        >  1:15189/1‑135 (MQ=255)
                           ccGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCCCCCCGGGTCGCGCGCCGTGCCg                                                                                                        >  6:131405/1‑135 (MQ=255)
                           ccGTGGTCGCCAGGCGTGGGCCCAGGGCGCGGCACAGGGCGATCTCGTGGGGGTCGAGTGCGCGCTTGCCATCGGCGCCTGCGTGGTGGCGGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCg                                                                                                        <  6:352737/135‑1 (MQ=255)
                               ggTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCg                                                                                                   <  8:131304/136‑1 (MQ=255)
                                ggcGCCAGGCGTGGGCCCAGGGCGCGGCACAGGGCGATTTCGTGGGGGGCGAGTTCGCGCTGGCCAGCGGCGCCTGCGTGGGGGCGGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCg                                                                                                   <  1:271506/133‑1 (MQ=255)
                                 tCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCg                                                                                                   <  2:32335/134‑1 (MQ=255)
                                  cGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCCTAAGGCGTTCCGCCCCCACGGGTCTCGAGCAGTGCCGCCgcgcg                                                                                                 >  4:78876/1‑134 (MQ=255)
                                  cGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCCCCCCGGGGCTCGAGCAGTGCCGCCgcgcgg                                                                                                >  7:362171/1‑134 (MQ=255)
                                      aGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCAGCGGCGCCTGCGTGGGGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAgg                                                                                            <  7:288770/136‑1 (MQ=255)
                                       ggCGTTGGCCCAGGGCGCGGCAGAGGGCGATTTCGTGGGGGGCGAGTTCGCGCTTGCCAGCGGCGCCTGCGGGGGGGCGGGCGCCATACGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAgg                                                                                            <  5:200403/135‑1 (MQ=255)
                                       ggCGTGGGCCCCGGGCGCGGCACAGGGCGATTGCGGGTGGGTCGCGTGCGCGCTGGCAAGCGGCGCCGGCGGGGGGGCTGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAgg                                                                                            <  3:237848/135‑1 (MQ=255)
                                        gggtgggCCCAGGGGGCGGCATTGGGCGTTTTGGTGGGGGTGTTGTTGGCGCTGGCCAGGGGCGCTGGCGTGGGGGGTGGCGCATTATGGGGGGCCGCCACCACGGGTCTCGCGCAGTGCCGTCTCGCTGTAgg                                                                                            <  8:338639/129‑1 (MQ=255)
                                          gtgggCCCAGGGCGCGGCACAGGGCGATTTCGTGGGGGGCGAGTTCGCGCTTGCCAGCGGCGCCGGCGGGGGGGCGGGCGCCATAAGGCGGCCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑a                                                                                        <  4:314544/132‑1 (MQ=255)
                                          gTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑a                                                                                        <  3:188518/135‑1 (MQ=255)
                                              gCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGGTGGCGGCCTAAGGCGTGCCCCCCCCCCGGGGCGCGGGCAGGGCCGCCGCGCTGTAGGGC‑GGGcc                                                                                    >  1:371669/1‑135 (MQ=255)
                                              gCCCAGGGCGCGGCACAGGGCGATTTCGGGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGGGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGcc                                                                                    <  3:327491/135‑1 (MQ=255)
                                                ccAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCGGGCGCCATAAGGCGTGCCCCCCCCCCCGGTGTAGCGCAGTGGCGCAGGGCTGTAGGGC‑CGGCCCa                                                                                  >  5:214654/1‑135 (MQ=255)
                                                    ggggcgGCACAGGGCGTTTTCGGGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGGGCCGCAACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCa                                                                               <  7:133767/131‑1 (MQ=255)
                                                           cacaGGGCGATTTCGTGTGGGGCGAGTTCGCGCTTGCCAGCGGCGCCTGCGTGGGGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATg                                                                       <  1:350118/135‑1 (MQ=255)
                                                                   gATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCCTAAGGCGTGCCGCCACCACGGGTCTCGCGCAGTGCCGCCTCGCTGTAGGGC‑CGGCCCCTCACCCGCAAGCCgtggtg                                                               >  5:211042/1‑135 (MQ=255)
                                                                   gATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCCTAAAGCGTGCCGCCCCCACGGGTCTCGGGCCGTGGCGACTCGCTGGAGGGC‑AGGCCCCACACCCGCATGCCgtggtg                                                               >  6:195754/1‑135 (MQ=255)
                                                                   gATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCGGGCGCCATAAGGCGTGCCGCCCCCACGGGGCTCGAGCAGGGCCGCCGCGCTGTAGGGC‑CGGGCCATCACCAGCCTGCCggggtg                                                               >  6:442926/1‑135 (MQ=255)
                                                                    tttgcGTGGGGGGCGAGTGCGCGGTTGCAAGCGCGGCCTGGGGGGGGGCGGGCCCAATACGGCGGGCCGCCACCCCGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGc                                                              <  3:299963/131‑1 (MQ=255)
                                                                        cGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCCTAAGGCGTGCCGCCCCCCCCGGTCTCGCGCAGTGCCGACTCGCTGTAGGGC‑ACGCCCAAAACCCCCCAGCCGGGGTGCATCa                                                          >  3:322239/1‑135 (MQ=255)
                                                                        cGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAATGCCTGCCGCCACCACGGGTCTCGCGCAGTGCCGCCTCGCGGTACGGC‑AGGCCCATCACCAGCCTGCCGTGGTTCATcc                                                          >  6:187074/1‑134 (MQ=255)
                                                                        cGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGGGGTGGCTGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAg                                                         <  8:124895/136‑1 (MQ=255)
                                                                          tgtgGGGCGAGTTCGCGCTTGCCAGCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCg                                                       <  8:119786/136‑1 (MQ=255)
                                                                             gggTCGAGTTCGCGCTTGCCATCGGCGCCTGCCTGGTGGCAGGCCCCAAAAGGCGGGGCCCCACACCGGGACTAGACAAGTTCCGAAACGCGTTAGGGCAAGCCCCAACA‑CACCCTGCCGTGGTGCATAAACGgca                                                    >  8:413558/1‑136 (MQ=255)
                                                                                     tgcgcgCTTGCCAGCGGCGCCGGCGGGGTGGCTGGCGCCATACGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGcatcat                                             <  7:327813/133‑1 (MQ=255)
                                                                                        gcgcTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGa                                          <  1:19027/135‑1 (MQ=255)
                                                                                         cgcTGGCCATCGGCGCCGGCGTGGTGGCTGGCGCCATACGGCGGGCCGCCCCCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGAca                                        <  7:76121/136‑1 (MQ=255)
                                                                                          ggttGCCATCGGCGCCGGGGTGGTGGCTGGCGCCATATGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACGAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGAca                                        <  5:172184/133‑1 (MQ=255)
                                                                                                     gcgcCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGttt                            <  7:378213/136‑1 (MQ=255)
                                                                                                            cGTGGTGGCGGGCGCCATACGGCGTGCCGCCACCACGGGTCTCGAGCAGGGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGgcc                     <  8:94384/136‑1 (MQ=255)
                                                                                                                gtggCGGGCGCCATAAGGCGTGCCGCCACCACGGGGCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGgccgcc                  <  6:201311/135‑1 (MQ=255)
                                                                                                                          ccATAAGGCGGGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAACATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCAGGCTgg       <  7:114809/136‑1 (MQ=255)
                                                                                                                              aaGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCAGGCTGGCgg    >  2:229791/1‑135 (MQ=255)
                                                                                                                                ggCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCCGGCTGGCGGTg  >  5:355383/1‑135 (MQ=255)
                                                                                                                                ggCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCAGGCTGGCGGTg  >  8:30659/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                   
AGCCACGCGCCGCCTCCAGGGCCTTGGCCGTGGTCGCCAGGCGTTGGCCCAGGGCGCGGCACAGGGCGATTTCGTGTGGGTCGAGTTCGCGCTTGCCATCGGCGCCTGCGTGGTGGCTGGCGCCATAAGGCGTGCCGCCACCACGGGTCTCGAGCAGTGCCGACTCGCTGTAGGGC‑AGGCCCATCACCAGCATGCCGTGGTGCATCAGCGGCAGCATCATCGACAGCAGGGTGGTTTCCTGGCCGCCGTGCAGGCTGGCGGTG  >  KT2440/1831621‑1831883

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: ‑