Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,507,078 A→C 44.5% E31A (GAG→GCG)  CJ019_05179 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,507,0780AC44.5% 6.3 / 2.5 22E31A (GAG→GCG) CJ019_05179hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (3/7);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCT  >  ACB143/5506943‑5507204
                                                                                                                                       |                                                                                                                              
gCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGa                                                                                                                                >  6:7643/1‑136 (MQ=255)
  tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCgcgc                                                                                                                              >  6:211729/1‑136 (MQ=255)
  tCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGc                                                                                                                              >  3:70558/1‑136 (MQ=255)
       cgcgAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGCGCCggcg                                                                                                                         >  7:9388/1‑135 (MQ=255)
        gctaataaTCCAACACCTATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAACCCCCCTGGCAACCCACGCCAGCGAGGGCTAATTCTTTCTCACCGCCAACCCCGCCGAGCAggccg                                                                                                                        <  7:46444/130‑1 (MQ=255)
        gccaataaTCCAACACCGATTTCCCGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCAAATCCCCCTGGCAACCCACGCCAGCGAGGGCTAATTCTTTCTCACCGCCAAGCCCGCCGAGCAggccg                                                                                                                        <  2:263920/130‑1 (MQ=255)
        gccaagaaTCCAACACCGATTTCCAGGATTCCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCTAATTCTTTCTCACCGCCAAGCCCGCCGAGCAggccg                                                                                                                        <  5:182931/133‑1 (MQ=255)
             taaTCCAACACCTATTTCCAGGATTGCCCGAATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAACCCCGCCGAGCAGGCCGGCCTg                                                                                                                   <  7:42036/135‑1 (MQ=255)
                             tCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGCGCACGCCCGCCTCCTCCAATACGGGgcag                                                                                                   >  6:207988/1‑136 (MQ=255)
                             tCCAGGATTGCCCGCATGAACAGACTCCTGCTCCCCCTCCTGCTGCAAATCCCCCTCCCAACCCCCGCCAGCGCGGGCTAATTCTTTCTCACCGCCACGCCCGCCGCGCCGGCCGGCCTGCTCGAAGGCTGGgca                                                                                                    <  8:9388/135‑1 (MQ=255)
                             ccccggATCGCCCGCATGCTCAGACTCCTGCCCACCCTCCTGCTGAAAACCCCCCTGCCAACCCACCCCCGCGCGGGCTAATTCTTTCCCCCCGCCACCCCCGCCGCGCCGGCCGGCCGGCTCGAAGGCTGGgca                                                                                                    <  4:221371/131‑1 (MQ=255)
                                          gCATGCTCAGACTCCTGCTCACCCTCCTGCTGATAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCCGCCCGCCCTCGGCAAGGCTGGCCAGCCCAACCCCATGc                                                                                      >  2:247684/1‑136 (MQ=255)
                                                                 ctcctGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCCGCCAGGTAGAACGGAGGGCCGCACAAACCCCAGCCGCCCGCCAGCCCCTGCTggca                                                                >  1:222963/1‑133 (MQ=255)
                                                                             aTCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGCGGCGGCCGGCCGGCGCGGCGGGGGGGCGGCACCACCCCCGGCCCCGCGCCGGCCGCGGCGGGGCAACCTGCGCCaac                                                   >  3:136054/1‑135 (MQ=255)
                                                                             aTCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCTCGACGGCTGGGCCGCCCACCCCCCTGCCGCCCGCCGGCCGCGGCGGGGAAACCTGCGcccac                                                   >  4:175776/1‑133 (MQ=255)
                                                                             aTCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCCGGCCCGCCCGCGCGCAGGCTGGGGCGCACCCCCCCAAGCCCCGCGCCGGCCGCGGGTGGGAAACCAGCGCCaac                                                   >  1:217044/1‑135 (MQ=255)
                                                                                         cccccCGCCCGCGAGGGCTAATTCTTTCCCCCCCCCCCCCCCCCCCCGCCGGCCGGCCGGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGcca                                       <  1:48791/135‑1 (MQ=255)
                                                                                             ccGCCAGCGCGGGCTATTTTTTTCCCCCCCCCACGCCCGCCGCGCCGGCGGGCCGGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCaa                                   <  3:64981/135‑1 (MQ=255)
                                                                                                     gAGGGCTAATTCTTCCCCCCCCCCCCCCCCGCCGCGCCGGCCGGCCGGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAgc                           <  5:128780/136‑1 (MQ=255)
                                                                                                              ttctttctCCGCCGCAGGCCCGGTGCGCCGCCCGGCGTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTg                   <  3:59375/135‑1 (MQ=255)
                                                                                                                ctttctCACCGCCAAGCCCGCCGAGGCGGCCGGGCGGCGGGAACGGTGGGGCGCACACCCCCCTGCCCCCCGCCCGCCCCGGGGGGAAAACCAGCGCCAACGCCGCCACGCCCCCAAAGAACAGCCCAACGTgggc                >  8:206381/1‑136 (MQ=255)
                                                                                                                ctttctCACCGCCAAGCCCGCCGAGCCGGCCGGCCGGCTCGCAGGCTGGGGCGCCCCACCCCCTGCCGCGCCCCTGCCGCGGGGGGCAAACCAGCGCCCACGCCGCCTCGCCCCCCAAGAAAAGCGCCACGTgcgc                >  4:234758/1‑136 (MQ=255)
                                                                                                                              cgcccgcccACCCGGCCGCCCGGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCt  <  6:3876/135‑1 (MQ=255)
                                                                                                                              cccccGCCGCGCGGGCCGGCGTGCGCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCt  <  7:62374/134‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
GCTCACCCGCGAAGAATCCAACACCGATTTCCAGGATTGCCCGCATGCTCAGACTCCTGCTCACCCTCCTGCTGCTAATCCCCCTGGCAACCCACGCCAGCGAGGGCGAATTCTTTCTCACCGCCAAGCCCGCCGAGCAGGCCGGCCTGCTCGAAGGCTGGGCAGCACAACCCGATGCCGCGCGCCTGCCGCTGCTGGAAAACCTGCGCCAAGGCCGCATCGCCACCAATGACAAGCGCAAGGTGCGCCTGAACAACCGCCT  >  ACB143/5506943‑5507204

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