Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,237,345 A→C 35.1% intergenic (+76/+25) dinG_3 → / ← CJ019_03036 3'‑5' exonuclease DinG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,237,3450AC35.1% 25.5 / 2.4 29intergenic (+76/+25)dinG_3/CJ019_030363'‑5' exonuclease DinG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/7);  new base C (3/8);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCGGCTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAGCACAGTGCCGAAACGCTTACCC  >  KT2440/3237251‑3237473
                                                                                              |                                                                                                                                
ggcggcggcTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAACCCACCCCCCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGt                                                                                         >  3:105454/1‑136 (MQ=255)
ggcggcggcTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACCCCGCCGCAAGGCGGCCCAAGCCCGTCAACCGGGCATTCCg                                                                                          >  1:52875/1‑135 (MQ=255)
        cTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCAGTCCGTCTGCCCgg                                                                                 >  3:218858/1‑136 (MQ=255)
        cTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCCCCCCGCCGCAAGGCGGCCCACGCACGTCCACCGGGCATTACGTCTTCCCgg                                                                                 >  2:362868/1‑135 (MQ=255)
        cTATGGCCTCGCGGGCTTACGCCTTTCCCCATCCGTCGCAACCCGGGGGGGGGAACAAACCACCCACCAAGTCCCCCTCATGCCCCCCCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCgg                                                                                 <  3:64614/136‑1 (MQ=255)
              cctcgCGGCCTTCCGCCTTTCCCGACCCCGCGCAACCCGGGGGGGCCACCAAACCACCCACCAAGTGCCCCTCATGCCACCCCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCg                                                                            <  5:132332/131‑1 (MQ=255)
                        tttCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATgg                                                                 >  3:314259/1‑136 (MQ=255)
                        tttCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATg                                                                  >  1:305044/1‑135 (MQ=255)
                           cGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCCCCCCGCCCCAAGGCGGCCCAATCACGTCAACCGGGCATTCCGTCTGCCCGGCGCCGGAGCAGGAttg                                                                 >  7:171996/1‑131 (MQ=255)
                           cGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATggg                                                                >  3:114804/1‑134 (MQ=255)
                           cGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTACCCTCATGGCACACCGCCGCAAGGCGGCCCCAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGCGCCGGATggt                                                                >  4:5378/1‑133 (MQ=255)
                            gCCTTTGACTATAAGTCGCAATCAGGGGGGGGCAACAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGAAAGGCGGCCCAAGCACGTCAACAGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCa                                                              <  4:157637/135‑1 (MQ=255)
                                   aCGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATg                                                       >  6:234181/1‑135 (MQ=255)
                                                      gggggggAGCAAACCACCCACAAAGTCCCCCTCATGGCACCCCGCCGCACGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAgcg                                   <  4:70717/136‑1 (MQ=255)
                                                      ggggggcAACAACCCACCCACCAAGCTCCCCTCCGGCCACCCCGCCCCATGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAgcg                                   <  3:208795/136‑1 (MQ=255)
                                                       gggggtACCAAACCACCCACCAAGTCCCCCTCATGCCACCCCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAgcg                                   <  8:171996/135‑1 (MQ=255)
                                                       ggggggATCAAACCACCCACAAAGTGCCCCACATGCCACCCCCCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAgcg                                   <  2:8458/135‑1 (MQ=255)
                                                       ggggggAACAAACCACCAACCAAGGACCCCTCATGGCACCCCGCCGCACGGCGGCCCAAGCACGTCACCCGGGCATTCCGGCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAgcg                                   <  8:277917/135‑1 (MQ=255)
                                                        gggggAACAAACCCCCCACCAAGTGCCCCTCTTCCACCACCTCCCCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGAGCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGcc                                 <  3:89502/136‑1 (MQ=255)
                                                         gcgcaccAACCCACCCACCAAGTACCCCTCATGGCACACCGCCCCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGcc                                 <  8:202090/129‑1 (MQ=255)
                                                            gATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGcctgt                              >  5:332167/1‑132 (MQ=255)
                                                            caacAAACCCCCCACCACGTACCCCGCCGGCCCCCCCCCCGCACGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAgg                              <  7:5148/132‑1 (MQ=255)
                                                                aacccacccacCAAGTACCCCTCATGCCACACCCCCGCACGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCa                         <  4:218858/133‑1 (MQ=255)
                                                                aaaCCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTc                           >  3:17185/1‑134 (MQ=255)
                                                                 aaCCACCCACAAAGGTCCCCTCATGCCACACCGCCGCACGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCa                         <  8:360654/135‑1 (MQ=255)
                                                                  accaccCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAg                        >  6:259129/1‑135 (MQ=255)
                                                                   ccacccacCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGACCAGCAACAGCGCCAGGCTCCAGc                       >  8:116321/1‑135 (MQ=255)
                                                                                       aTGGAACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCAGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAGCACAGTGCCGAAACGCTTAccc  <  4:264848/136‑1 (MQ=255)
                                                                                        tgccacaCCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAGCACAGTGCCGAAACGCTTAc    <  6:135462/130‑1 (MQ=255)
                                                                                        tGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAGCACAGTGCCGAAACGCTTAccc  <  2:69628/135‑1 (MQ=255)
                                                                                              |                                                                                                                                
GGCGGCGGCTATTGTCTAGCGGTCTTTCGCCTTTGACGATCAGTCGCAATCAGGGGGGGGGATCAAACCACCCACCAAGTTCCCCTCATGGCACACCGCCGCAAGGCGGCCCAAGCACGTCAACCGGGCATTCCGTCTGCCCGGGGCCGGAGCAGGATGGGCAGGTAATGGCGCAGGAACAGCAACAGCGCCAGGCTCCAGCACAGTGCCGAAACGCTTACCC  >  KT2440/3237251‑3237473

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: