Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,544,167 A→G 12.5% E250G (GAA→GGA)  CJ019_05222 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,544,1670AG12.5% 83.8 / 4.8 32E250G (GAA→GGA) CJ019_05222hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (2/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGCATTGTCGATCCTCGGCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACTGTTCTGCGTGGGCCTGGCGGTGCTGGGCA  >  ACB122/5544033‑5544298
                                                                                                                                      |                                                                                                                                   
cgcATTGTCGATCCTCGGCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGaa                                                                                                                                    >  7:150752/1‑136 (MQ=255)
        cGATCCTCGGCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGt                                                                                                                            >  3:151725/1‑136 (MQ=255)
             cTCGGCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCCACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGcgcg                                                                                                                       >  5:106775/1‑136 (MQ=255)
                cgCCTGCACACCTCTCCCCCGGCGAGCTAGCCGCCGGTCGCCCTGGAGTGCAGCAACGCCAACGCGGCCATGCTGTCTGGCGAGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGAGCGCgg                                                                                                                      <  4:167178/133‑1 (MQ=255)
                 gCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCaa                                                                                                                   <  5:55336/136‑1 (MQ=255)
                 gCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCaa                                                                                                                   <  6:57447/136‑1 (MQ=255)
                      cAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGcc                                                                                                               <  1:119057/135‑1 (MQ=255)
                      cAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGc                                                                                                                <  1:115714/134‑1 (MQ=255)
                          ccTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCa                                                                                                          <  4:31772/136‑1 (MQ=255)
                                   cGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCAt                                                                                                 <  6:26690/136‑1 (MQ=255)
                                       gATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGc                                                                                             >  1:166479/1‑136 (MQ=255)
                                        aTCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGc                                                                                             >  3:40318/1‑135 (MQ=255)
                                                       gAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCgcggc                                                                             <  7:177095/136‑1 (MQ=255)
                                                       gAACTGGAGTTCAGAAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGATGTGGCTGGAACTGCTGGTGCGCGGCAAAGAACGCCAACTACGCATCGGCCAGGGCGTAGCgcggc                                                                             >  8:177095/1‑136 (MQ=255)
                                                        aaCTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGGACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCgcggcg                                                                            >  6:39338/1‑136 (MQ=255)
                                                        aaCTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGGACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCgcggcg                                                                            >  6:39337/1‑136 (MQ=255)
                                                             gAGTGCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTgc                                                                         <  3:64506/134‑1 (MQ=255)
                                                                 tCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTAc                                                                                                     <  5:121870/102‑1 (MQ=255)
                                                                 tCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTAc                                                                                                     >  6:121870/1‑102 (MQ=255)
                                                                     caacgccaacgcGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCCACCGGTg                                                               >  8:37015/1‑136 (MQ=255)
                                                                             acgcGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTgcg                                                       <  7:14689/136‑1 (MQ=255)
                                                                                      tGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCCACCCGTGCGGCTGCGCGGCTgggg                                              >  7:118816/1‑135 (MQ=255)
                                                                                                cGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGGACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGgcg                                    <  5:39338/136‑1 (MQ=255)
                                                                                                cGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGGACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGgcg                                    <  5:39337/136‑1 (MQ=255)
                                                                                                cGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGgcg                                    >  4:87805/1‑136 (MQ=255)
                                                                                                       ctgGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACt                               <  1:4950/134‑1 (MQ=255)
                                                                                                                   cTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACTGTTCTGCGTGGGcc                 <  2:166479/136‑1 (MQ=255)
                                                                                                                    tGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACt                               <  5:105020/121‑1 (MQ=255)
                                                                                                                    tGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGGGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACt                               >  6:105020/1‑121 (MQ=255)
                                                                                                                       cAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTTCTGGCGCAACAGTTCTGCGTGGGCCTGGc             >  3:100058/1‑136 (MQ=255)
                                                                                                                        aaTGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCCCCCCGTGCGGCGGCGCGGGGGGGCGGTGGGGGCGCCACGGGTGTGCGTGGGCCGGGCg            >  7:16434/1‑136 (MQ=255)
                                                                                                                        aaTGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACTGTTCTGCGGGGGCCTGGCg            <  4:40318/136‑1 (MQ=255)
                                                                                                                                  cTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCCACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGGGCAACTGTTGTGCGTGGGCCGGGCGGGGCTGGGcc  >  8:11052/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CGCATTGTCGATCCTCGGCCTGCAAACCTTTACCACGGCGATCTACCAGCAGTTCGAACTGGAGTTCAGCAACGCCAACGCGGCCATGCTGTCTGCCGTGCTGCTGGCGATGTGCCTGGCAATGCTGTGGCTGGAACTGCGGGTGCGCGGCAAAGCCCGCCATGTACGCATCGGCCAGGGCGTAGCGCGGCGTGCGCAACCGGTGCGGCTGCGCGGCTGGGCGGTGCTGGCGCAACTGTTCTGCGTGGGCCTGGCGGTGCTGGGCA  >  ACB122/5544033‑5544298

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