Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,097,064 C→G 40.6% intergenic (+28/‑162) CJ019_04792 → / → CJ019_04793 Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,097,0640CG40.6% 21.9 / 8.8 29intergenic (+28/‑162)CJ019_04792/CJ019_04793Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base G (9/3);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGC  >  ACB122/5096931‑5097189
                                                                                                                                     |                                                                                                                                
cACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGccccg                                                                                                                                 >  7:398585/1‑134 (MQ=255)
    gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCTCAACGACTTCATCATGTGGGCACGCACTGACAAGGTGGCGTAAACGGCTGAAGGCCGGAAGTAAACACCAGACGGAAGGGGCCAGAGGCCCCttttt                                                                                                                             <  4:178841/135‑1 (MQ=255)
    gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGGGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                             <  5:55385/135‑1 (MQ=255)
     tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGgtttt                                                                                                                             >  5:57046/1‑134 (MQ=255)
      gACCTGGTGGATCCACTCTATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCTGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGTGCCAGAGGCCCCTTTTTTc                                                                                                                           <  6:149100/135‑1 (MQ=255)
      gACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCttttt                                                                                                                             >  8:36630/1‑133 (MQ=255)
      gACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                           <  3:179457/135‑1 (MQ=255)
      gACCTGGTGGAGCCACTCGAGCAACGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACGGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                           <  1:155878/135‑1 (MQ=255)
         cTGGTGGAGCCACTCGATCAAACGCCGGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGACGGACGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                        <  6:378545/135‑1 (MQ=255)
          tggtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTGTTTTCATgg                                                                                                                       >  2:54791/1‑135 (MQ=255)
                   cACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAtgt                                                                                                             <  4:462758/135‑1 (MQ=255)
                     cTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCATATGTgc                                                                                                           >  3:152384/1‑135 (MQ=255)
                                gCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCCTATGTGGGCGGTTCGccc                                                                                              >  3:352584/1‑135 (MQ=255)
                                      aCCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTCATGGGGGATATGTGGGCTGGTCGCCCCCaaaa                                                                                        >  8:298287/1‑135 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaa                                                                                 >  8:34550/1‑135 (MQ=255)
                                             cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTAATGGGGGATTTGTTCGCTGGTGGCCCGCCAAACCGCCac                                                                                 >  6:42251/1‑134 (MQ=255)
                                                cTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTTCTGGGGCATATTTTCGGTGTTCGCCCGCAAAACCGCCACACg                                                                              >  8:189506/1‑135 (MQ=255)
                                                     tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGGTTTTTATTGGGGGTCTGTGAGCGGTTGGCGCGCCAAACCGCCCAaccacccc                                                                         >  5:293019/1‑133 (MQ=255)
                                                             gcacgcacTGCCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCCACAAACGACCGGGaaaaaaa                                                                 >  7:53871/1‑135 (MQ=255)
                                                             gcacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCATATGTGCGCTGTTGGCCCGCAAAACCGACAAACGACCGGGaaaaaaa                                                                 >  7:61408/1‑135 (MQ=255)
                                                                 gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTTATGGGGCCTATTTTCGCTGGTGGGCCCCCAAAACCCCAAACCACCGGGAAAAAAAcaa                                                              >  2:348787/1‑132 (MQ=255)
                                                                   acTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGGATATTTGGGCTGTTCGCCCGCAAAACCGACAAACGACCGCGAAAAAAACCcaaa                                                           >  1:54194/1‑134 (MQ=255)
                                                                           ggTGGCGGAAAGGTCTGAGGGCCGTAAGTAAACACAACATGGACGGCGCCAGCGCCCCATTTTTACATGGGTCATTTGGGCGTGTTCCCCCCGCAAAACCGCCAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                   <  8:220559/135‑1 (MQ=255)
                                                                                            aaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTGTTTGATGGGGCATATGTGCGCTGTTGGCCCGCCAAACCGACCAACGCCCGGGAAAAAAACCAAACTGTCATCCAGCCTTGTGTATcaggc                                  >  4:7494/1‑131 (MQ=255)
                                                                                              ggCCGTAAGAAAACACCACAAGGATGGGGCCAGCGGCCCCTTTTTACTTGGGGCATATGGGCGCGGTGCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGccc                                <  3:120126/135‑1 (MQ=255)
                                                                                                    aaGTAAACACCAGAAGGAAGGGCCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                          <  7:309248/135‑1 (MQ=255)
                                                                                                                 aaggaaggGGCCAGAGGCCCCGTTTTTCATGGGGGATATTTGCGCTGGTGGCCCGCAAAACCGCCAAACCACCCGGAAAAAAAACAAACAGTCAACCAGACATGTGTATACTGCCCCCCCGTTTTGAATCAACa              >  4:281020/1‑134 (MQ=255)
                                                                                                                     acggggCCAGCGGCCCGTTTTTTCATGGGGCATAGGTGCGCTGAGCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcgc         <  8:128799/133‑1 (MQ=255)
                                                                                                                         ggCCAGAGGCCCGGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACACACGCCCGGGAAAAAAACCAAACTGGGATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGATGCTACCAACGCGccg     >  8:223283/1‑134 (MQ=255)
                                                                                                                            acgagGCCCGTTTTTACTGGGGCAATAGGGGCGGGGTCCGCCCGAAAACCCCCCAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgc  <  2:18955/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                
CACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGC  >  ACB122/5096931‑5097189

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: