Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 150,804 A→C 48.2% W7G (TGG→GGG)  CJ019_00156 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122150,8040AC48.2% 1.6 / 7.8 14W7G (TGG→GGG) CJ019_00156hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base C (4/3);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

TCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTC  >  ACB122/150683‑150929
                                                                                                                         |                                                                                                                             
tcgtcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCGCAAGCTGCAGGGCGGCCCAGGCGGGCGTGCCGCGCAGCGGCCGGGCGCCGCCCAGCAGCCCGGCACCCACgcgcg                                                                                                                 >  3:167494/1‑135 (MQ=255)
                        ttAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCACCTCAAGCTGGACGTGGGCCAACGCGGGCGTGCCCACCAGCGGCCCGGCGCCGCCCACCAGCCCCGCCCCCACCCCCGTGACAATGCGAGTATTGCttaacc                                                                                         >  4:66271/1‑132 (MQ=255)
                         tAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAACCGGCAGGTCGGCCAACGCGGGCGTGCCGCGCCGCGTCACGGCGCCGCCCAGCAGCCCCGCCCCCACCCGCTTCACAGTGCGTGTATGGCGTCAGcc                                                                                        >  6:97150/1‑136 (MQ=255)
                                        cGAGGCCGGACGCCTGCACCGCAAGCAGCAGCTGGGCCAAGGCGGGCGGGCCGCGCAGCGGCAGGGCGCAGCCAAGCACCCCCGCACGCCCGCGGGTAACCCTGGGCGTATGGGTTAACCCAAAACCAGGTCgggg                                                                         >  4:168461/1‑133 (MQ=255)
                                                                  ggccggTGGGCCAAGGGGGGGGGGCCGAGCAGGGTCAGGGCGCCGCCAAGCAGCCCGGCACGAACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc                                               <  2:32456/132‑1 (MQ=255)
                                                                  gcacgctCGCCCAAGGCGGGCGTGCCGAGCCGCGTCAGGGCGCAGCCAACCCCCCCGCAACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc                                               <  7:126180/130‑1 (MQ=255)
                                                                    ctgggcGGCCATGGGGGGCGGCCCGTCCTGCGTCAGGGCCCAGCCAAGCAGCCAGCCACGCAGCCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc                                               <  2:49806/129‑1 (MQ=255)
                                                                    cggggcGGACAAGGCGGGCGTGCCGAGCCGGGTCAGGGCGCTGCCAAGCAGCCAGGACCGAACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACtcgctcgc                                             <  4:9417/131‑1 (MQ=255)
                                                                    cgggTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCTGGGCGCAGCAAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACtcgctcgc                                             <  3:66271/134‑1 (MQ=255)
                                                                      gtgcGGCAATGGCGGGCGGCCGGAGCCGCGGCAGGGCGCAGCAAAGCAGCCTGGCACGAACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCgcgc                                           <  2:1168/133‑1 (MQ=255)
                                                                      gggcGGCCAAGGCGGGCGTGCCGCGCAGCGGCAGGGCGCAGCCAAGCAGCCAGGCACGAACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCgcgc                                           <  6:144920/133‑1 (MQ=255)
                                                                                         gCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCt                        >  2:185622/1‑136 (MQ=255)
                                                                                                    tctgggCGCAGCCAAGCAGCCCGGCACGCACGCGCTTCACAGTGCGAGTAGGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTAcg             <  3:158330/133‑1 (MQ=255)
                                                                                                               ccAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCGCGCGTTGTACGCCGCGGCTGTc  >  1:69349/1‑136 (MQ=255)
                                                                                                                         |                                                                                                                             
TCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTC  >  ACB122/150683‑150929

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: