Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,266,383 T→G 49.1% intergenic (+121/+47) cysK_1 → / ← CJ019_01185 Cysteine synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,266,3830TG49.1% 2.6 / 3.8 29intergenic (+121/+47)cysK_1/CJ019_01185Cysteine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base G (10/3);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGTCAGGCTCATGATCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGCGTTGCG  >  KT2440/1266268‑1266513
                                                                                                                   |                                                                                                                                  
tgtgTCAGGCTCATGATCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCAt                                                                                                                 >  6:82343/1‑135 (MQ=255)
tgtgTCAGGCTCATGATCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATAGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCCCACTGTTAGAGccca                                                                                                                 >  2:236079/1‑133 (MQ=25)
            aTGATCCAGAGGCTGGCACGGAAACTGTGGGAGAGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTAGAGCCATCTCTCAGAatt                                                                                                      >  8:212053/1‑134 (MQ=255)
               aTCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCTCGGGCGCGCCCGCTCCCACAGTGTTAGAGCCAGCTCTCAGAttttttt                                                                                                  >  7:225828/1‑133 (MQ=255)
               aTCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGGGCGCCCGCGCCCCCAGAGTTCGCGCCATCTCTCACAatttttt                                                                                                  >  7:316119/1‑133 (MQ=25)
               aTCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTAt                                                                                                  >  8:383868/1‑135 (MQ=255)
               aTCCAGAGGCTGGCACGGAAACTGTGGGAGAGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCGCCCACAGTGTTCGCGCCATCTCTCAGAttttt                                                                                                    >  3:53444/1‑133 (MQ=255)
                     agGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCGCCCCCCGGGTTCGCGCCCTCTCGCAGGtgtttttgggg                                                                                              >  6:281914/1‑127 (MQ=18)
                       gCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTGGCGGGCGCGCCCGCGCCCCCAGGGGTCGAGCCCGCGCGCAGCGTTTTATggcggggg                                                                                          >  3:125430/1‑131 (MQ=18)
                       gCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCGCCCCCAGGGTTGGTGCCCTCGCTCAGCTTTTTATGGGGTCAg                                                                                          >  1:44303/1‑135 (MQ=21)
                       gCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCCCGCCCCCCGGGGTGGGGGCCGCTCGGAGATTTTTATAGCGTCAg                                                                                          >  5:175484/1‑135 (MQ=21)
                             aCGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGGGCGCCCCCGGCCACACTGTTGGCGCCCTCGCTCAGATTTTTATGGGGTTCGTCGa                                                                                      >  8:266216/1‑133 (MQ=21)
                               ggAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCGCCGCCACTGTTCGAGCCAACGCTCAGAATTTTATAGCGGCAGTCGATccg                                                                                  >  8:377329/1‑133 (MQ=25)
                                        gggAGCGGGCACGCCCGCGAAGCTGGCGACGCGGTGGAGGGCACGGGCTTCGCCCGTGTGCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTgg                                                                         <  6:384146/135‑1 (MQ=255)
                                                gCACGCCCGCTACGCCGGCGACGCGGTGGATGGCCCGGGCTGCGCCGGGGTGCGCGGGCGCGCCCGCGCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCActtctt                                                                 <  4:89930/135‑1 (MQ=255)
                                                   aGCCCGGGAAGCCGGCGACGCGGTGGATGGCACGGGCGTCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTa                                                              <  8:297157/135‑1 (MQ=255)
                                                                              tGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATcgcg                                                                                              >  8:303772/1‑76 (MQ=255)
                                                                              tGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATcgcg                                                                                              <  7:303772/76‑1 (MQ=255)
                                                                                gCACGGGCTGCGCCCGGGGTCGCGGGGGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGt                                 <  4:82162/135‑1 (MQ=255)
                                                                                       cTTCGCCCGTGTTCGCGGGCGCGCCCGCGGCCCCAGGGTTGGAGGCCTCTCTCCGCTTTTTATAGGGGGCGGGGGTCAGGCGCCGGGGGCCGTGGTGGACGTATTATCGCTACCGCCCGGGCCCGGGGCGTTTCg                          >  6:200833/1‑135 (MQ=255)
                                                                                        ggcGCCGGTGTGGGGGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACgg                         <  5:200833/133‑1 (MQ=255)
                                                                                           gCCGGTGTGCGGGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGcc                       <  2:57624/134‑1 (MQ=255)
                                                                                           gCCCGTGTTCGCGGGCGCGCCCGCGCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGACa                      <  8:137808/135‑1 (MQ=255)
                                                                                                   tCGCGGGCGCGCCCGCTCCCCCAGGGTTCGCGCCCTCGCGCAGAGTTTTATCGCGGCCGTCGCTCAGGCGCCGGGGCCCGGCTTGTACTTCTTCTCGCTTCCGCACCGGCACGGGTCGGTACGGCCGACGTCCt               >  6:120982/1‑134 (MQ=255)
                                                                                                         gcgcgcCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGc        <  8:215016/135‑1 (MQ=255)
                                                                                                          cgcgcCCGCCCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGCg       <  1:71362/135‑1 (MQ=255)
                                                                                                               ccGCGCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGCGTTGCg  <  6:254736/135‑1 (MQ=255)
                                                                                                               ccGCGCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATAGCGTCAGTCGATAAGGCGCCGTGGCACTTATTTTTCGTCATCTAGCAAACGCCAGTGGACCGGTCGTTAAGGCCCACCTCAGACAGGGCGTTggg  >  8:250990/1‑133 (MQ=255)
                                                                                                               ccGCCCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGCGTTGCg  <  5:47146/135‑1 (MQ=255)
                                                                                                                   |                                                                                                                                  
TGTGTCAGGCTCATGATCCAGAGGCTGGCACGGAAACTGTGGGAGCGGGCAAGCCCGCGAAGCAGGCGACGCGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTTCGAGCCATCTCTCAGATTTTTATCGCGTCAGTCGATCAGGCGCCGTGGCACTTCTTGTACTTCTTCTCGCTACCGCACGGGCACGGGTCGTTACGGCCAACGTCCTTCAGGGCGTTGCG  >  KT2440/1266268‑1266513

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: