Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,440,348 T→G 54.2% A19A (GCT→GCG CJ019_03266 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,440,3480TG54.2% ‑3.5 / 6.2 20A19A (GCT→GCGCJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (8/3);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTG  >  KT2440/3440228‑3440478
                                                                                                                        |                                                                                                                                  
aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCgcgcc                                                                                                                     >  7:140512/1‑135 (MQ=255)
aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCCGCGCCCCGGc                                                                                                                      >  3:28060/1‑135 (MQ=255)
                            cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCgat                                                                                          >  8:158657/1‑135 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCACCGCCTATAAGGCGCAGCGCGCGGATCAgg                                                                    >  7:25957/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCGGCCTGGGTGTGGCCGGCCACCGCCTATAAGGCGCCGCTCGCGGATCAg                                                                     >  3:115886/1‑135 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCCGGCAGCCGGGGGGTGGGAGGCCCACCCCAATAAGGCGCCGCTCGCGGAACAg                                                                     >  3:90514/1‑135 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCGGCCCGGGGGTGGCAGGCCCACGCCGATAAAGAGCCGCTCGCTGCACAGGc                                                                   >  1:198952/1‑135 (MQ=255)
                                                        gTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCAGGCGGCGCCCGCGCAGCCGGGGTGTGGGAGGCCCACCCCAATAATCCGCCGCTCGCTGATCAGGCCGAtt                                                              >  8:169172/1‑134 (MQ=255)
                                                               aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCCGGCCGCCGGGGGGTGGGCGGCCACCCCCCATAATGCGCCGCTCGCGGCTCAGGGCGATGACTTCGGc                                                      >  7:261460/1‑136 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGCGCCGCCGGGGGGTGGCGGGGCACCGCCGATAAGGCGCAGCTCGCGGATCACGCCGCGGACTTCGGcgaccac                                                >  5:133802/1‑129 (MQ=255)
                                                                       cgggggCCCCCCGCGGGGGCGGCGGGGCGCGGCGGCGGCCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                               <  1:149296/132‑1 (MQ=255)
                                                                        ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCAGCCCGGGCGGCGCCCCCGCCGCCAGGGGGTTGGGCGCCCCCGCCGATTAGTCGGACCTCGCGGGTCATGCCGCAGACATCGGCCCGCAGTTg                                             >  8:74410/1‑136 (MQ=255)
                                                                                          ggggggggCGCGGCGGCGGGCGCCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                            <  3:177455/129‑1 (MQ=255)
                                                                                                ggCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGGCCGGGGCGCCGGGGGGTGGCAGGCCACCCCCCATAAGGCGCAGCTCGCGGCTCAGGCCGCGGCCGGCGGCACGCCGTTGGCGGAGAAGGCTCgcgcgccccgg                     >  8:256503/1‑131 (MQ=255)
                                                                                                  cgcTGCTGCTGGCCGCCAGCGCTAGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                             >  7:244389/1‑110 (MQ=255)
                                                                                                  cgcTGCTGCTGGCCGCCAGCGCTAGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                             <  8:244389/110‑1 (MQ=255)
                                                                                                          cgggCCGGCTGGGCTGGCGTCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCTCACGGTCTGGAGCgg            <  6:115585/133‑1 (MQ=255)
                                                                                                              ccgccCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa       <  7:63704/136‑1 (MQ=255)
                                                                                                               gcgctgcgcTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa       <  1:23684/130‑1 (MQ=255)
                                                                                                                    gcgcTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGCGGACTTAGGCAAGCAGTTGGCGGAGAAGGCTCGCGCGCCCGGTCTGGAGCGGGAGAAGGCTg  >  2:152188/1‑135 (MQ=255)
                                                                                                                        |                                                                                                                                  
ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTG  >  KT2440/3440228‑3440478

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: