Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,670,685 A→C 41.8% 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,6850AC41.8% 10.3 / 2.9 34intergenic (+22/‑109)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/2);  new base C (10/5);  total (27/7)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATT  >  KT2440/1670571‑1670817
                                                                                                                   |                                                                                                                                    
aGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCAGCCCCAACGTCTCAg                                                                                                                  >  7:7186/1‑135 (MQ=255)
        agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTcc                                                                                                          >  6:130767/1‑135 (MQ=255)
               tGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCCGCCGTTTGCCTCCTGt                                                                                                   >  7:212081/1‑135 (MQ=255)
                       ccTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCAGCCCCCACCTCTCCGCCGTGTCCCTCCTGTCTATccg                                                                                            >  6:281394/1‑132 (MQ=255)
                        cgggCCAACGCCCTGACGGCTGTACTGGGCTTTGCCACCCAGCACGCAACCCGCGGCAAGCCGCTTTAATTGCAGGGGGTGT‑GGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGaa                                                                                          <  4:253124/133‑1 (MQ=255)
                            ccAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCCGCCCCAACGTCGCAGCCGTTTCCCGCCTGTCTATCAGAAACAt                                                                                      >  6:279723/1‑135 (MQ=255)
                              accGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTTCAGCGCAGCCCCAACATCTCAGCCTTTTGCCTCCTTTGTATAATCAACATcc                                                                                    >  6:238270/3‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑TGGGCGGCCGCGCAGCCCCCACGGCGCAGCCGTTTACCGGCTTTCGATCAGAAACATACTATTc                                                                               >  8:122542/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑TGGGCGGCCGCGCACCCCCCAAGGCTACGCCTTTTCCCACCAGTCGATCACAAACATCCTaat                                                                                >  2:185601/1‑132 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGGAGCGCAGCCCCCACGTCGCACCCTTTTGCCGCCGGTCTATCAGAAACATCCTAt                                                                                 >  5:104602/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCCCCCCCAACCTCGCAGCCATTTACCTCCGGTCTATCAGAAACagacaat                                                                                 >  5:62552/1‑127 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  7:281254/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCCGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  4:225444/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCCGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                               >  2:275952/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAtt                                                                                >  6:280491/1‑134 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCCGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  6:141775/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTACTGTCTATCAGAAACATCCTAt                                                                                 >  8:185001/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCCCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                               >  4:232330/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCGGGCCCGCAGCCCCAACGCACCAACCGTTTACCTCATGTCGATAAGTAACAACCTAATc                                                                               >  1:258737/1‑135 (MQ=255)
                                               aCAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                   >  6:321980/1‑135 (MQ=255)
                                                cAGGGCTTTGTCCGCCTGCACGGAACCCGCGGCAAGCCGCTGTAATTGCAGGGGGTTT‑GGGGGTGCACGGCAGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                   <  3:170105/134‑1 (MQ=255)
                                                 aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGAGTTT‑GGGGGTGGCGCGCCGCCCCCACGGCTCCGCCGTTTGCCGCCGGGGTATCAGCAACATCCGATTAAGCCTGCCGGgcgc                                                                 >  8:102163/1‑135 (MQ=255)
                                                        tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCa                                                           >  7:258167/1‑134 (MQ=255)
                                                        tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGGCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCa                                                           >  2:10257/1‑134 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGGAGCAGCGCCGCCCCCACCGCTCAGCCTTTTCCCTCCGGTCGATCAGAAACATCCTATTCAGCCTGCCCGGCGCCAAGCATaa                                                        >  6:91641/1‑134 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCCGCGCAGCCCCCACGGCGCAGCCTTTTCCCGCCTGTCTATCAGAAACATCCTATTCCGCCTGCCGGGCGCCCAGCCTaaa                                                       >  2:98044/1‑135 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCCACCCCCACCTCTCAGCCTTTTACCGCCGGTCGTTCCTCAACATCCTATTCAGCCTGCCCGGCGCCCAGcgataa                                                       >  5:268620/1‑130 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCGCAGCCGTTTCCCTGCTGTCGATAATAAACATTCTATTCCGCCTGCCGGGCGCCAAGCATaaa                                                       >  2:259761/1‑135 (MQ=255)
                                                                    cgcaaaCCGCGGCAAGCCGCTGAATTTGCTGGGCGTTG‑GGGGCGGCCGCGCACCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                              <  7:10128/131‑1 (MQ=255)
                                                                     ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAATCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtc                                             >  7:119471/1‑135 (MQ=255)
                                                                                cAACCCGCTGTATTTGCAGGGCGTTG‑GGGGGGCCCGCGCACCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAgg                                  <  4:274665/135‑1 (MQ=255)
                                                                                             ttGCATGGCGTTT‑GGGGGGCCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGc                      <  2:21952/134‑1 (MQ=255)
                                                                                             atGCTTGGGGTTTGGGGGCTGCCGCCCCCCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAACATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATg                       <  6:290061/133‑1 (MQ=255)
                                                                                                   ggCGGTG‑GGGGCGGCCGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcc               <  2:134668/135‑1 (MQ=255)
                                                                                                                 gcagcccagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAtt  <  3:234455/135‑1 (MQ=255)
                                                                                                                   |                                                                                                                                    
AGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTT‑GGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATT  >  KT2440/1670571‑1670817

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: ‑