Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,096,422 C→A 29.0% intergenic (+35/‑155) 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
*ACB1435,096,4220CA29.0% 21.4 / 3.6 33intergenic (+35/‑155)CJ019_04792/CJ019_04793Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/12);  new base A (5/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  ACB143/5096287‑5096556
                                                                                                                                       |                                                                                                                                       
tGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGACGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTc                                                                                                                                         <  3:249373/136‑1 (MQ=255)
    cTGTTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTAATCATGTGTGCAGGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                                      <  4:235263/135‑1 (MQ=255)
     tggtggAGCCACTCGATCAAAGGCCTGCACCGCAGCGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTGATggg                                                                                                                                    >  1:124062/1‑136 (MQ=255)
     tggtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACTGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATtgg                                                                                                                                    >  1:298006/1‑136 (MQ=255)
      tgttgAGCCCCTCGAACAAATGCCAGCACCGCAACGACGTCCTCCTGTGTGCACGCACTGACAATGTTGCGGTAACGGCTGAAGGCCGGAATTAAACACCATAATGAATGGGCCAGAGGCCCCTTTTTTCATgggg                                                                                                                                   <  4:208705/132‑1 (MQ=255)
      ggtggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGATTTTTAATgggg                                                                                                                                   >  6:51797/1‑136 (MQ=255)
         ggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTAATGGGGCa                                                                                                                                 >  2:9586/1‑135 (MQ=255)
           atccACTCGCTCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGTAACGGCTGCAGGCCGGAAGTAAACACCAGAAGGTAGGGGCCAGAGGCCCCTTTTTTCATGGGGCatat                                                                                                                              <  3:99421/134‑1 (MQ=255)
                           gCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                              <  4:307242/136‑1 (MQ=255)
                            ccTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTAGCCCGc                                                                                                             >  1:267696/1‑136 (MQ=255)
                                 aCCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCCTGGGGCATATGTGCGCTGGTCGCCCGCAAA‑Ac                                                                                                       >  6:76147/1‑136 (MQ=255)
                                  ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCCTTTTTTATTGTGCATATTTTCGCTGTTTCGCCGCAAA‑Acc                                                                                                      >  6:193082/1‑136 (MQ=255)
                                        cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATTGGGCATATTTTCGCTGTTAGCCCGGAAA‑Acacccaa                                                                                                 >  5:247409/1‑129 (MQ=255)
                                        cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACCGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATTGGGCATATGTGCGCTCTTAGCCCCGAAA‑ACCGACaa                                                                                                 >  6:210032/1‑135 (MQ=255)
                                                         cacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGGCTATTTTCGGGGTTCGCCCGCCAA‑ACCCCCAAACGACCGGGACAAAAAc                                                                                >  2:300895/1‑135 (MQ=255)
                                                            gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTAAAGGGGGATATTTGCGCTGTTCGACCGCAAA‑ACCCACAAACCCCCCGGATAAAAAACaa                                                                             >  8:273224/1‑135 (MQ=255)
                                                                 gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAAATGt                                                                        >  8:289826/1‑135 (MQ=255)
                                                                 gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACAGTc                                                                       >  8:43678/1‑136 (MQ=255)
                                                                   cAAGGGGGCGTAAACGGCGGACGGCCGTAAGTAAACACCAGAAGGACGGGGCCAGAGGCCCCTTTTTGCATGGGGCTTATGTGCGCGGTGCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                                     <  6:115496/136‑1 (MQ=255)
                                                                     cggTGGCGGAATCGGCTGAAGGCCGGAAGTAACCACCAGAAGGACGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTACGCCCGCAAA‑ACCCACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                                   <  5:254969/135‑1 (MQ=255)
                                                                     aGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAAGTGCGCTGTTAGCCCGCAAA‑ACCGACAAACGCCCCGGAAAAAAACCAAACTGGCATc                                                                    >  1:262274/1‑135 (MQ=255)
                                                                                    cTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATTGGGCAGATTTGGGCGGTTCGCCCGCCAA‑AACGACAAACGCCCGGGGAAAAAACCAAACGGGCCTACAGACGTGTGTTTCCt                                                    >  3:223886/1‑136 (MQ=255)
                                                                                              taaGTAACCACCAGAAGCAAGGGGCCAGAGGACCCTTTTTTCATGGGGCATATGGGCTCTGTTCGCCCGCAAA‑ACCGACAAACGACCTGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          <  5:49786/135‑1 (MQ=255)
                                                                                               aaGAAAACACAAGAATTAAGTGTCCATAGTCCCCTTTTTTAATGGGTCATATTTGAGCTTTTAGAACGCAAA‑ACAGAAAAACGAACGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          <  5:45048/135‑1 (MQ=255)
                                                                                                       acCAGAAGGACGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                                 <  6:279804/136‑1 (MQ=255)
                                                                                                              ggagtgggCCCCAGGCCCATTTTTAATTGGGGCATATGTGCGCGGTGCCCCCCCAAA‑ACCGACACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                          <  1:9586/131‑1 (MQ=255)
                                                                                                              gcagggggCCCGCGGCCCGTTTTTTCTTGGGGCATTTGGGCGCTGTCCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                          <  5:76147/132‑1 (MQ=255)
                                                                                                                    ggCCAGAGGCCCGGTTTTTTATTGGGGGGATGTTCGCTGGTCGGCCGCCAA‑ACCGCCCAACGACCGGGAAAAAAACCAAACAGTGATCCAGACTTTTGTTTCCTGCCCCAGCTGTTTTTAGCTACAAACCCGccag                    >  2:309481/1‑134 (MQ=255)
                                                                                                                      cccctcgCCCCTTTTTTTATTGGGCAAATGGGGCGGGTTCCCCCCCAAACCCCCCCCACCGCCCGGCAAAAACACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTgg                   <  6:200649/130‑1 (MQ=255)
                                                                                                                                  tttttCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          <  7:206218/132‑1 (MQ=255)
                                                                                                                                  tttttCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          >  8:206218/1‑132 (MQ=255)
                                                                                                                                  tttttCATGGGCCATAGGTGCGCTGTTCGCCCGAAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc       <  5:110084/135‑1 (MQ=255)
                                                                                                                                  tttttAAAGGGGCATATGTGCGCTGTTACCCCGCAAA‑ACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCAAGACATGTGTATACTGCCCCGGCTTTTTGAATCTACAAACCCGCCAGGCGCAGACTTTtca      >  3:243603/1‑136 (MQ=255)
                                                                                                                                    tttCATGGGGCTATTGTCCGCTGACCCCCCCAACA‑ACCGCAATACGACCTGGACAAAAACCAAACTGTCATCCAGAATTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcac     <  2:220307/135‑1 (MQ=255)
                                                                                                                                       cATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcactca  <  6:120280/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAA‑ACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  ACB143/5096287‑5096556

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: ‑