Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,805,741 A→T 52.2% intergenic (+161/‑588) dsbG → / → CJ019_04513 Thiol:disulfide interchange protein DsbG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,805,7410AT52.2% ‑2.2 / 3.3 15intergenic (+161/‑588)dsbG/CJ019_04513Thiol:disulfide interchange protein DsbG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (3/3);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGT  >  KT2440/4805607‑4805869
                                                                                                                                      |                                                                                                                                  
gcCGTGAAGCAGGCGACGTGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGCGAACCCGCTCCTACGGGGCGCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTa                                                                                                                                    <  1:127697/135‑1 (MQ=255)
     gAAGCAGGCGACGGGGGGCATGCCACCGGCACCGCCGGTGTGCGCGGGTCAACCCGCTCCTACGGGGCGCGCGCTAACGGCGGGGTTTTCCGCCAGGCCCGTGAGCGCAAGGCCTGCGCGGTCGGCCGTAGGAGCg                                                                                                                            <  8:137536/136‑1 (MQ=255)
      aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTTGGcgcgc                                                                                                                           >  8:68025/1‑135 (MQ=255)
      aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTAGGGGCgg                                                                                                                           >  8:423978/1‑136 (MQ=255)
      aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCCTacggcgc                                                                                                                              >  2:16490/1‑129 (MQ=255)
               aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTGGGGGCGGGGTTTAcc                                                                                                                    >  3:113745/1‑135 (MQ=255)
                            cACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTTGGCGCCGGGTTTACCCCGCAAAGGCCCg                                                                                                       >  4:371430/1‑135 (MQ=38)
                                         ggTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTTGGCGCCGGGTTTTCCCCGGAAAGGGCCGGCCCGGGGGGCCa                                                                                          >  4:68438/1‑135 (MQ=255)
                                         ggTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTGGGGGCGGGGTTTCCCCCGAAACGGCCGGGCCGGCggggcca                                                                                         >  4:348301/1‑131 (MQ=255)
                                                                                         ggttggcgcCCGGCCGGTGAGCGCACGCCCGGCGCGGGCGGCCGTAGGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                         <  1:152053/130‑1 (MQ=17)
                                                                                         gTTTTGCGCCAGGCCGGTGATCGAACGGCCTGGGCGGGCGGCCGTAGGCGCGGGTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                         <  2:267829/135‑1 (MQ=17)
                                                                                                  cAGCCCGGTGAGCGAACGCCCGGCGGGGGCGGCGGAAGCGGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                               <  8:364365/136‑1 (MQ=14)
                                                                                                                 atggCCTGCGGGGGCGGCGGTTGTCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCggggg                 <  6:184971/133‑1 (MQ=17)
                                                                                                                        gcggGGGCGGCCGTTGGTGGGGTTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGac          <  3:68438/135‑1 (MQ=18)
                                                                                                                               ggcccGTTGTAGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGt  <  7:72786/133‑1 (MQ=255)
                                                                                                                                 gCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                    >  7:110376/1‑84 (MQ=255)
                                                                                                                                 gCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                    <  8:110376/84‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
GCCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGT  >  KT2440/4805607‑4805869

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: