Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,491,323 T→G 61.9% intergenic (‑685/+36) CJ019_03306 ← / ← CJ019_03307 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,491,3230TG61.9% 0.2 / 6.2 34intergenic (‑685/+36)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/8);  new base G (15/6);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 8.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCG  >  KT2440/3491202‑3491454
                                                                                                                            |                                                                                                                                    
cTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCGGCCCGCGAACAcc                                                                                                                        >  7:42080/1‑135 (MQ=255)
 tACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGGCGCCGCCCCGCCAACAccc                                                                                                                       >  3:10378/1‑134 (MQ=21)
      cGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGCACACCCGCGCa                                                                                                                  >  7:55332/1‑135 (MQ=255)
      cGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGGACCCCCgcgcc                                                                                                                  >  3:110621/1‑134 (MQ=21)
         ccGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGAACACCCGCgccgcc                                                                                                               >  2:80330/1‑135 (MQ=255)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCGCGCACCCCCGCgccgcccg                                                                                                             >  7:96975/1‑133 (MQ=18)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACCCCCGCGCCGCCgg                                                                                                             >  7:251506/1‑135 (MQ=21)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACCCCCGCGCAGCCgg                                                                                                             >  8:198288/1‑135 (MQ=21)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGCACACCGGCGCCCCCgg                                                                                                             >  2:18331/1‑135 (MQ=255)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGCACACCCGCgccgcccg                                                                                                             >  4:111753/1‑133 (MQ=21)
           ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGCACACCCGCGCAGcccg                                                                                                             >  6:111557/1‑133 (MQ=255)
           ggtggtGTCGCCACGATGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGGCGCCGGGCCGCGCACCCCCGCgccgcccg                                                                                                             >  7:141430/1‑133 (MQ=18)
                tgtCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGAACACCCGCGCCGCCGGGTTCa                                                                                                        >  8:32880/1‑135 (MQ=21)
                     ccACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACACCCGCGCACCCGGGTTCATGCa                                                                                                    >  3:298096/1‑134 (MQ=18)
                     ccACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTTCCCGCCGAAACCCGCCCACCCGCtttagaaca                                                                                                    >  1:185165/1‑128 (MQ=17)
                        ggTTGAAGGAGCTGTTTAACCCGCTATTTGGCGGGACCTGCGGGCAATTAACCGGTGTCTTTCACCGCCACCGCCGCCAGCGCATTCGGGGGCTGCGTGCCCGCTAACACCGGCGCAGCCGGTTTCATCCAACTc                                                                                                <  6:109413/134‑1 (MQ=255)
                                  gCTGGTTTCCCCGCTATTAGGCCGGACCTGCCGCCAACTATCCCTGGTCATTCCCCGCCACCCCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATc                                                                                      <  7:32880/135‑1 (MQ=37)
                                  gCTGGTTTACCCGCCAATAGGCCGGACCGGCCGGCCACTATCCCTGGACATTCACCGTCACCGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATc                                                                                      <  5:111557/135‑1 (MQ=255)
                                    tgtttgaCCCGCTAATTGGCGGCCCTGGCGGGCCAATCTCCGTGTTCATACCCGGCCCCCGCCCCCCCCGCTTTCGGGGGCAGCCGCCCCGGGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGa                                                                                    <  8:131359/129‑1 (MQ=255)
                                      gTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCCCGCACACCCGCGCCCCCCCGGTCATCCAACTAAACCACCCTCGccc                                                                                  >  7:145319/1‑135 (MQ=17)
                                       tttATCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGAACACCGGCGCCGCCGCGTTCCTCCAACTCACCCCCCATCGAccg                                                                                 >  4:228622/1‑134 (MQ=21)
                                              gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCATCCTGCCCGCGCACACCGGCGCAGCCGGTTTGCTACAACTAAAAGGCCAAAGACCTTTAccg                                                                           >  8:321860/1‑132 (MQ=18)
                                                       ccGGACCTGCCGGCCAATACCCCTGGTCATTCCCCGACACCGCCCCCCCCGCATGCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTctct                                                                 <  8:96975/135‑1 (MQ=38)
                                                        gggtccTGCCGGCCAAATAGCCGTGGTCTATCCCGGTCCCCGCCCCCCCCCATTGCGGGGGCCGCCGCCCCGCCAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTctct                                                                 <  2:293653/131‑1 (MQ=255)
                                                                                     atccaccgccaccgcCGCCTCCGCCTTGGCGGGCAGCCTGCCGCGTAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTg                                     <  2:71526/132‑1 (MQ=255)
                                                                                       tCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCGCGCACACCGGCGCCGCCGGGTTGATCCAACTCCCCCACCCTCGCCCGGGACGCCAACCTCTCTGGCAACGATTGGCGTGGCTGGCCCATTGtg                                   >  2:232076/1‑134 (MQ=255)
                                                                                         aCCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGCACACCCGCGCAGCCGGTTTCATACAACTCAAACCCCAACGACCTGTTAGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTTTCGt                                 >  8:179103/1‑135 (MQ=255)
                                                                                             gcacaGCCGCCACCGCATGCGCGGGCAGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGAt                             <  4:110621/134‑1 (MQ=255)
                                                                                                  gccgccCCCGCATGCGGGGGCCGCCGCCCGGCTAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                        <  6:284290/135‑1 (MQ=255)
                                                                                                        cccGCATTCGGGGGAAGCCGCCCCGCCAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                  <  8:55332/134‑1 (MQ=255)
                                                                                                        aCCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                  <  7:209449/135‑1 (MQ=255)
                                                                                                        aCCGCATTCGCGGGCAGCATGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                  <  4:298096/135‑1 (MQ=255)
                                                                                                              atgcgcgGGCCGCCTCCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTg             <  8:216362/132‑1 (MQ=255)
                                                                                                                         gCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCg  <  4:282372/135‑1 (MQ=255)
                                                                                                                            |                                                                                                                                    
CTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCG  >  KT2440/3491202‑3491454

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: