Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,880,375 T→G 44.3% 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,3750TG44.3% 6.9 / 6.3 37E92A (GAG→GCG) CJ019_04574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base G (12/5);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCGAATCAGGC  >  ACB122/4880245‑4880504
                                                                                                                                  |                                                                                                                                   
ctcctcGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCCGCCCCACCCGCCCCGCA‑cc                                                                                                                                >  1:263375/1‑135 (MQ=255)
          gCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCCGCCGCTCC‑CCCATCCC‑ggt                                                                                                                      >  3:279632/1‑132 (MQ=255)
            aCTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCaa                                                                                                                   <  4:70272/135‑1 (MQ=255)
                           cTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGCATCAGCACCCCACCCGGCGCCGGCGCCCCCGGCCCCGCC‑CCCCACCC‑GGGCAACGCGGCCAAAAAAgc                                                                                                    >  8:355309/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCCCCCCCCCCCGCGCCGGCCCCCCCGGCCGCGCC‑CCCATCCC‑CGTCAACGCGGCCAAAAAAgccg                                                                                                  >  8:17054/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCGGCCCCACCCGCCGCGAA‑CCCCTCCC‑GGTCAACGCCGCATAAAAAgccg                                                                                                  >  2:160427/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCCGCACCCCCCGCCCCGCA‑ACCATCCC‑GGTCAACGCGGCAAAAAAAcccg                                                                                                  >  6:272043/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCCCCGGCCCCCCCCGCCCCGGACCCCACCCC‑GTGCACGCCCGCaaaaaaaccg                                                                                                   >  2:170722/1‑132 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCGGCCGCGCA‑CCCCTCCC‑GGTCAACGCGGCAAAAAAAgccg                                                                                                  >  1:64414/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCCGCGCCCGCCCCCCCGGCCGCGCC‑CCCCTCCCGGTGCAAC‑CCGCAAAAAAAcccc                                                                                                  >  2:115765/1‑134 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCAACCGGCGCCGGCCCCACCGGCCGCGCA‑CACCTCCC‑GGTCAACGCGGCAAAAAAAgccg                                                                                                  >  1:70177/1‑135 (MQ=255)
                             cgtcgATATCCGACGTATCGCCGCTGATCACCACGGCACCCACCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCACCACCGGCCCCGAA‑CCCCACCC‑CTTCCACTCGGCAAAAAAAgccg                                                                                                  >  6:135403/1‑135 (MQ=255)
                                      ccGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCCGCCCCACCCGCCGCACA‑CCCATCAC‑CTTCAACCCGGCAAAAAAAGCCCGCCgctgca                                                                                         >  4:173166/1‑134 (MQ=255)
                                        gACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCCGCGCCGGCCCCACCGGCCGCGCC‑CCCCTCCC‑GGTCAACGCGGCCAAAACAGCCCGCCCCtgcggc                                                                                       >  3:372523/1‑135 (MQ=255)
                                          cGTATCGCCGCTCATCCCCACGCCACCCAGCACCTGGCCGTCCTGATCACGAATCACCACCCCACCCGGCGCCGGCACCCCCGGCCGCTCA‑CCCAGCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTgcgc                                                                                     <  4:209378/135‑1 (MQ=255)
                                           gTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCCCCCGGCGCCGGCCCCACCGGCCGCGCC‑CCCCTCCC‑CGTCCACGCGGCAAAAAAAGCCGGCCGCTGCTgcgcc                                                                                    >  8:243958/1‑134 (MQ=255)
                                                                gCACCCAGCACCTTGCCGCCCTCATCACGAAGCAGCCCCCCCCCCGGCGCCGGCCCCACGGGCCGCACA‑CCCAGCCC‑TTTCCACGCGGAAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAg                                                               <  7:158206/135‑1 (MQ=255)
                                                                     cAGCACCTTGCCGTCCTGAGCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCtt                                                          <  6:173147/135‑1 (MQ=255)
                                                                         aCCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCCCCGGCCGCACC‑CCCATCCC‑GGTCAACGCCGCAAAAAAAGCCGGCCCCTGCTGCGCCTCCACCGCCAGCCCCCGCGCGCCCATGccc                                                      >  3:67591/1‑135 (MQ=255)
                                                                           cTTGCCGTCCTGATCACCAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCAACCC‑TTTCAACGCGGCAAAAAAACCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAg                                                    <  5:272043/135‑1 (MQ=255)
                                                                                cGTCCTGATCACTAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTAAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                               <  3:173166/135‑1 (MQ=255)
                                                                                cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCCCCACCGGCCGCTCC‑CCCATCCC‑GTTCAACGCCGCAAAAAAAGCCGCCCGCTGCTGCGCCTCCCCCGCCACCCAGCCCGCCCCCATGCCCAgcgcca                                               >  8:305567/1‑133 (MQ=255)
                                                                                cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCCCCCGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCAGCCAGCGCCCGCAACCGCGAGCCCATGGCCAGGGcc                                                >  8:187871/1‑133 (MQ=255)
                                                                                cGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCC‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGCGCCCTTGCCCAGGGCaa                                               >  1:36931/1‑135 (MQ=255)
                                                                                cGTCCTGATCACGAAGCAGCACCCCCCCCGGCGCCGGCACCCCCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCaa                                               <  5:149642/135‑1 (MQ=255)
                                                                                                     ccccacccGGCTCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCTCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAg                          <  2:75837/135‑1 (MQ=255)
                                                                                                      cccacccGGCGCCTGCACCACCGGCCGATCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCTCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTTCCAGt                         <  2:70177/135‑1 (MQ=255)
                                                                                                      cccacccGGCGCCGGCCCCACCGGCCGCTCA‑CCCATCCC‑GTTCCACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                         <  7:266556/135‑1 (MQ=255)
                                                                                                      cccacccGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                         <  8:136105/135‑1 (MQ=255)
                                                                                                      cccacccGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGt                         <  2:64414/135‑1 (MQ=255)
                                                                                                           ccGGCGCCGGCACCACCGGCCGCTGC‑CCCATCCC‑CGTCAACCCGGCAAAAAAACCCCGCCCCGGCTGCGCATACCCCGCCCCCCACCCCCCGCCCCTGCCCCCGGCCAAAACACCCCACGCGTTGCCAGGTGCCa                    >  4:66075/1‑135 (MQ=255)
                                                                                                                     caccacCGGCCGCGCA‑CCCATCCC‑GTTCAACGCGCCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg          <  6:292987/135‑1 (MQ=255)
                                                                                                                     caccacCGCCCGCGCA‑CCCAACCC‑GTGCAACGCGGCAAAAAACGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCg          <  1:213557/135‑1 (MQ=255)
                                                                                                                          ccGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCGAATCa     <  3:66075/135‑1 (MQ=255)
                                                                                                                            ggCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCCGCAAGCGCGAGCCCTTTCCCCGGGCAATCGCACCCCAGGCTTTGCCACTGGCCCCTTCACGCCGAAACAgg   >  4:176965/1‑135 (MQ=255)
                                                                                                                             gCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCCGGCGTTGCCAGTGGCCACTTCAGGCCGAATCAgg   >  8:221430/1‑134 (MQ=255)
                                                                                                                             gCCGCGCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACTCCAGGCTTTGCCAGTGGCCACTTCAGGCCGAATCAGGc  <  4:227061/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
CTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCCCACGGCACCCAGCACCTTGCCGTCCTGATCACGAATCAGCACCCCACCCGGCGCCGGCACCACCGGCCGCTCA‑CCCATCCC‑GTTCAACGCGGCAAAAAAAGCCGGCCGCTGCTGCGCATCCAGCGCCAGCAAGCGCGAGCCCTTGCCCAGGGCAATCGCACCCCAGGCTTTGCCAGTGGCCACTTCAGGCCGAATCAGGC  >  ACB122/4880245‑4880504

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