Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,096,964 G→C 36.8% intergenic (+70/‑120) 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
*KT24405,096,9640GC36.8% 28.2 / 8.2 29intergenic (+70/‑120)CJ019_04792/CJ019_04793Putative pterin‑4‑alpha‑carbinolaminedehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/2);  new base C (5/6);  total (20/9)
Fisher's exact test for biased strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGCATGAAATCC  >  KT2440/5096834‑5097096
                                                                                                                                  |                                                                                                                                    
cGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTAATGGGGCATATGTGCGCTGTTGGCCCGCAAAACCGACaa                                                                                                                                  >  1:115393/1‑135 (MQ=255)
        tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGCCCggg                                                                                                                         >  3:108865/1‑136 (MQ=255)
                gcacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGGGTATTTGGGCGGTTAGCCCGCAAAAACCCCAAACGCCCGGGaaaaaaa                                                                                                                  >  5:302204/1‑135 (MQ=255)
                gcacgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCGTTTTTCATGGGGGATATTTGCGGTGTTCGCCCGCCACACCCCCAAACCCCCCGGcaaaaaa                                                                                                                  >  2:123876/1‑135 (MQ=255)
                    gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGTTTTTGATGGGGCATATTTTCGGTGGTAGCCCGCCAAAACCACAAAACACCCGGGAAAAAACCaa                                                                                                              >  3:5420/1‑135 (MQ=255)
                         gACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTg                                                                                                          >  7:94312/1‑134 (MQ=255)
                               ggggCGTAATCGGCTTATGGCGGTATTTAACCACCTGATGGATGGGGCCAGAGGCCCTTTTTTTCATGGGGAATTTGGGCGCTGTGCGCCCGCAAACCCTACAAACGACCGGGAAAAAAACCAAACTGTCATcc                                                                                                    <  7:245436/132‑1 (MQ=255)
                                 ggCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGCGGCCCGTGTTTTGATGGGGCATATGTTCGGTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCAACCAg                                                                                                  >  6:179914/1‑134 (MQ=255)
                                     gAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGGGTTTTGATGGGGGATATTTGGGCTGTTCGCCCGCAAAAACGACAAACGACCGCGAAAAAAACCAAACTGTCATCCaacatag                                                                                             >  2:13558/1‑129 (MQ=255)
                                                       aaGTAAACACCAGAAGGAAGGGGCCAGAGGCCCGTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTTTCCTGCCCGAGCtt                                                                           >  7:56647/1‑135 (MQ=255)
                                                                    aaggaaggGGCCAGAGGCCCGTTTTTTAATGGGGGGGAATTGCGCTGGTGGCCCCCAAAACCCACCAACACCCCCGAAAAAAACCAAACTGTCATCCAGACATGTGTATCAGGGCCCACCTGTTTGAAGcgagaa                                                              >  6:279729/1‑130 (MQ=255)
                                                                              ccAGAGGCCCGGTTTTTGATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGCCCAACGCCCGGGAAAAAAACCAAAAAGTCATCCCGACGTGTGTATGCTGCCCGCGCGGTTTGAAGCTACAAACGCGCCTgg                                                    >  2:84115/1‑135 (MQ=255)
                                                                                         ttttttCATGGGGCATATGTTCGCAGTTAGCCCGCAAAACCCACAAACGCCCCGGAAAAAAACCAAACTGTCATCCACACATGTGTATACAGCCCGAGCGTTTTGAAACAACCAACCCGCCTGGCGCACACttt                                          >  8:50976/1‑134 (MQ=255)
                                                                                         ttttttCAGGGGGCATATGGGCGCTGTGCCCCCGAAAACCCCACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc                                        <  4:108865/136‑1 (MQ=255)
                                                                                         gtttttCATGGGGCATATGTGCGCGGGTAGCCCGCAAAACCGACAAACGCCCGGGAAAAAAACCAAACAGTCATCCAGACGTGTGGATCCTGCCCGCGCTTTTTGAAGCTACACACGCGCCGGGCGCAGACgttt                                         >  2:8591/2‑135 (MQ=255)
                                                                                                   gggCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACtt                              >  4:79201/1‑136 (MQ=255)
                                                                                                                cggTTCCCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAacgac                  <  7:129732/133‑1 (MQ=255)
                                                                                                                cTGTGCGCCCCAAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCAGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCTCACTTGCGAGGAacgac                  <  1:60749/135‑1 (MQ=255)
                                                                                                                 tGTTCGCCCGCAAAACCGACACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAACCGCGCCTGGCGCAGACTGTTCACTCACACTTGCGAGGAccga                   >  7:75110/1‑133 (MQ=255)
                                                                                                                 tGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAACCTGTCATCCAGACTTCTCTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTATTCACTCACACTTGCGAGGAacga                   >  5:30280/1‑133 (MQ=255)
                                                                                                                 tGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAacga                   >  7:87464/1‑133 (MQ=255)
                                                                                                                 tGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTGTTCACTCACACTTGCGAGGAacga                   >  1:356214/1‑133 (MQ=255)
                                                                                                                   ttCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTGTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAacga                   >  1:56288/1‑131 (MQ=255)
                                                                                                                    gcGCCCGCAAAACCCAAACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACgaga              <  6:115208/134‑1 (MQ=255)
                                                                                                                      gCCCACAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGAt             >  8:110366/1‑134 (MQ=255)
                                                                                                                       cccGAAAAACCCACACACCACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGc           <  8:85713/135‑1 (MQ=255)
                                                                                                                        ctcaaaaaCACACAAACGACAGGAAAAAAAACCAAACTGTCTTCAAGACTTGTGTATCCTGCCCGAGATTTTTGAAGCTACAAACGCGGCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGCa          <  1:169651/131‑1 (MQ=255)
                                                                                                                             aaaCCGAACAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTTCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGCATGaaa     <  7:63132/135‑1 (MQ=255)
                                                                                                                               aCCCACACACCACCGGTAAAAAAACCCAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGCATGAAATcc  <  4:5420/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCACACTTGCGAGGAACGACGAGATGCATGAAATCC  >  KT2440/5096834‑5097096

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: