Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,491,323 T→G 47.1% 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,3230TG47.1% ‑0.6 / 3.6 32intergenic (‑685/+36)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base G (3/9);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCA  >  KT2440/3491194‑3491448
                                                                                                                                 |                                                                                                                              
gCCCGCTCCTACATCGACCGGGGGTGGCGCAACGTTGTCGGCGCGGGTTAACCCGCTAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATCCACCGCCCCCGCCGCCACCGCATTCGCGGGCAGCCGGCCcgc                                                                                                                          <  4:114602/136‑1 (MQ=18)
             tCGCCCGGGGGTGTCGCCACGTTGTAGGAGCGGGTTTACCCGCTAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATCCACCGCCACCCCCCCCCCCGCATCCGCGGGCAGCCGCCCCGCGAACACCGgcgc                                                                                                             <  2:33647/135‑1 (MQ=21)
                 ccGGGGGCGGCGCCACGGGGGCGGCGCGGGTTACCCCGCGAATAGGCCGGACCGGCCGGCCACTCTCCCGGGCCAACCCCCGCCCCCCCCCCCCCCCCATCCCCGGGCCGCCGCCCCGCCAACACCGGCGCAGCCg                                                                                                        <  4:352091/136‑1 (MQ=255)
                  cGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGCACACCCGCGCAGCCgg                                                                                                       >  4:131180/1‑136 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCCCGCACACCGCCGCCGCCgg                                                                                                       >  5:161793/1‑135 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCCCGCACACCCGCGCAGCCggg                                                                                                      >  4:207967/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCgg                                                                                                       >  8:289971/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTACCCGGCACCCCGCCCGCGCACACCCGCGCCCCCgg                                                                                                       >  2:274261/1‑135 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCAGGCCACTATCACTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCCCGCAGCCACCCCGCCCACACACGCGCAccccgg                                                                                                      >  4:5929/1‑133 (MQ=18)
                          cccccAGGTTTTAGGCGCTGGTTAACCGGCAAATAGCCCGGACCTCCCGCCCAAAACCCGTGGACATTCCCCGCCACCCCCCCCCCCGCATGCGCGGGCAGCCGCCCCGCCAACACCGGCGCAGCCGGTTTCATc                                                                                                <  6:333738/132‑1 (MQ=255)
                                ggTTGTCGGAGCTGGTTTCCCCGCTAATTGGCGGGACCGGCCGGCCACTATCCCTGGTAATTCACCGTCACAGCCGCCCCCGCATGCGCGGGCAGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTc                                                                                          <  7:217806/134‑1 (MQ=255)
                                             ggTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGCACACCGGCGCCGCCGGGTTGGTCCCACTCACAGAcccacggc                                                                             >  7:346293/1‑129 (MQ=18)
                                                ttACCCGCTAATAGGCCGGACCGGCCGGCAACTATCCCGGGTCATACACCGCCACCGCCGCCACCCCATGCGCGGGCAGCCTCCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTg                                                                          <  3:309462/135‑1 (MQ=255)
                                                      gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCAGCCCGCGAACACCGGCGCCGCCGGGTTACTACAACTACACGCCCCTCGACCGGTACGc                                                                     >  8:122721/1‑134 (MQ=18)
                                                      gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCACCCAGCCCGCGAACACCCGCGCCCCCGGGTTCATCCAACACCACCACCCACGACCCGTACGCa                                                                    >  6:202428/1‑135 (MQ=18)
                                                      gcgAACAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGAACACCCGCGCCGCCCGGTTCCTCCAACTCCACCCCCCTCCACCGGTACGc                                                                     >  2:235709/1‑134 (MQ=18)
                                                           tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACCCCCGCGCACCCGGGGTCCTCCCACTCAACGAACAAGGCCCGGTGCGCAAAccgg                                                              >  4:308096/1‑133 (MQ=18)
                                                              gCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGtct                                                             <  7:289971/134‑1 (MQ=255)
                                                                              aCTATCCCTGGGCATTCCCCGGCCCAGCCCCCCCCGCATGCGCGGGCAGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCtt                                            <  8:366279/135‑1 (MQ=255)
                                                                                           aTTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCCGCGCAGCCCGTTTCATCCAACACAACCACCATCGACCTGTACCCAAACGTCTCTTCCATCCATTGCCTTGAATGGCCAATTGt                              >  4:230879/1‑136 (MQ=255)
                                                                                            ttCACCGGCACCGCCTCCACCGCATTCGGGGGAAGCGTGCCCGCGAACACCGGCTCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGt                              <  8:12993/135‑1 (MQ=255)
                                                                                                ccGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTg                         >  3:363011/1‑136 (MQ=255)
                                                                                                   gcacccccgccCCGGAATTCGCGGGCAGCCCGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGAt                       <  6:91153/129‑1 (MQ=255)
                                                                                                       cgccgccCCCCCATGCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  3:153366/135‑1 (MQ=255)
                                                                                                        gccgccCCCCCTTTCGCGGGAAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  6:103276/135‑1 (MQ=255)
                                                                                                        gcccccACCGAAGGCGCGGCCAGCGCCCCCGCCAACACCGGCGCAGCCGGTTTCATCCAACTCACCGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  2:128240/135‑1 (MQ=255)
                                                                                                         ccgccACCGCATGCGGGGGAAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATg                 <  1:245718/135‑1 (MQ=255)
                                                                                                          cccccccGCATTCGGGGGCCCCCCCCCCCCCACACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATg                 <  7:171141/130‑1 (MQ=255)
                                                                                                                    ttCGCGGGAAGCCGCCCGGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  3:76890/136‑1 (MQ=255)
                                                                                                                     gcgggggCAGCGTGCCCGGCAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  6:69535/131‑1 (MQ=255)
                                                                                                                     ccgcgGGCAGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  3:131180/134‑1 (MQ=255)
                                                                                                                       gcgGGCAGCCGCCCCGCGAACAGCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCGATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGca  <  3:230879/136‑1 (MQ=255)
                                                                                                                                 |                                                                                                                              
GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCA  >  KT2440/3491194‑3491448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: