Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,605,842 G→A 10.2% G575D (GGC→GAC)  CJ019_02463 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,605,8420GA10.2% 111.6 / 2.2 39G575D (GGC→GAC) CJ019_02463hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (2/2);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAGCCACCGTGTACAAGATGCACGCCACCA  >  KT2440/2605717‑2605976
                                                                                                                             |                                                                                                                                      
ggACCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCaa                                                                                                                                 <  5:83214/133‑1 (MQ=255)
ggACCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAAt                                                                                                                               <  6:159707/135‑1 (MQ=255)
ggACCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAAt                                                                                                                               >  3:240302/1‑135 (MQ=255)
ggACCTCCGCAAGATGAACTTCAAGGGCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAAt                                                                                                                               <  8:23424/135‑1 (MQ=255)
  aCCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGACATTCAACTGt                                                                                                                             >  1:141164/1‑135 (MQ=255)
  aCCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGACATTCAAATGt                                                                                                                             >  7:264725/1‑135 (MQ=255)
         cAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCg                                                                                                                      >  8:188056/1‑135 (MQ=255)
          aaGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGACATTCAAATGTTCCTCCgg                                                                                                                     <  2:141164/135‑1 (MQ=255)
          aaGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGACATTCAAATGTTCCTCCgg                                                                                                                     <  8:264725/135‑1 (MQ=255)
           aGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCggg                                                                                                                    <  3:240611/135‑1 (MQ=255)
               gAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGCTGGCATTCAAATGTTCCTCCGGGGCCa                                                                                                                >  3:261307/1‑135 (MQ=255)
               gAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCa                                                                                                                <  6:204867/135‑1 (MQ=255)
                    tCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCt                                                                                                           >  8:176078/1‑135 (MQ=255)
                    tCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCt                                                                                                           >  2:247636/1‑135 (MQ=255)
                               gAAAAGGGACAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGAGCGGTTCTCCCCCCGCAACAGTCCAGCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGcccc                                                                                                >  2:264937/1‑135 (MQ=255)
                                aaaaGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCa                                                                                               >  6:96277/1‑135 (MQ=255)
                                aaaaGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCa                                                                                               >  3:273436/1‑135 (MQ=255)
                                       ccAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCg                                                                                        >  6:100836/1‑135 (MQ=255)
                                          gAAGTTTGGCAATGCGGCGATGCACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGCGGATCGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCt                                                                                     >  4:273310/1‑135 (MQ=255)
                                                   cAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGAccg                                                                            >  8:162/1‑134 (MQ=255)
                                                   cAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCt                                                                            >  1:25957/1‑135 (MQ=255)
                                                   cAATGCGGCGATGGACGCGGTGATTAAAAATGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCCGGCGTCTGTGGGGccg                                                                            >  3:200335/1‑134 (MQ=255)
                                                           cgATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGCCCTGTCCGGcc                                                                    >  7:45392/1‑135 (MQ=255)
                                                                        gATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGTCCATTCCTACATGAGCCCCACCGTGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAAc                                                       <  1:244131/135‑1 (MQ=255)
                                                                        gATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAAc                                                       <  1:38599/135‑1 (MQ=255)
                                                                                    tAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCAt                                             <  1:59338/133‑1 (MQ=255)
                                                                                    tAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCAt                                             <  1:188050/133‑1 (MQ=255)
                                                                                     aaGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCCCCGGGCGTCTGTGGGGCCTGTCCGGCCTGAATGAAGAAACCGTAGTGCATgtg                                          >  3:120826/1‑133 (MQ=255)
                                                                                      aGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGGCCGGTCCGGCCGGAATGAAGAAACCTTCGTGCATGAGc                                         >  3:108494/1‑135 (MQ=255)
                                                                                              gTTCTCACCCCGCACCAGTCCACCGGGGAGGGCATTCAAATGTTCCTCCGGGGCCATTCATACATGAGCCCCCCCGGGCGTCGGTCGGACCGGGCCGGCCGGAATGAAGCAACCTACGTGCAGGGTCTGGTCCAc                                 >  6:122660/1‑135 (MQ=255)
                                                                                                  tcACCCCGCCCCAGCCCACAGGGGAAGGCCATCAAATGTGCCTCCGGGGCCATTCCTACCTGAGCCCCACCGGGCTTCTGTCGGACATGTCCGGCCGGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAg                             <  7:286462/135‑1 (MQ=255)
                                                                                                     ccccGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAGcc                           <  2:38268/134‑1 (MQ=255)
                                                                                                     ccccGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAGc                            <  6:59992/133‑1 (MQ=255)
                                                                                                            ccAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTGGGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCCGCCAcccggt                    >  7:105050/1‑130 (MQ=255)
                                                                                                            ccAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCCGCCTGAATGAAGAAACCATCGTGCATGAGCTGGTCCCCGCAGCCACCctgt                    >  7:45286/1‑134 (MQ=255)
                                                                                                            ccAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGCCCATTCCTACATGAGCCCCACCGGCCGTCAGTCGGCCCTGTCCGGCCTGAAAGAATAAAACCACCTGCAAGAACAGGTCCCCGCAGCCAccctct                    >  2:159629/1‑130 (MQ=255)
                                                                                                              aGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCa                                                                                                                <  1:156899/40‑1 (MQ=255)
                                                                                                              aGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCa                                                                                                                >  2:156899/1‑40 (MQ=255)
                                                                                                                    ccGGGTATGGCATTCAAATGTTCCTCGGGGGCCATTCCTACATGACCCGCACCGGGGGTCTGACGGACCTGTCGGGCCTGATTGAAGATACCTTCGTGCATGAGCTGGACCACCCTGCCACCGTGTACAAGATGc           <  6:302209/135‑1 (MQ=255)
                                                                                                                             gCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAGCCACCGTGTACAAGATGCACGccacca  >  4:128576/1‑135 (MQ=255)
                                                                                                                             gCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCCGACCTGGCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACCCAGCCACCGTGTACAAGATGCACCccccca  >  7:253018/1‑135 (MQ=255)
                                                                                                                             |                                                                                                                                      
GGACCTCCGCAAGATGAACTTCAAGGTCATCGAAAAGGGCCAGAAGTTTGGCAATGCGGCGATGGACGCGGTGATTAAAACTGCTAAGGGCGCGGTTCTCACCCCGCACCAGTCCACCGGGGATGGCATTCAAATGTTCCTCCGGGGCCATTCCTACATGAGCCCCACCGGGCGTCTGTCGGACCTGTCCGGCCTGAATGAAGAAACCTTCGTGCATGAGCTGGTCCACGCAGCCACCGTGTACAAGATGCACGCCACCA  >  KT2440/2605717‑2605976

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: ‑