Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,880,375 T→G 33.5% E92A (GAG→GCG)  CJ019_04574 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,880,3750TG33.5% 34.5 / 4.0 47E92A (GAG→GCG) CJ019_04574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base G (13/3);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.78e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCG  >  ACB122/4880246‑4880496
                                                                                                                                  |                                                                                                                         
tcctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCTCCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCAccc                                                                                                                      <  1:371694/135‑1 (MQ=255)
tcctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCAccc                                                                                                                      <  3:41308/135‑1 (MQ=255)
tcctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCAccc                                                                                                                      <  3:279682/135‑1 (MQ=255)
 cctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCCGCCGCGGAcccc                                                                                                                     >  3:407120/1‑134 (MQ=255)
      aTCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCCGCGCCGGCACCCCCGGCCGCGCCCCCAAcccc                                                                                                                >  6:358885/1‑134 (MQ=255)
         gCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGATTCAGCACCCCACCCGGCGCCGGCACCACCCGCCCCGCCCCCCTccccgt                                                                                                              >  7:423022/1‑131 (MQ=255)
            cTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCACCCCCCGCCGCGCCCCCCGCCCGGTCCAc                                                                                                          >  4:279682/1‑135 (MQ=255)
                  gCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCCGCACCCCCCGCCCCGCACCCCCCCCCTTAAAccccgc                                                                                                     >  6:28457/1‑129 (MQ=255)
                  gCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCCGCGCCGGCCCCCCCGGCCGCGCCCCCCTCCCGGTCAAcgccgg                                                                                                     >  7:217012/1‑131 (MQ=255)
                    acacTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGCTCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCGGCCGCTCCCCCATCCCGTGGCACGCGGCaa                                                                                                   >  7:272163/1‑134 (MQ=255)
                    acacTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCCCCCCCCGCCCCGCCCCCATCCCCTTCAACGCGGCaa                                                                                                   >  4:264606/1‑134 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCCCCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCACCCCCCGCCGCTCCCCCCTCCCGGTGAACGCCGCAAAAAAAaccg                                                                                          >  1:347211/1‑135 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCCGCACCCCCCCCGGCGCCGGCCCCCCCGGCCCCGCCCCCAACCCCTTCAACCCGGCAAAAAAAgcag                                                                                          >  8:202133/1‑133 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCCCCACCGGCCGCGCCCCCCTCCCGGTCCACGCGGCAAAAAAAgccg                                                                                          >  3:280091/1‑135 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCCGCCCCACCCGCCGCTCACCCCACCCGGTAAAACCGGCCAAAAAAgccg                                                                                          >  5:36711/1‑135 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCCCCGGCCCCACCCGCCGCGCCCCCCTCCCGGTCAACGCCGCAAAAAAAcccg                                                                                          >  2:66546/1‑135 (MQ=255)
                            cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGAACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCACCCGCCCCCGCCACCACCGCCCGGCCACCCGACCGGTCAACCCGGCaaaaaaaaacc                                                                                          >  2:260672/1‑131 (MQ=255)
                            cgtcgATATCCGACGTATAGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCCACCCGTTCAACGCGGCAAAAAAAgccg                                                                                          >  3:187537/1‑135 (MQ=255)
                             gtcgATATCCGACGTATCGCCGCTGAACCCCCCGGCACCCAGCCCCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGTAACCACGTTCCTCTCACCCATCCCGTTCAACGCGGCAAAAAAAgccgg                                                                                         <  6:179731/135‑1 (MQ=255)
                                            aTCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCCCCCGCCGCGCACCCATACCGTTCAACGCGGCAAAAAAAGCCGGCCCCAGCTTCCCAt                                                                           >  4:182931/1‑134 (MQ=255)
                                              cgccgcTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCACCGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATcc                                                                         <  8:69035/134‑1 (MQ=255)
                                                   cTGATCCCCACGGCACCCACCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGTCACCACCGGCCGCGCACCCATCCCGTTCAACGCGGAAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGcc                                                                   <  6:128843/135‑1 (MQ=255)
                                                          ccACGGCACCCCGCCCCTTGCCGCCCTGATCACCACACACCACCCCACCCGGCGCCGGCCCCACCGGCCGCTCACCCATCCCGTTCAACGCGCCAAAAACAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAgc                                                            <  8:423022/135‑1 (MQ=255)
                                                          ccACGGCACCCAGCACCTTGCCGCCCTGAGCACCAATCAGCACCCCCCCCGGCGCCGGCACCCCCGGCCGCTCACCCATCCCGTTCACCGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAgc                                                            <  1:17582/135‑1 (MQ=255)
                                                               gCACCCAGCCCCTCGCCGCCCTGCTCACGAAGCAGCACCCCCCCCGGCGCCGGCACCCCCGCCCGCGCACCCATCCCGTTCAACGCGGCAAAAAACCCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAg                                                       <  5:70793/135‑1 (MQ=255)
                                                                     aGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCCCCCCCGGCCGCGCCCCCCTCCCCGTGAACGCCGCCAAAAAAGCCCGCCGCGGGGGCGGCTCCCGCCCCCCCCAGCGCGGGCCCCTg                                                 >  5:356080/1‑135 (MQ=255)
                                                                       cACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTACAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGcc                                               <  5:89756/135‑1 (MQ=255)
                                                                             gCCGGCCTGACCACGAATCCGCCCCCCCCCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGc                                         <  1:337557/135‑1 (MQ=255)
                                                                               cGTCCTGATCACGAATCCGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                       <  5:358885/135‑1 (MQ=255)
                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCCGCCGCGCACCCATCCCCTTCAACGCCGCAAAAAAACCCCGCCGCGGCAGCGCAGACAGCGCCCGCCACCGCGAGCCCATGCCCCGGGcca                                       >  6:103957/1‑133 (MQ=255)
                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCCCCGGCCGCTCCCCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCCGCGCCAGCAAGCCCCAGCCCTTGCCCAGGGCaa                                       >  3:223586/1‑135 (MQ=255)
                                                                               cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                       <  7:435987/135‑1 (MQ=255)
                                                                               cGTCCTGATCACCAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGAAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                       <  1:66546/135‑1 (MQ=255)
                                                                                         aCGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCa                             >  2:46486/1‑135 (MQ=255)
                                                                                                aGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCCACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCCTCCAGCGCCCGCAACCGCGCGCCCATGCCCAGGGCAATCGCCCCCCCGGCtttt                      >  8:17036/1‑134 (MQ=255)
                                                                                                 gCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTg                      <  4:357964/134‑1 (MQ=255)
                                                                                                      cccccccGGCGCCGGCCCCCCCGGCCGCGCCCCCCCCCCGTACACCGCGCCAAAACACCCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                 <  2:442479/135‑1 (MQ=255)
                                                                                                           ccGGCGCCGGCACCCCCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCa            <  7:197312/135‑1 (MQ=255)
                                                                                                               cgcCGGCACCACCGGCCGCTCACCCAACCCGTTCCACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTc        <  6:36711/135‑1 (MQ=255)
                                                                                                               cgcCGGCACAACCGGCCGCTCACCCATCCCGTACAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTc        <  6:438727/135‑1 (MQ=255)
                                                                                                                gccgccACCCCCGGCCGCTCACCCAACCCTTTCAACGCGCCAAAAAACCCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCa       <  8:166573/135‑1 (MQ=255)
                                                                                                                gcCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGCCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCa       <  5:103957/135‑1 (MQ=255)
                                                                                                                     caccccCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  3:264606/135‑1 (MQ=255)
                                                                                                                     caccacCGGCCGCTCACCCATCCCGTTCACCGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  4:407120/135‑1 (MQ=255)
                                                                                                                     caccacCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  8:272163/135‑1 (MQ=255)
                                                                                                                     caccacCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg  <  8:380093/135‑1 (MQ=255)
                                                                                                                     caccacCGCCCGCTCCCCCCTCCCGTTGACCGCGCCAAAAAACCCCGGCCGCTGCTGCGCATGCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTAGCCAGTGGCCACTTCAGGCCg  <  2:451055/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                         
TCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCACCCATCCCGTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCG  >  ACB122/4880246‑4880496

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: