Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,670,685 A→C 27.7% intergenic (+22/‑109) thrC_2 → / → CJ019_01569 Threonine synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,670,6850AC27.7% 36.1 / 2.7 40intergenic (+22/‑109)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/4);  new base C (5/6);  total (28/12)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGT  >  KT2440/1670554‑1670812
                                                                                                                                    |                                                                                                                               
aCCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTgcagcg                                                                                                                              >  2:171794/1‑135 (MQ=255)
      gCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGccc                                                                                                                        >  3:274141/1‑135 (MQ=255)
          ccTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGGAGCGCAGCCCCAAc                                                                                                                    >  4:15361/1‑135 (MQ=255)
               tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtc                                                                                                                 >  3:96702/1‑133 (MQ=255)
                           agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCCCGCCCCCCCAACCTCTCACCCatttaacg                                                                                                   >  6:220771/1‑131 (MQ=255)
                                tGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGt                                                                                              >  3:178450/1‑135 (MQ=255)
                                ggCACGGGGCTGGAAAACTAGCTCAAGGCTGTACAGGGATTTGACAGCCTGCATGCAAACCGCGCCAAGCCGTTGTATTTGCTTGTCTTTTTGGGCGGCCTCGCTGCCCCAACGTCTCAGCCTTTTCCCTCCTGTc                                                                                             <  8:79998/135‑1 (MQ=255)
                                      gtcggggCAAACCACCTGATGCTTGTACAGGGTTTTGTCACCCTGCCGGCAACCCGCGTAACGCGTCTTTATTTGGATGTGTTTTGGGGGGGCTGGGCACCCAAAACGTCTCAGCCTTTTCCCTCCTGTCTATCa                                                                                        <  2:223529/130‑1 (MQ=255)
                                          tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGaa                                                                                     >  2:185674/1‑134 (MQ=255)
                                          tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGaa                                                                                     >  4:35841/1‑134 (MQ=255)
                                             ccAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCAGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTGTATCAGAAACAt                                                                                 >  3:271703/1‑135 (MQ=255)
                                               caCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCCCAGCCCCCACGTCTCAGCGTTTTTCCTCCTGTGTATAAGAAACAtac                                                                               >  7:232739/2‑133 (MQ=255)
                                               aaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCACCGCCTTTTACCTCCTGTCTATAAGAAACATc                                                                                >  3:68269/1‑134 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGGCGCAGCCGTTTGCCTCCGGGCTATGAGAAACCTCCTAtt                                                                           >  7:114372/1‑134 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCCGCCTTTTCCCGCCTGTCTATCAGAAACAGCCTATTCa                                                                         >  7:21394/1‑136 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAtt                                                                           >  8:40690/1‑134 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAtt                                                                           >  8:207967/1‑134 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                            >  4:33005/1‑133 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                          >  5:72384/1‑135 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATAAGAAACATCCTATTc                                                                          >  5:49719/1‑135 (MQ=255)
                                                    ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGCCGCGCCGCCCCCCCCTGTCAGCCCTTTCCCGCCGGGCTATCAACAACAACCGATTc                                                                          >  5:274474/1‑135 (MQ=255)
                                                      tGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTAAGc                                                                       >  7:30057/1‑136 (MQ=255)
                                                                  aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCACCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGcc                                                           >  7:167056/1‑136 (MQ=255)
                                                                       tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGCCGCCCCCCCCCCCCGGGGCAGCCGTTTACCTCCCGTCGATAACAAACAAACGATTAAGCCGGCCCGGCGcccaa                                                        >  1:272711/1‑131 (MQ=255)
                                                                         tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTGCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCAt                                                     >  7:167369/1‑135 (MQ=255)
                                                                         tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCAt                                                     >  8:98529/1‑135 (MQ=255)
                                                                            cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCGATCAGAAACAACCTATTCAGCCTGCCGGGCCCCAAGCATaaa                                                  >  6:246717/1‑135 (MQ=255)
                                                                            cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCATTTCCCGCCTGTGTATGAGAAACATCCAATTCCCCCAGCCGGGCGCCAAGCATaaaa                                                 >  7:252113/1‑135 (MQ=255)
                                                                            cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCAGCAGCGCAGCCCCAACGTCTCAGCCATTACCCTCATGTCTATCAGAAACATCCTAATCACCCAGCCGGGCGCCAAGCATaa                                                   >  4:25128/1‑134 (MQ=255)
                                                                                     cGGCACCGGCGGCAAGCCGTTGTAATGGCAGGGGGTTGGGGGGGGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                         <  4:121929/135‑1 (MQ=255)
                                                                                     cGGCACCCGGGCCACGCCGCTGTAATGCCAGGGGGTGGGGGGCGGCGGCGCCGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                         <  1:171794/135‑1 (MQ=255)
                                                                                      ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCGCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtc                                        >  6:47351/1‑135 (MQ=255)
                                                                                                ggaacgCCGCTTTAATTGCATGGGGTTGGGGGCGCCAGCGCCCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                               <  6:49719/130‑1 (MQ=255)
                                                                                                       ggCTGTAATTGCTTGGGGTGTGGGGCGGCCGCGCCGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                       <  8:252113/135‑1 (MQ=255)
                                                                                                            aaaTGGCAGGGCGTGTGGGGCTCCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCAt                   <  7:93576/134‑1 (MQ=255)
                                                                                                             aaTGGCAGGGCGTTGGGGGCGGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCAt                   <  2:11353/134‑1 (MQ=255)
                                                                                                               ttGCATGGCTTTGGGGGTTCCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCa                <  3:207217/135‑1 (MQ=255)
                                                                                                               ttGCAGGGGGTTGGGGGCGGCCGCGCCCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCa                <  4:271703/135‑1 (MQ=255)
                                                                                                                       cGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCGGTCTTTCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcccg        >  3:118012/1‑134 (MQ=255)
                                                                                                                         tttGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt      <  4:178450/135‑1 (MQ=255)
                                                                                                                             gggCTGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATgtgt  <  3:33005/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                               
ACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGT  >  KT2440/1670554‑1670812

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: