Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,440,343 A→C 100% S18R (AGC→CGC)  CJ019_03266 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,440,3430AC80.9% 12.1 / ‑1.1 26S18R (AGC→CGC) CJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base C (17/3);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 3.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  KT2440/3440211‑3440468
                                                                                                                                    |                                                                                                                             
gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                             >  5:128638/1‑135 (MQ=255)
                 aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCACCCCGGCCGGCGCCCGGGc                                                                                                            >  2:121317/1‑135 (MQ=255)
                 aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGc                                                                                                            >  5:33406/1‑135 (MQ=255)
                 aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAACCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCCCCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCCTGCGGCGCCCGGGc                                                                                                            >  6:211104/1‑135 (MQ=255)
                          gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCGGGCGGCGCCCGCGCCGCCTgcgg                                                                                                   >  5:449409/1‑132 (MQ=255)
                                      aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCCGCCGGGGGGTggccggccca                                                                                       >  3:13964/1‑133 (MQ=255)
                                        gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGGGGCCAACg                                                                                     >  5:366223/1‑135 (MQ=255)
                                             cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCACCGccga                                                                                 >  6:219023/1‑132 (MQ=255)
                                                       tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCTGCGGGTGGCCGGCCACCGCCTAAAAGGCGCAGcg                                                                     >  7:1115/1‑135 (MQ=255)
                                                              gagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCACGCCCCCGCCGATAATGCGCCgcgcgcggg                                                               >  3:258480/1‑132 (MQ=255)
                                                               agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCGGGGTGTGGCACGCCAACCCCAATAAAGCGCCGCGCGCGGa                                                               >  6:220371/1‑134 (MQ=255)
                                                               agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCAGGCCCACCCCAATAAGGCGCCGCTCGCTgc                                                               >  4:373420/1‑133 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCACCCCCATAAGGCGCAGCCCGCGGaggag                                                           >  3:113890/1‑131 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGTGGGGGCGGCCACCGCCGAGAAGGCGCAGCGCGCGGGTCAg                                                           >  2:3975/1‑135 (MQ=255)
                                                                    tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGGAGGCCCACGCCGATAAGGCCCCGCGCGCGGAGCCGGc                                                         >  4:319434/1‑135 (MQ=255)
                                                                                     gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGGCCGGGCCGCCGGGGTGGGGCCGGCCCACCCCGATAAGGCGCCGCGCGCGGAGCAGGCCGCGGCCTTCGGcca                                          >  3:144455/1‑131 (MQ=255)
                                                                                       gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCACCCCCGATAATGCGCCGCGCGCGGCTCCGGGCGCTGCCGTCGGGACGCAg                                      >  2:292698/1‑135 (MQ=255)
                                                                                         ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCGATAAGGCGCAGCTCGCTGAGCAGGCCGAGGACTTCGGCCAGCAGtt                                    >  4:168465/1‑135 (MQ=255)
                                                                                                          gcGCCGGGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  7:73747/135‑1 (MQ=255)
                                                                                                           ggCCGGGGCGCGGGGGCGGGCCGCGTGCGCGCGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:211104/134‑1 (MQ=255)
                                                                                                           cgggggggCGCGGCGGCGGGCCGCCGGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  6:128638/129‑1 (MQ=255)
                                                                                                           cgcgggggCGCTGGGGCGGGGCGCCCGCGCGGGCGCCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  5:220371/129‑1 (MQ=255)
                                                                                                           cgcgggggCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  3:224472/129‑1 (MQ=255)
                                                                                                           cGCCGGGGCGGGGCGGCGGCCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCACGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  4:406673/135‑1 (MQ=255)
                                                                                                                tGGCGCTGCTGCTGGCCGCCCGCCCAGGGGGCGGCCCGGCCGCCGGGGGGTGGGAGGCCCCCGCCGATAATGCGGACCGCGCGGCTCAGGGGGCGTACTTAGGGAAGCACTTGGGGGGGGAGGTTCGCGCGCCCg             >  4:116564/1‑135 (MQ=255)
                                                                                                                     ctgctgctgGCCGCGAGCGCTGACGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg        <  2:9498/135‑1 (MQ=255)
                                                                                                                           gcggCCGCGCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  1:121317/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  KT2440/3440211‑3440468

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: