Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,449,130 A→C 26.3% W97G (TGG→GGG)  CJ019_05126 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,449,1300AC26.3% 63.6 / 8.3 50W97G (TGG→GGG) CJ019_05126hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base C (12/2);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 5.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCACACGGGCGGTTTCCTGCACCACCAGCG  >  KT2440/5449002‑5449264
                                                                                                                                 |                                                                                                                                      
cacgcTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGACGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCAGGcccg                                                                                                                                  >  1:284699/3‑133 (MQ=255)
               gCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATGCGGTCCAGTGCCGTGGGCGATGGGGTGGAGCCCCCCCCCCCCGCCCCGGCCCCcgccgcccg                                                                                                                    >  1:333214/1‑129 (MQ=255)
                ccGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCgcgc                                                                                                                  <  4:212264/135‑1 (MQ=255)
                  gTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCg                                                                                                               <  7:313010/136‑1 (MQ=255)
                   tAATGCAACCCCCCACTGCTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGAGGCAAAGCGCATCCGTTCAAGTGCAGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTAGCGCCCg                                                                                                               <  1:421394/135‑1 (MQ=255)
                         aaCCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGcc                                                                                                          <  1:115256/134‑1 (MQ=255)
                         aaCCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGcc                                                                                                          <  2:101444/134‑1 (MQ=255)
                              accacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGGATGTCCCCCCCCCCCCCCGCCGCGGCCCCGCGCCCGCGCCCCGCGCCCacccc                                                                                                    >  2:437216/1‑131 (MQ=255)
                              accacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCCCCCCCCCCCCCGGCCCGGTCCCCGCGCTCGCGCCCGGCGccccccac                                                                                                    >  5:35180/1‑130 (MQ=255)
                              accacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCgccccgcccc                                                                                                   >  8:205146/1‑136 (MQ=255)
                              accacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCCCCCCCCCACGCCGGCTCCCCGCCCTCGCGCCCGGCGCCCAGca                                                                                                     >  2:272809/1‑134 (MQ=255)
                              accacTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCGGGTGACCGCCCCCCCCGGCCGGGTCCCCCCGCGcgcgccccgcgcccccccc                                                                                                    >  1:57041/1‑130 (MQ=255)
                                  cTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCAGGCCGGGTCCCCGCGCACCCGCCCCGCGCCCAGcaccac                                                                                                 >  4:202135/1‑134 (MQ=255)
                                     gTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCCCCCCCCCCGGCCGGGTACCCGCGCTCGCGCCCCGCgccccgcccccccaat                                                                                             >  4:367592/1‑129 (MQ=255)
                                        aaCACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCCCCGCCCCCCCCCGCCCGGTCCCCGCGCTCGCGCCCCGCGCCCCGCACCCCCAAGGCa                                                                                          >  5:206451/1‑135 (MQ=255)
                                        aaCACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCa                                                                                          >  4:259350/1‑135 (MQ=255)
                                           acTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCCTGGGCGGGTCACCCCCACCCCAGGCCCGGTCCCCGCGCTCCCGCCCCGCGCCCCGCCCCCCAATGGCaaac                                                                                       >  2:49829/1‑132 (MQ=255)
                                             tGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGTTCACCGCCAACACAGGCCGGGTCCCAGCGATAGCGCACGGCGCACCCACCCACAATGGCATACGc                                                                                     >  6:29657/1‑135 (MQ=255)
                                             tGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCCCCGCCCCCCCCCGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATcccc                                                                                     >  3:173078/1‑133 (MQ=255)
                                             tGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGc                                                                                     >  4:86054/1‑135 (MQ=255)
                                                cAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCGGGTCCCCGCCCCCCCCGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCCGCACCCCCATGGCCTCCGCCCgg                                                                                 >  7:227702/1‑136 (MQ=255)
                                                  ggCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCACGCCCGGTCCCCCCGCTCGCCCCCGGCGCCCAGCCCCACCATTGCATACGCCCCGc                                                                                >  1:163115/1‑135 (MQ=255)
                                                  ggCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCt                                                                               >  8:282156/1‑136 (MQ=255)
                                                  ggCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCt                                                                               <  7:282156/136‑1 (MQ=255)
                                                  ggCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGGTCCAGTGCCGTGGCCGTGGGCGGGTCACCCCCCCCCCACGCCGGGTCCCCGCCCGCGCGCCCCGCGCCCCGCCCCCCCAAGGCCACCGCCCGGc                                                                                >  1:372213/1‑135 (MQ=255)
                                                          gCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGGCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAggc                                                                       <  8:269968/136‑1 (MQ=255)
                                                                gCCCGCCGCGCGGCAAAGCGCATCCGTTCCCGTGCCGTGGCCGTGGGCGGGTCACCGCCACCCCCGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAggccggc                                                                   <  2:446902/134‑1 (MQ=255)
                                                                gCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCCGCACCACAATGGCATCCGCCCGGCGCGCCAGGCCGGccc                                                                 >  8:152343/1‑136 (MQ=255)
                                                                 ccAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGccc                                                                 <  1:377907/135‑1 (MQ=255)
                                                                  cAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTGCCCGCCCCCCCCCGCCCGGGCCCCCCCCTCGCGCCCGGGGCCCCGCCCCCCCATGGGCTGCGCCCCGCTCGCCCCgcccgcccc                                                                >  4:382151/1‑134 (MQ=255)
                                                                  cAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCCCCGCCCCCCCCCGCCGGGGCCCCGCGCGCGCGCCCGGCGCCCCGCCCCCCACTGGCCTCCGCCCCGCGCGCCCGgcccgccccg                                                               >  7:61624/1‑134 (MQ=255)
                                                                  cAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCGGGTCACCGCCCCCCCAGGGCGGGTCCCCGCGCTCGCGCCCCGCGGCCAGCCCCCCCATGGCATCCGCCCGGCTCCCCAGGCCGGCCCa                                                                >  3:109063/1‑135 (MQ=255)
                                                                  cAGCCGCGAGGCAATGAGCATCCGTTCCACTGACGTGGACGTGGGCTGGTCACCGCCACCCAAGGCCCTGTCCCCGCCCACCCCCCCGGCCACCAGCACCAAAATGCCAACCCCCCCGCAACCCACGCCGGCCCa                                                                >  4:378006/1‑135 (MQ=255)
                                                                      cgcgcgGCAAAGCGCAGCGGTTCCCGGGCGGTGGCGGTGGGCTGGTCACCGCCACCCCCGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCg                                                            <  2:411874/135‑1 (MQ=255)
                                                                        cgcgGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCCCCCCCGGCCGGGTCCCCGCGCGCGCGCCCGGCGCCCCGCCCCACAATGGCATCCGCCCGGCTCGCCCGGCCGGCCCCGTCGCTc                                                         >  8:145967/1‑136 (MQ=255)
                                                                           ggCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGGGGGTCCCCCCCCCCCCCGGCCGGGGCCCCCCGCGCGCGCCCCGCGCCCAGCCCCACCAAGGCCTCAGCCCCCCGCGCCCCGGCGGCCCCGTCGcgccc                                                       >  2:308901/1‑131 (MQ=255)
                                                                             cAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGc                                                    <  7:145967/136‑1 (MQ=255)
                                                                              aaaGCGCAGCCGTGCCAGGGCCGGGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCGCACGGc                                                    <  1:308901/135‑1 (MQ=255)
                                                                                         gTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCggcggc                                          <  3:367592/134‑1 (MQ=255)
                                                                                         gTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCggcggc                                          <  2:57041/134‑1 (MQ=255)
                                                                                                    gtggccgtggGCTGGTCACCGCCACCCCAGGCCGGGGCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGc                              <  1:272809/135‑1 (MQ=255)
                                                                                                     tggccgtggGCGGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGca                             <  6:250245/135‑1 (MQ=255)
                                                                                                        cggtggGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCacac                          <  2:70717/133‑1 (MQ=255)
                                                                                                                cTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCCGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCGCCCCGCCCCCGGCGGCGCCCGTTTCACCACCCcggcggg                   >  5:409050/1‑134 (MQ=255)
                                                                                                                    tCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCACACGGGCGGTTTc                <  3:137363/134‑1 (MQ=255)
                                                                                                                      ccctccACCCCAGGCAGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGTCATCCGCCCGGCGCGCCAGGCCGGACAAGTCTCTCCCGGCAAGCGGCTGATCACTTTCACGCCCACGGGAGGGTTACTg             >  6:141290/5‑135 (MQ=255)
                                                                                                                      aCCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCCCTTTCCAGCACACGGGCGGTTTCCTg             >  4:247827/1‑135 (MQ=255)
                                                                                                                      aCCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCACACGGGCGGTTTCCTGc            >  8:71909/1‑136 (MQ=255)
                                                                                                                        cGCCACCCCAGGCCGGGTCCCCGCGCGCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTCCCAGCACACGGGCGGTTTCCTGcac          <  7:205146/136‑1 (MQ=255)
                                                                                                                                cAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCCCACGGGCGGTTTCCTGCCCCACCCGc   >  1:79071/1‑135 (MQ=255)
                                                                                                                                 aGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCACACGGGCGGTTTCCTGCACCACCAGCg  <  5:208208/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CTCGCTGGGCGGCGTGCCGTAATGCAACCCACCACTGGTCAACACTGGCAGGCCCGAAGCCTTGGCCAGCCGCGCGGCAAAGCGCATCCGTTCCAGTGCCGTGGCCGTGGGCTGGTCACCGCCACCCCAGGCCGGGTCCCCGCGCTCGCGCCCGGCGCCCAGCACCACAATGGCATCCGCCCGGCTCGCCAGGCCGGCCCAGTCGCTCACGGCCAGCGGCGGCTCACTTTCCAGCACACGGGCGGTTTCCTGCACCACCAGCG  >  KT2440/5449002‑5449264

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: