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

TGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  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)
                                 aCCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCCTGGGGCATATGTGCGCTGGTCGCCCGCAAAAc                                                                                                       >  6:76147/1‑136 (MQ=255)
                                  ccGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCCTTTTTTATTGTGCATATTTTCGCTGTTTCGCCGCAAAAcc                                                                                                      >  6:193082/1‑136 (MQ=255)
                                        cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTAATTGGGCATATTTTCGCTGTTAGCCCGGAAAAcacccaa                                                                                                 >  5:247409/1‑129 (MQ=255)
                                        cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACCGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATTGGGCATATGTGCGCTCTTAGCCCCGAAAACCGACaa                                                                                                 >  6:210032/1‑135 (MQ=255)
                                                         cacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGGCTATTTTCGGGGTTCGCCCGCCAAACCCCCAAACGACCGGGACAAAAAc                                                                                >  2:300895/1‑135 (MQ=255)
                                                            gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTAAAGGGGGATATTTGCGCTGTTCGACCGCAAAACCCACAAACCCCCCGGATAAAAAACaa                                                                             >  8:273224/1‑135 (MQ=255)
                                                                 gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAATGt                                                                        >  8:289826/1‑135 (MQ=255)
                                                                 gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACAGTc                                                                       >  8:43678/1‑136 (MQ=255)
                                                                   cAAGGGGGCGTAAACGGCGGACGGCCGTAAGTAAACACCAGAAGGACGGGGCCAGAGGCCCCTTTTTGCATGGGGCTTATGTGCGCGGTGCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAt                                                                     <  6:115496/136‑1 (MQ=255)
                                                                     cggTGGCGGAATCGGCTGAAGGCCGGAAGTAACCACCAGAAGGACGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTACGCCCGCAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                                   <  5:254969/135‑1 (MQ=255)
                                                                     aGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAAGTGCGCTGTTAGCCCGCAAAACCGACAAACGCCCCGGAAAAAAACCAAACTGGCATc                                                                    >  1:262274/1‑135 (MQ=255)
                                                                                    cTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATTGGGCAGATTTGGGCGGTTCGCCCGCCAAAACGACAAACGCCCGGGGAAAAAACCAAACGGGCCTACAGACGTGTGTTTCCt                                                    >  3:223886/1‑136 (MQ=255)
                                                                                              taaGTAACCACCAGAAGCAAGGGGCCAGAGGACCCTTTTTTCATGGGGCATATGGGCTCTGTTCGCCCGCAAAACCGACAAACGACCTGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          <  5:49786/135‑1 (MQ=255)
                                                                                               aaGAAAACACAAGAATTAAGTGTCCATAGTCCCCTTTTTTAATGGGTCATATTTGAGCTTTTAGAACGCAAAACAGAAAAACGAACGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          <  5:45048/135‑1 (MQ=255)
                                                                                                       acCAGAAGGACGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                                 <  6:279804/136‑1 (MQ=255)
                                                                                                              ggagtgggCCCCAGGCCCATTTTTAATTGGGGCATATGTGCGCGGTGCCCCCCCAAAACCGACACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                          <  1:9586/131‑1 (MQ=255)
                                                                                                              gcagggggCCCGCGGCCCGTTTTTTCTTGGGGCATTTGGGCGCTGTCCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAAcg                          <  5:76147/132‑1 (MQ=255)
                                                                                                                    ggCCAGAGGCCCGGTTTTTTATTGGGGGGATGTTCGCTGGTCGGCCGCCAAACCGCCCAACGACCGGGAAAAAAACCAAACAGTGATCCAGACTTTTGTTTCCTGCCCCAGCTGTTTTTAGCTACAAACCCGccag                    >  2:309481/1‑134 (MQ=255)
                                                                                                                      cccctcgCCCCTTTTTTTATTGGGCAAATGGGGCGGGTTCCCCCCCAAACCCCCCCCACCGCCCGGCAAAAACACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTgg                   <  6:200649/130‑1 (MQ=255)
                                                                                                                                  tttttCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          <  7:206218/132‑1 (MQ=255)
                                                                                                                                  tttttCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          >  8:206218/1‑132 (MQ=255)
                                                                                                                                  tttttCATGGGCCATAGGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc       <  5:110084/135‑1 (MQ=255)
                                                                                                                                  tttttAAAGGGGCATATGTGCGCTGTTACCCCGCAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCAAGACATGTGTATACTGCCCCGGCTTTTTGAATCTACAAACCCGCCAGGCGCAGACTTTtca      >  3:243603/1‑136 (MQ=255)
                                                                                                                                    tttCATGGGGCTATTGTCCGCTGACCCCCCCAACAACCGCAATACGACCTGGACAAAAACCAAACTGTCATCCAGAATTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcac     <  2:220307/135‑1 (MQ=255)
                                                                                                                                       cATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcactca  <  6:120280/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  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: