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

GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCA  >  KT2440/3491194‑3491448
                                                                                                                                 |                                                                                                                              
gCCCGCTCCTACATCGACCGGGGGTGGCGCAACGTTGTCGGCGCGGGTTAACCCGCTAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATCCACCGCCCCCGCCGCCACCGCATTCGCGGGCAGCCGGCCcgc                                                                                                                          <  4:114602/136‑1 (MQ=18)
             tCGCCCGGGGGTGTCGCCACGTTGTAGGAGCGGGTTTACCCGCTAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATCCACCGCCACCCCCCCCCCCGCATCCGCGGGCAGCCGCCCCGCGA‑ACACCGgcgc                                                                                                             <  2:33647/135‑1 (MQ=21)
                 ccGGGGGCGGCGCCACGGGGGCGGCGCGGGTTACCCCGCGAATAGGCCGGACCGGCCGGCCACTCTCCCGGGCCAACCCCCGCCCCCCCCCCCCCCCCATCCCCGGGCCGCCGCCCCGCCA‑ACACCGGCGCAGCCg                                                                                                        <  4:352091/136‑1 (MQ=255)
                  cGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGC‑ACACCCGCGCAGCCgg                                                                                                       >  4:131180/1‑136 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCCCGC‑ACACCGCCGCCGCCgg                                                                                                       >  5:161793/1‑135 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCCCGC‑ACACCCGCGCAGCCggg                                                                                                      >  4:207967/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCgg                                                                                                       >  8:289971/1‑135 (MQ=255)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTACCCGGCACCCCGCCCGCGC‑ACACCCGCGCCCCCgg                                                                                                       >  2:274261/1‑135 (MQ=21)
                   ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCAGGCCACTATCACTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCCCGCAGCCACCCCGCCC‑ACACACGCGCAccccgg                                                                                                      >  4:5929/1‑133 (MQ=18)
                          cccccAGGTTTTAGGCGCTGGTTAACCGGCAAATAGCCCGGACCTCCCGCCCAAAACCCGTGGACATTCCCCGCCACCCCCCCCCCCGCATGCGCGGGCAGCCGCCCCGCCA‑ACACCGGCGCAGCCGGTTTCATc                                                                                                <  6:333738/132‑1 (MQ=255)
                                ggTTGTCGGAGCTGGTTTCCCCGCTAATTGGCGGGACCGGCCGGCCACTATCCCTGGTAATTCACCGTCACAGCCGCCCCCGCATGCGCGGGCAGCCGGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTc                                                                                          <  7:217806/134‑1 (MQ=255)
                                             ggTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCGCGC‑ACACCGGCGCCGCCGGGTTGGTCCCACTCACAGAcccacggc                                                                             >  7:346293/1‑129 (MQ=18)
                                                ttACCCGCTAATAGGCCGGACCGGCCGGCAACTATCCCGGGTCATACACCGCCACCGCCGCCACCCCATGCGCGGGCAGCCTCCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTg                                                                          <  3:309462/135‑1 (MQ=255)
                                                      gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCAGCCCGCGA‑ACACCGGCGCCGCCGGGTTACTACAACTACACGCCCCTCGACCGGTACGc                                                                     >  8:122721/1‑134 (MQ=18)
                                                      gcgAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCACCCAGCCCGCGA‑ACACCCGCGCCCCCGGGTTCATCCAACACCACCACCCACGACCCGTACGCa                                                                    >  6:202428/1‑135 (MQ=18)
                                                      gcgAACAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGA‑ACACCCGCGCCGCCCGGTTCCTCCAACTCCACCCCCCTCCACCGGTACGc                                                                     >  2:235709/1‑134 (MQ=18)
                                                           tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGC‑ACCCCCGCGCACCCGGGGTCCTCCCACTCAACGAACAAGGCCCGGTGCGCAAAccgg                                                              >  4:308096/1‑133 (MQ=18)
                                                              gCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGtct                                                             <  7:289971/134‑1 (MQ=255)
                                                                              aCTATCCCTGGGCATTCCCCGGCCCAGCCCCCCCCGCATGCGCGGGCAGCCGGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCtt                                            <  8:366279/135‑1 (MQ=255)
                                                                                           aTTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCCGCGCAGCCCGTTTCATCCAACACAACCACCATCGACCTGTACCCAAACGTCTCTTCCATCCATTGCCTTGAATGGCCAATTGt                              >  4:230879/1‑136 (MQ=255)
                                                                                            ttCACCGGCACCGCCTCCACCGCATTCGGGGGAAGCGTGCCCGCGA‑ACACCGGCTCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGt                              <  8:12993/135‑1 (MQ=255)
                                                                                                ccGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTg                         >  3:363011/1‑136 (MQ=255)
                                                                                                   gcacccccgccCCGGAATTCGCGGGCAGCCCGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGAt                       <  6:91153/129‑1 (MQ=255)
                                                                                                       cgccgccCCCCCATGCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  3:153366/135‑1 (MQ=255)
                                                                                                        gccgccCCCCCTTTCGCGGGAAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  6:103276/135‑1 (MQ=255)
                                                                                                        gcccccACCGAAGGCGCGGCCAGCGCCCCCGCCA‑ACACCGGCGCAGCCGGTTTCATCCAACTCACCGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAAt                  <  2:128240/135‑1 (MQ=255)
                                                                                                         ccgccACCGCATGCGGGGGAAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATg                 <  1:245718/135‑1 (MQ=255)
                                                                                                          cccccccGCATTCGGGGGCCCCCCCCCCCCCACACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATg                 <  7:171141/130‑1 (MQ=255)
                                                                                                                    ttCGCGGGAAGCCGCCCGGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  3:76890/136‑1 (MQ=255)
                                                                                                                     gcgggggCAGCGTGCCCGGCA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  6:69535/131‑1 (MQ=255)
                                                                                                                     ccgcgGGCAGCCGGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGaa     <  3:131180/134‑1 (MQ=255)
                                                                                                                       gcgGGCAGCCGCCCCGCGA‑ACAGCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCGATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGca  <  3:230879/136‑1 (MQ=255)
                                                                                                                                 |                                                                                                                              
GCCCGCTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGA‑ACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCA  >  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: ‑