Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 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
*ACB1223,440,3430AC93.3% 24.7 / ‑5.8 22S18R (AGC→CGC) CJ019_03266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (17/2);  minor base G (0/2);  total (18/4)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440211‑3440473
                                                                                                                                    |                                                                                                                                  
gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCcgcg                                                                                                                                 >  7:348610/1‑136 (MQ=255)
            tATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGcccc                                                                                                                      >  2:152239/1‑135 (MQ=255)
                          gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCAGGgt                                                                                                        >  5:353787/1‑135 (MQ=255)
                          gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCAGGgt                                                                                                        >  3:268435/1‑135 (MQ=255)
                          gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCAGccgcggg                                                                                                        >  5:347843/1‑130 (MQ=255)
                                 tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGGGTgggcgg                                                                                                >  8:140689/1‑132 (MQ=255)
                                        gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACg                                                                                          >  3:303990/1‑135 (MQ=255)
                                        gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCGGGCCACCg                                                                                          >  3:207450/1‑135 (MQ=255)
                                                ggCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCCGGCACCCGGGGTGTGGCAGGCCACCCCCTATAAg                                                                                  >  4:55189/1‑135 (MQ=255)
                                                        aCAAGCGAGCGCTCGAGGTGCTGCAGTGTGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCCCCGCGGGGTGGCCGGCCCACCCCGATAAAGCGCAGcg                                                                          >  3:221637/1‑134 (MQ=255)
                                                               agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGCCCCCCCCGAGAAGGCGCCGCTCGcggc                                                                    >  3:38397/1‑131 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCCCCTGGGTGTGGCAGGCCCACCCCGAAAAGGCGCCGCTCGCGGGTCac                                                                >  2:84182/1‑134 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCCGGCCCACCCCTATAAGGCGCCGCGCGCGGATCAg                                                                >  6:115535/1‑135 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCGGGGTGGGGCCGGCCCACCCCGATAAGGCGCAGCTCGCTGATCAgg                                                               >  8:262967/1‑136 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGGAGGCCCACGCCTAGAAGGCGCAGCGCGCTGATCAg                                                                >  3:172593/1‑135 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGCCGGCCCCCGGGCCGCCGGGGGGTGGGTGGCCCACGCCCAGAAGGCGCCGCACGCTGAtgat                                                                >  4:202769/1‑132 (MQ=255)
                                                                  gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCGGCCCGCCGGGGGGGGGCAGGCCCACGCCCATCAGGCGCAGCTCGCTGAGCAg                                                                >  5:265843/1‑135 (MQ=255)
                                                                    tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCGGCTCAGGc                                                              >  3:199180/1‑135 (MQ=255)
                                                                    tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGGGGCGCCCGGGCCGCCGGGGGGTGGGCGGGCCCCGCCGATAAGGCGCAGCTCGCGGCTCAGGc                                                              >  1:310629/1‑135 (MQ=255)
                                                                        ggggctgcCGTGCGCGCTTGGGCCACCCGCCTGGGCGCCGTGGCGCGGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATg                                                         <  7:262967/133‑1 (MQ=255)
                                                                                                           ggggggggCGGTGCGGCGCGCCGCCGGCGCGGGCGCCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:161523/129‑1 (MQ=255)
                                                                                                           ggcgggggCGCGGCGGCGGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:31527/129‑1 (MQ=255)
                                                                                                                          gcgggCCGCGGGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg       <  7:140689/133‑1 (MQ=255)
                                                                                                                               ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  7:14844/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  ACB122/3440211‑3440473

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: