Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,851,109 T→A 9.1% V35E (GTG→GAG)  CJ019_04549 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,851,1090TA 9.1% 127.5 / 4.5 44V35E (GTG→GAG) CJ019_04549hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base A (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAAC  >  KT2440/4850979‑4851237
                                                                                                                                  |                                                                                                                                
gATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGc                                                                                                                                <  6:84821/133‑1 (MQ=255)
  ttAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCt                                                                                                                            >  2:210340/1‑135 (MQ=255)
       ccAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCtttttt                                                                                                                       <  5:65534/135‑1 (MQ=255)
               aaGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAACTGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  2:257350/135‑1 (MQ=255)
               aaGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  8:42889/135‑1 (MQ=255)
                aGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgcg                                                                                                              <  1:13017/135‑1 (MQ=255)
                        cAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATctct                                                                                                      >  4:71165/1‑135 (MQ=255)
                         aaaTGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTa                                                                                                     <  7:63898/135‑1 (MQ=255)
                               aCCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCGCCTTTTTTACATCGGCGCAATCTCTACAggg                                                                                                >  2:172020/1‑134 (MQ=255)
                                ccGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  1:183735/118‑1 (MQ=255)
                                ccGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               >  2:183735/1‑118 (MQ=255)
                                ccGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGt                                                                                              >  7:249858/1‑135 (MQ=255)
                                ccGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTCCTCCGTTTTTACATCGGCGCAATCTCTACAGGGGt                                                                                              >  1:294325/1‑135 (MQ=255)
                                      cgccCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTActgct                                                                                        >  7:93317/1‑135 (MQ=255)
                                         ccAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTActgctg                                                                                       >  8:156294/1‑133 (MQ=255)
                                                             cACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTgc                                                                 <  8:72284/135‑1 (MQ=255)
                                                                tgaatgaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGgtgt                                                                                                                                  >  6:305931/1‑67 (MQ=255)
                                                                tgaatgaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGgtgt                                                                                                                                  <  5:305931/67‑1 (MQ=255)
                                                                  aatgaaGCAGAAGGGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  5:38695/84‑1 (MQ=255)
                                                                  aatgaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGACATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               >  6:38695/1‑84 (MQ=255)
                                                                   atgaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGcgctcg                                                            <  3:71165/134‑1 (MQ=255)
                                                                     gaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               >  1:285147/1‑81 (MQ=255)
                                                                     gaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  2:285147/81‑1 (MQ=255)
                                                                     gaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCgctcgctc                                                         >  7:78790/1‑135 (MQ=255)
                                                                     gaaGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCgctcgctc                                                         >  3:258384/1‑135 (MQ=255)
                                                                          aGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGAc                                                     >  4:105340/1‑134 (MQ=255)
                                                                          aGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGAc                                                     >  2:248180/1‑134 (MQ=255)
                                                                               cGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgcg                                                                                                              >  5:236710/1‑72 (MQ=255)
                                                                               cGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgcg                                                                                                              <  6:236710/72‑1 (MQ=255)
                                                                                    gCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAggg                                                                                                <  6:41617/81‑1 (MQ=255)
                                                                                    gCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAggg                                                                                                >  5:41617/1‑81 (MQ=255)
                                                                                    gCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTg                                          >  6:253293/1‑135 (MQ=255)
                                                                                     ccGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGc                                         <  4:258384/135‑1 (MQ=255)
                                                                                         agaaATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTc                                     >  2:47035/1‑135 (MQ=255)
                                                                                                   aaCCTTCGAGCGTCTTCGCCATGGACGGGTGAGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTNCCGTCNCANTAAAg                            >  8:77679/1‑134 (MQ=255)
                                                                                                   aaCCTTCGAGCGTCTTCGCCATGGACGGGTGAGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAg                            >  8:72529/1‑134 (MQ=255)
                                                                                                           aGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAAc                   >  7:47019/1‑135 (MQ=255)
                                                                                                           aGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAAc                   <  7:156294/135‑1 (MQ=255)
                                                                                                             cGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAAccc                 <  1:172020/135‑1 (MQ=255)
                                                                                                              gTCTTCGCCATGGACGGGTGAGCTCCTTTTTTACATNGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCg                <  7:77679/135‑1 (MQ=255)
                                                                                                              gTCTTCGCCATGGACGGGTGAGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCg                <  7:72529/135‑1 (MQ=255)
                                                                                                                  tCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               <  4:298195/36‑1 (MQ=255)
                                                                                                                  tCGCCATGGACGGGTGTGCTCCTTTTTTACATCGgc                                                                                                               >  3:298195/1‑36 (MQ=255)
                                                                                                                     ccATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGTTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCaa         <  8:179728/135‑1 (MQ=255)
                                                                                                                     ccATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCaa         <  1:210340/135‑1 (MQ=255)
                                                                                                                            cGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGaacaac  >  4:9961/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                                
GATTAATCCAGCCACAAGGAGTAACAAATGCACCGCGCCGCCCAGGATCACACCAGAACAGCACTGAATGAAGCAGAAGCGTTCGCCGAAGAAATGCGCAACCTTCGAGCGTCTTCGCCATGGACGGGTGTGCTCCTTTTTTACATCGGCGCAATCTCTACAGGGGTACTGCTGTTAGCAGGCATGCTGATTTTGCGCTCGCTCTGACCCAACTCCCTGCCGTCACACTAAAGCAAGTAAACCCGCACTCAAGAACAAC  >  KT2440/4850979‑4851237

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