Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,670,690 A→C 34.8% intergenic (+27/‑104) 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,6900AC34.8% 17.9 / 2.4 35intergenic (+27/‑104)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/3);  new base C (10/3);  total (28/7)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCA  >  KT2440/1670569‑1670819
                                                                                                                          |                                                                                                                                 
tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGtc                                                                                                                        >  8:96099/1‑133 (MQ=255)
       tGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGGGCCGCGCCGCCCCCACGGCTCAGCCgtt                                                                                                               >  5:20692/1‑135 (MQ=255)
            agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCACCTCTCAGCCGTTTccat                                                                                                          >  4:229824/1‑133 (MQ=255)
                   gccGGGCCCGGCCAACCCCCGTGACGGCTTTACTGGGTTTGGCCCCCCCGCAGGCCAACCGCGGAACGCCGCTTAAATGCCAGGGGGTTGGGGGCGGCCGCGCACCCCCCACGTCTCAGCCTTTTCCCTCCTGTc                                                                                                    <  6:108342/133‑1 (MQ=255)
                              cccaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAAcct                                                                                        >  8:129617/4‑133 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAATCCGCTGTAATTGCATGGCGTTTGGGGCTTCACCGCAGCCCCACCGTCTCAGCCTTTTCCCTACAGTCTATCAGAAACAtacaat                                                                                   >  2:283961/1‑128 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGGTCCGCGCCGCCCCCCCGTGTCAGCCGGTTTCCTCCGGTCGATCGGGAACAGACTAtt                                                                                  >  5:42921/1‑134 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCACCCTTTTCCCTCCGGTCTATAAGAAAAAACCTATTc                                                                                 >  8:344574/1‑135 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGGCGCAGCCTTTTACCTCCTGTCTATCAGAAACATCCTAt                                                                                   >  8:339054/1‑133 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCGTTTCCCTCCGGTCTATCAGAAACATCCTAt                                                                                   >  2:93435/1‑133 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACCTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                   >  2:344750/1‑133 (MQ=255)
                                     ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACCGCTCAGCCTGTTCCCTCCTGTCGATCAGAAACATCCTAtt                                                                                  >  1:63850/1‑134 (MQ=255)
                                                  aCAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTGTCAGCCTTTTCCCTCCAGGGGTTAAGAAACATCCTATTAACGCTGCCCGgc                                                                     >  6:144594/1‑135 (MQ=255)
                                                   cAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGGAGCGCCGCCCCCACCGGTCAGCCGTTTGCCTCCGGTCTATCAGAAACCACCAATTAAGCCGTCCCGgcg                                                                    >  5:122259/1‑135 (MQ=255)
                                                    aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCAACCTCTACGCCGTTTCCCTACTGTCTATCAGAAACAAACTATTCAGCCTGCCGGgcgc                                                                   >  6:359884/1‑135 (MQ=255)
                                                    aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTATCAGCCTTTTCCCTCCGGTCTATCAGCAACATCCTATTCAGCCTGCCGGgcgc                                                                   >  8:369954/1‑135 (MQ=255)
                                                      ggCTTTGTCAGCCAGCACGGCAACCGCGGAAAGCCGCTTTAATTGCAGGGCGTTGGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCa                                                                 <  4:87780/135‑1 (MQ=255)
                                                       gCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAAAGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCaa                                                                >  4:91922/1‑135 (MQ=255)
                                                         tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGTGCCGCGCCGCCCCAACGTCTCAGCCGTTTCCCGCCAGGATTTAAGTAACAACCGATAAAGCCTGCCCGgcggccag                                                               >  1:353712/1‑129 (MQ=255)
                                                           tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCGATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCa                                                             >  7:52548/1‑134 (MQ=255)
                                                             tCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaa                                                          >  1:35025/1‑135 (MQ=255)
                                                              cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTTCATCGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCCCCCTTTTCCCGCCTGTCTATAAGAAACAGCCAATACAGCCTGCCGGGCGCCAAGCCTaaa                                                         >  8:11211/1‑135 (MQ=255)
                                                              cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaaa                                                         >  4:76558/1‑135 (MQ=255)
                                                              cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGGCTCCGCCGTTTCCCTACGGTCTATCAGAAACAACCGATTCCGCCAGCCGGGCGCCAACCATaaa                                                         >  7:342534/1‑135 (MQ=255)
                                                              cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATaaa                                                         >  3:361778/1‑135 (MQ=255)
                                                                   agcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTACCTCCTGTCTATAAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                    >  3:87780/1‑135 (MQ=255)
                                                                        ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGTGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtc                                               >  4:276731/1‑135 (MQ=255)
                                                                                 ggaacgCCGCTGTAATGGCGGGCCTGGTGGGCCTGCGGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                                      <  2:273858/130‑1 (MQ=255)
                                                                                   cacgCCGCTGAAATGGATTGGGGTTGGGGCCTGCTGCGACCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAgg                                    <  7:101721/132‑1 (MQ=255)
                                                                                   cAAGCGGCTTTATTTCCAGGGCGTTGGGGGCTGAAGCGATGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAgg                                    <  2:43456/135‑1 (MQ=255)
                                                                                   cAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCa                                      >  2:325371/1‑133 (MQ=255)
                                                                                              tatggCAGGGCTTTGGGGGCGCCCGCGCCCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCAt                          <  6:122259/130‑1 (MQ=255)
                                                                                               aTTTCATGGGGTGTGGGGCGCGACCGCAGCCCAAACGTCTCTGACTTTTCCCTCCTGTGACTCAGTAACTTCCTATTAAGGGTGCAGGGCGGCAAGCATAATTTAGCCGGTCCTCCCATCTGGAGCTGGTCATg                         <  8:63244/134‑1 (MQ=255)
                                                                                                        cGTTTGGGGCTGCAGCGCCGCCCCAACGGCTCAGCCTTTTGCCTCCGGTCTATCAGAAACATCCTATTCAGCCGGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCCTCAGGCGCGGCTCATGCACCCGccc                >  2:18414/1‑134 (MQ=255)
                                                                                                        cGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCCCCAGCCCt               >  7:363123/1‑135 (MQ=255)
                                                                                                                     agcgcagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCa  >  8:328240/1‑135 (MQ=255)
                                                                                                                          |                                                                                                                                 
TCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCA  >  KT2440/1670569‑1670819

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: