Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,802,244 T→A 14.8% intergenic (‑706/‑25) CJ019_04509 ← / → dsbD_2 hypothetical protein/Thiol:disulfide interchange protein DsbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,802,2440TA14.8% 59.3 / 5.8 27intergenic (‑706/‑25)CJ019_04509/dsbD_2hypothetical protein/Thiol:disulfide interchange protein DsbD
Reads supporting (aligned to +/- strand):  ref base T (15/8);  new base A (2/2);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

GCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGAAACCAGACTTCCTGCCGGTTGACCAAGCCTTCGT  >  KT2440/4802121‑4802372
                                                                                                                           |                                                                                                                                
gccTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTg                                                                                                                       <  5:175697/135‑1 (MQ=255)
                  aacaacaaAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATgc                                                                                                      <  6:65210/134‑1 (MQ=255)
                                       ccAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTcac                                                                                <  4:258807/135‑1 (MQ=255)
                                         aGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACt                                                                             >  8:45433/1‑136 (MQ=255)
                                           aaCAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTc                                                                            >  4:25355/1‑135 (MQ=255)
                                             cAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCAATCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  <  8:173679/127‑1 (MQ=255)
                                             cAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCAATCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  <  8:169366/127‑1 (MQ=255)
                                             cAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCAATCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  >  7:173679/1‑127 (MQ=255)
                                             cAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCAATCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  >  7:169366/1‑127 (MQ=255)
                                             cAAACTCTTACTACTCCTGACCAACCCCCAACACCCGCCTTAAACTTCATTAAATCAAACCCAGCCAACATGGCCACATTCACAGATCTGCCAAGGCGAACACATTCTCGTTCACCTGCTCATCCAAACACTCCt                                                                          >  6:101768/1‑135 (MQ=255)
                                                   ctTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTggc                                                                    >  1:94838/1‑135 (MQ=255)
                                                               gACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGc                                                         >  7:199905/1‑134 (MQ=255)
                                                                   cacccCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGcgcg                                                                                                   >  6:54522/1‑88 (MQ=255)
                                                                   cacccCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGcgcg                                                                                                   <  5:54522/88‑1 (MQ=255)
                                                                        ccAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGCCCGGCCCGTTGCAAGCCAACCCg                                              >  7:95558/1‑136 (MQ=255)
                                                                            cacCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAAccc                                               <  7:105710/131‑1 (MQ=255)
                                                                            cacCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAAccc                                               >  8:105710/1‑131 (MQ=255)
                                                                             acCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  <  2:162093/95‑1 (MQ=255)
                                                                             acCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTc                                                                                  >  1:162093/1‑95 (MQ=255)
                                                                             acCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCg                                          >  6:216176/1‑135 (MQ=255)
                                                                                               ttAATCAACCCCAGCCAGCATGGCCCCATTCCCAGAACTGCCAAGGAGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCAACAGTCCGGCCCGTTGCAAGCCAACCCGTTAGATGTTAAACCAGCCTTCCt                       >  7:257976/1‑136 (MQ=255)
                                                                                                      aCCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGAAACCAGACTTCCTGCCGGt                 >  5:77574/1‑135 (MQ=255)
                                                                                                         ccagcctgCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCAGACTCCTACTGGGCGGCCCGTTGCTAGCCATCCCGTTCGATGTGAAACGAGACTTCCTGCCGGTTGa              <  5:101768/135‑1 (MQ=255)
                                                                                                              cagcATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGAAACCAGACTTCCTGCCGGTTGACCAAg         <  3:25355/135‑1 (MQ=255)
                                                                                                               agcATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGAAACCAGACTTCCTGCCGGTTGACCAAGcc       >  8:34374/1‑136 (MQ=255)
                                                                                                                    ggCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCACCCCGTTCGATGTGAAACCAGACTTCCTGCCGGTTGACCAAGCCTTCGt  >  8:277863/1‑136 (MQ=255)
                                                                                                                    ggCCCCATTCCCAGATCAGCCAAGCCGAAGACATGCGCGTTCTCCTGCTCTTCCACACACTCCTACTGGCCTCCCCGATGCACGCCAACCCGTACGATGTGAAACCAGACTTCCTGCCGGTTGACCAAGccaacg   >  4:229175/1‑131 (MQ=255)
                                                                                                                           |                                                                                                                                
GCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTCATTTAATCAACCCCAGCCAGCATGGCCCCATTCCCAGATCTGCCAAGGCGAAGACATGCGCGTTCTCCTGCTCTTCCTCACACTCCTACTGGCCGGCCCGTTGCAAGCCAACCCGTTCGATGTGAAACCAGACTTCCTGCCGGTTGACCAAGCCTTCGT  >  KT2440/4802121‑4802372

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