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

AGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATT  >  KT2440/1670571‑1670817
                                                                                                                   |                                                                                                                                    
aGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCAg                                                                                                                  >  7:7186/1‑135 (MQ=255)
        agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTcc                                                                                                          >  6:130767/1‑135 (MQ=255)
               tGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCGTTTGCCTCCTGt                                                                                                   >  7:212081/1‑135 (MQ=255)
                       ccTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACCTCTCCGCCGTGTCCCTCCTGTCTATccg                                                                                            >  6:281394/1‑132 (MQ=255)
                        cgggCCAACGCCCTGACGGCTGTACTGGGCTTTGCCACCCAGCACGCAACCCGCGGCAAGCCGCTTTAATTGCAGGGGGTGTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGaa                                                                                          <  4:253124/133‑1 (MQ=255)
                            ccAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCAACGTCGCAGCCGTTTCCCGCCTGTCTATCAGAAACAt                                                                                      >  6:279723/1‑135 (MQ=255)
                              accGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTTCAGCGCAGCCCCAACATCTCAGCCTTTTGCCTCCTTTGTATAATCAACATcc                                                                                    >  6:238270/3‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCGGCCGCGCAGCCCCCACGGCGCAGCCGTTTACCGGCTTTCGATCAGAAACATACTATTc                                                                               >  8:122542/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCGGCCGCGCACCCCCCAAGGCTACGCCTTTTCCCACCAGTCGATCACAAACATCCTaat                                                                                >  2:185601/1‑132 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGGAGCGCAGCCCCCACGTCGCACCCTTTTGCCGCCGGTCTATCAGAAACATCCTAt                                                                                 >  5:104602/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCCCCCCAACCTCGCAGCCATTTACCTCCGGTCTATCAGAAACagacaat                                                                                 >  5:62552/1‑127 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  7:281254/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  4:225444/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                               >  2:275952/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAtt                                                                                >  6:280491/1‑134 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                 >  6:141775/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTACTGTCTATCAGAAACATCCTAt                                                                                 >  8:185001/1‑133 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCCCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                               >  4:232330/1‑135 (MQ=255)
                                   ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGGCCCGCAGCCCCAACGCACCAACCGTTTACCTCATGTCGATAAGTAACAACCTAATc                                                                               >  1:258737/1‑135 (MQ=255)
                                               aCAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                   >  6:321980/1‑135 (MQ=255)
                                                cAGGGCTTTGTCCGCCTGCACGGAACCCGCGGCAAGCCGCTGTAATTGCAGGGGGTTTGGGGGTGCACGGCAGCCCAAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                   <  3:170105/134‑1 (MQ=255)
                                                 aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGAGTTTGGGGGTGGCGCGCCGCCCCCACGGCTCCGCCGTTTGCCGCCGGGGTATCAGCAACATCCGATTAAGCCTGCCGGgcgc                                                                 >  8:102163/1‑135 (MQ=255)
                                                        tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCa                                                           >  7:258167/1‑134 (MQ=255)
                                                        tGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGGCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCa                                                           >  2:10257/1‑134 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGAGCAGCGCCGCCCCCACCGCTCAGCCTTTTCCCTCCGGTCGATCAGAAACATCCTATTCAGCCTGCCCGGCGCCAAGCATaa                                                        >  6:91641/1‑134 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGGCGCAGCCTTTTCCCGCCTGTCTATCAGAAACATCCTATTCCGCCTGCCGGGCGCCCAGCCTaaa                                                       >  2:98044/1‑135 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCACCCCCACCTCTCAGCCTTTTACCGCCGGTCGTTCCTCAACATCCTATTCAGCCTGCCCGGCGCCCAGcgataa                                                       >  5:268620/1‑130 (MQ=255)
                                                           cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCGCAGCCGTTTCCCTGCTGTCGATAATAAACATTCTATTCCGCCTGCCGGGCGCCAAGCATaaa                                                       >  2:259761/1‑135 (MQ=255)
                                                                    cgcaaaCCGCGGCAAGCCGCTGAATTTGCTGGGCGTTGGGGGCGGCCGCGCACCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                              <  7:10128/131‑1 (MQ=255)
                                                                     ggCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAATCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtc                                             >  7:119471/1‑135 (MQ=255)
                                                                                cAACCCGCTGTATTTGCAGGGCGTTGGGGGGGCCCGCGCACCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAgg                                  <  4:274665/135‑1 (MQ=255)
                                                                                             ttGCATGGCGTTTGGGGGGCCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGc                      <  2:21952/134‑1 (MQ=255)
                                                                                             atGCTTGGGGTTTGGGGGCTGCCGCCCCCCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAACATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATg                       <  6:290061/133‑1 (MQ=255)
                                                                                                   ggCGGTGGGGGCGGCCGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcc               <  2:134668/135‑1 (MQ=255)
                                                                                                                 gcagcccagcCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGAtt  <  3:234455/135‑1 (MQ=255)
                                                                                                                   |                                                                                                                                    
AGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATT  >  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: