Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,670,685 A→C 36.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,6850AC36.8% 21.8 / 2.5 38intergenic (+22/‑109)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/2);  new base C (12/3);  total (32/6)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCC  >  KT2440/1670567‑1670804
                                                                                                                       |                                                                                                                       
gTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGCCCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTGGGGGCAGCCGCGCACCCCCCCCGtc                                                                                                         >  1:133302/1‑135 (MQ=255)
  tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCCGCCAGCACGGCAACCGCGGCACGCCGCTGTAATTGCATGGCGTTTTGGGCTGGCGCGCCCCCCCAAccga                                                                                                         >  8:113149/1‑130 (MQ=255)
  tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCTGCAGCGCAGCCCCCACggc                                                                                                         >  3:261904/1‑131 (MQ=255)
  tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCACCGtctc                                                                                                       >  7:235751/1‑135 (MQ=255)
  tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGtc                                                                                                         >  8:103920/1‑133 (MQ=255)
  tCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGCCTTGGACGCCACGCCCGGCAACCCCGGCAGGCCGGTGTAATTGCTTGCCGTTTGGGGGTGCAGCGCAGCCCCACCggggc                                                                                                       >  1:49441/1‑131 (MQ=255)
    aGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAg                                                                                                     >  7:188409/1‑135 (MQ=255)
         tGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCtt                                                                                                 >  8:52320/1‑134 (MQ=255)
            agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACCTCGCCGCCgtttac                                                                                             >  8:275166/1‑133 (MQ=255)
              agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGGCTCAGCCGTTTGcct                                                                                           >  3:53424/1‑135 (MQ=255)
                             tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCATCCCCAACGTCGCAGCCTTTTACCTCCGGGCGATCAGaaa                                                                            >  8:286546/1‑135 (MQ=255)
                             tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGaa                                                                             >  1:49308/1‑134 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCTGCAGCGCAGCCCCAACGTCTCCTCCGTTTCCCACCTGGCTATCAGAAACAACCTAt                                                                    >  1:231952/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCACGTCTCAGCCGTTTCCCTCCGGTCTTTACGAAACCTCCTAt                                                                    >  7:134269/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                  >  3:149628/1‑135 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                    >  2:213115/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGTCTCAGCCTTTTACCGACTGTCTATAAGAAACCTCCTAt                                                                    >  1:66147/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGGCTCCGCCTTTTCCCTCCTGTCTATCAGCAACATCCTAt                                                                    >  8:195277/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGGCTCACCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                  >  1:157814/1‑135 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGTCGCAGCCGTTTCCCTCCTGTGGATCAGAAACATCCTATTc                                                                  >  5:107052/1‑135 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACCTCTCCGCCGTTTACCTCCAGTCTATCAGAAACATACTAATc                                                                  >  6:251776/1‑135 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCATCCTTTTACCTCCAGTGTATCAGAAACATCCTAtta                                                                  >  6:35543/1‑134 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                    >  8:171614/1‑133 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                  >  4:22003/1‑135 (MQ=255)
                                       ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCAAt                                                                    >  4:112851/1‑133 (MQ=255)
                                            aGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCt                                                             >  4:253346/1‑135 (MQ=255)
                                              gCTGTACCGGGCTTTGGCACCCAGCACGGAACCCGCGCAACGCCGGTTTAATTGCAGGGGGTTGGGGGCGGCCGCGCACCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGc                                                           <  8:134269/135‑1 (MQ=255)
                                                     aGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCAACGTCTCCGCCTTTTCCCTCCGGTCTATCAGCAACATCCTATTCAGCCGGCCGGgcgc                                                    >  6:166431/1‑135 (MQ=255)
                                                          tttGTCAGCCAGCACGGCAACCGGGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGc                                               >  4:122880/1‑135 (MQ=255)
                                                          tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGc                                               >  5:53671/1‑135 (MQ=255)
                                                               cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGGGTCCGGCGTTTGCCGCCGGTGTATCAGAAACCGACGATTCCACCTGCCGGGCGCCAAGccaaaa                                          >  2:265718/1‑130 (MQ=255)
                                                               cagccagcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCGGGCGCGCCGCCCCCACGGGTCAGCCGTTTTCCGCCCGTGTTTTCGAAAAAACCGATTACGGCCGCCGGGGGCCCAGCAaaaa                                          >  5:254270/1‑135 (MQ=255)
                                                                    agcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCCAGCATAAATTAGc                                     >  4:214215/1‑135 (MQ=255)
                                                                         cGGCACCCGCGGCAAGCCGGTGTAATTGCTGGGCGTTGGGGGCGGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                 <  8:235751/135‑1 (MQ=255)
                                                                              aCCGCGTCAAGCGGGTTTAATTGATTGGCTTTGGGGGCTGCTGCGAAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGTtcctcc                            <  4:109735/135‑1 (MQ=255)
                                                                                           gcggTAATGGCAGGGCGTGTGGGGCGGCAGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc               <  8:251846/132‑1 (MQ=255)
                                                                                                 aTTGCCTGGGGTTTGGGGCTGCCGCGCTGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATg          <  8:98963/134‑1 (MQ=255)
                                                                                                         ggttttGGGGCGGCCGCGCTCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGcc  <  7:103920/131‑1 (MQ=255)
                                                                                                                       |                                                                                                                       
GTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCC  >  KT2440/1670567‑1670804

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: