Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,026,996 C→T 12.1% intergenic (+335/‑75) yabJ_3 → / → CJ019_03742 2‑iminobutanoate/2‑iminopropanoate deaminase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,026,9960CT12.1% 91.1 / 4.5 34intergenic (+335/‑75)yabJ_3/CJ019_037422‑iminobutanoate/2‑iminopropanoate deaminase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base T (2/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

ACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGT  >  KT2440/5026862‑5027101
                                                                                                                                      |                                                                                                         
acaGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGCCGCTCGCCAGGCATCACCATCGTAAGCGTCCGGCGAACACCCCATcc                                                                                                           >  2:359049/1‑135 (MQ=255)
  aGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGa                                                                                                         >  3:197916/1‑135 (MQ=255)
   gggggTCGCCGACGGGCTCACTCAGTGGGGCTTGGTGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCGCGTTGGACGCTCGCCAGGAATGACCATCGTAAGGGGGCGGCGAACTCACCTTCCGaa                                                                                                        <  8:381060/132‑1 (MQ=255)
   gCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGaa                                                                                                        >  2:100811/1‑135 (MQ=255)
    cTGGTCGCCGACCGGCTCACTCAGTGGGGCTTTGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAAc                                                                                                       <  7:55842/135‑1 (MQ=255)
    cTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAAc                                                                                                       <  7:211161/135‑1 (MQ=255)
          gCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGCCGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTACGAACTCCAg                                                                                                  >  5:64688/1‑134 (MQ=255)
          gCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAg                                                                                                  >  8:193618/1‑134 (MQ=255)
          gCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGc                                                                                                 <  4:31722/135‑1 (MQ=255)
                  gCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAggg                                                                                         <  5:80001/135‑1 (MQ=255)
                  gCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAggg                                                                                         <  7:193618/135‑1 (MQ=255)
                             ggggCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTcc                                                                              <  1:319776/135‑1 (MQ=255)
                               ggCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCTGAACTCCAGCATTAAGGGCGATGGTGTCcgc                                                                            >  8:252771/1‑135 (MQ=255)
                               ggCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCcgc                                                                            >  6:362241/1‑135 (MQ=255)
                               ggCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAAACACCATCGTAAGCGTCCGGCGAACTCACCTTCTGAACTCCAGCATTAAGGGCGATGGTGTCcgc                                                                            >  8:258995/1‑135 (MQ=255)
                                    ggAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAtt                                                                       >  3:329436/1‑135 (MQ=255)
                                       gaTTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTTGACGCTCGCCAGGAATCACCATAGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTACGCGTAtttt                                                                     >  6:167458/1‑134 (MQ=255)
                                       gaTTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAtttt                                                                     >  3:112041/1‑134 (MQ=255)
                                            aCCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGc                                                                <  4:197916/134‑1 (MQ=255)
                                             ccGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCgg                                                              >  3:297069/1‑135 (MQ=255)
                                                 ggACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGAcgc                                                          <  8:176463/135‑1 (MQ=255)
                                                        cAAAAACCGGGGGCGTTGGAACGCTCGGTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCAcc                                                   <  6:305485/135‑1 (MQ=255)
                                                          aaaaCCGGGGTCGTTGGAACGCTCGTTGGAAGCTCGCCAGGAATCAACATCGTAAGCGTCCGGCGAACTCACCTTCAGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCTCTCAccc                                                  >  6:361428/1‑134 (MQ=21)
                                                                    tCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGCTCACCCTTATCGCCa                                         >  1:192803/1‑133 (MQ=255)
                                                                                  gtgggACGCGCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCa                         <  1:359049/132‑1 (MQ=255)
                                                                                        cGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAg                   <  6:64688/135‑1 (MQ=255)
                                                                                        cGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAg                   <  5:362241/135‑1 (MQ=255)
                                                                                        cGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAg                   <  3:96102/135‑1 (MQ=255)
                                                                                               cAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTAt            >  2:317588/1‑135 (MQ=40)
                                                                                                        ccATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCg   >  3:198468/1‑135 (MQ=33)
                                                                                                         cATCGTAAGCGTCCGGCGAACTCACCTTCTGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGt  <  7:252771/135‑1 (MQ=17)
                                                                                                         cATCGTAAGCGTCCGGCGAACTCACCTTCTGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGt  <  7:258995/135‑1 (MQ=17)
                                                                                                         cATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGt  <  1:73702/135‑1 (MQ=32)
                                                                                                         cATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGt  <  1:100811/135‑1 (MQ=32)
                                                                                                          aTCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGt  >  1:47144/1‑134 (MQ=31)
                                                                                                                                      |                                                                                                         
ACAGCTGGTCGCCGACCGGCTCACTCAGTGGGGCTTGGAGATTCACCGTGGACTGGCAAAAACCGGGGTCGTTGGAACGCTCGTTGGACGCTCGCCAGGAATCACCATCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGT  >  KT2440/5026862‑5027101

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