Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,094,808 A→C 51.0% L186R (CTG→CGG)  CJ019_01995 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,094,8080AC51.0% ‑6.0 / 2.8 18L186R (CTG→CGG) CJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base C (5/4);  total (10/8)
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

CCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAA  >  KT2440/2094686‑2094938
                                                                                                                          |                                                                                                                                   
ccATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCCCGCGCGCCGG‑Gc                                                                                                                        >  1:170050/1‑135 (MQ=255)
                          gCCCGGGCCAACAGTTTCTCAAGATACGTCATCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTTCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCCCGCCGGCGCCGGG‑GGGCCCGTTTTCAGGGccagg                                                                                              >  4:251312/1‑132 (MQ=255)
                                                            tctattCATCCACCTGCGGCGCCGGCTGCTTGCTCAACTCCTGCACGCCGCGCGCGCCGGCCAGGCGCGCCGG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                           <  8:180593/132‑1 (MQ=255)
                                                              tctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCAG‑GCCGGGCTGGCCGTTTTGCAGGGCCAGTTGGCCGAGGTTGGTGAAGAACTCTATATgccgcagc                                                          >  4:57226/1‑135 (MQ=255)
                                                               ctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCCCGCCAG‑GCCGGGCTGGCCGTTTTTCAGGGCCGGTTTGCCGAGGTTGTTTAATAACTCTATATgcagcagca                                                         >  4:103628/1‑135 (MQ=255)
                                                                             ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTG‑GCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGa                                          >  8:143458/1‑136 (MQ=255)
                                                                             ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCGCGCCGG‑GCCGGGGGGGGCTGTTTTTAGGGCCTGTTGGGCGCGGTTGTTGAAGAACTCTAAATGGAGCAGCACCCGGATATGgggg                                           >  3:35134/1‑135 (MQ=255)
                                                                               cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTG‑GCCGGGCTGGCCGTTTTGCCGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCGATATGCAGCCGCACCAGGATATGGCGGagt                                        >  8:134797/1‑134 (MQ=255)
                                                                                 ccGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCGGGCGCGCCTG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                       <  3:8377/135‑1 (MQ=255)
                                                                                      tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGG‑GCCGGGCTGGGGGTTTTGCAGGGGCGGTTGGGCGCGGGTGTTGTAGAACTCTATATGCAGCAGCCCCAGGATATGGCGGCTTTACTTGa                                 >  7:129518/1‑136 (MQ=255)
                                                                                      tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCTG‑GCCCGGCGGGCCGTTTTTCAGGGCCGGGTGGCCGCGGGTGTTGAAGAACTCGATATGCAGCAGCACCCGGTAATGGCGGATTTACTTgc                                 >  8:224612/1‑135 (MQ=255)
                                                                                      tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGG‑GCCGGGGGGGGCGGTTTTGAGGGCCGGTTGGGCGCGGTTGTTGTAGCACTCGTTATGCAGCAGCCCCAGGATATGGCGGCTTTTCttt                                  >  2:6/1‑134 (MQ=255)
                                                                                        cttcctcAAATCCTGCCCGCCCGGCTCGCGGCCCCCGCGCGCCGG‑CCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc                                <  3:103628/129‑1 (MQ=255)
                                                                                                         cgccgcGCGCGACGGCCAGGCGCGCCGG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa               <  4:37199/135‑1 (MQ=255)
                                                                                                           ccgcGGTCGCCGCCCAGGCGCGCGGG‑CCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa             <  5:211052/135‑1 (MQ=255)
                                                                                                                   gccggCCCGGCGCGCCTG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag    <  8:161878/133‑1 (MQ=255)
                                                                                                                    ccggCCCGGCGCGCCGG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag    <  3:140142/133‑1 (MQ=255)
                                                                                                                     cggCCCGGGGCGCCGG‑CCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaagaa  <  7:75528/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                                   
CCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTG‑GCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAA  >  KT2440/2094686‑2094938

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