Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,692,266 G→T 100% E320* (GAA→TAA)  cheB_1 → Protein‑glutamatemethylesterase/protein‑glutamine glutaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,692,2660GT77.0% 55.7 / 18.7 35E320* (GAA→TAA) cheB_1Protein‑glutamatemethylesterase/protein‑glutamine glutaminase
Reads supporting (aligned to +/- strand):  ref base G (7/1);  new base T (13/14);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGCCGCGACGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGATTTACCGATCGAAGG  >  KT2440/1692139‑1692396
                                                                                                                               |                                                                                                                                  
gggCCGCGACGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGCCCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCACGGCTGCTGCGGCTATTGATGCCGCCGGGTAAATccgc                                                                                                                            >  8:87744/1‑136 (MQ=255)
      cgACGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCt                                                                                                                       <  5:136517/135‑1 (MQ=255)
         cGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTggc                                                                                                                    >  1:6578/1‑134 (MQ=255)
         cGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGa                                                                                                                    >  1:11311/1‑135 (MQ=255)
                cAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCAtt                                                                                                             <  5:101639/135‑1 (MQ=255)
                           gCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCtt                                                                                                  >  3:67698/1‑135 (MQ=255)
                                       gcagGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCACGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAAttt                                                                                        >  6:114838/1‑133 (MQ=255)
                                         agGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGCCTGAAtttttg                                                                                     >  1:77785/1‑133 (MQ=255)
                                             cTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCaaa                                                                                <  5:255393/135‑1 (MQ=255)
                                                 cTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATg                                                                             <  1:251710/134‑1 (MQ=255)
                                                    aCCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGAc                                                                           <  7:44576/133‑1 (MQ=255)
                                                     ccATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGAc                                                                           <  3:260834/132‑1 (MQ=255)
                                                     ccATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGAc                                                                           <  3:49617/132‑1 (MQ=255)
                                                       aTTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAAt                                                                      <  1:5064/135‑1 (MQ=255)
                                                          gCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTa                                                                   >  1:206673/1‑135 (MQ=255)
                                                                     ggCCCGCCGCGCGGTGTACGGCTTGCCCATGGCTGCTGCGGCTATTGTGGCCGCGGGTTAAAGCCCCCCGCTGGCTCGCATTGCCGGGCGCTTGAGTGAATTTTTCGCAAAATGAGGATTTGTATTCAAGCAcgc                                                        <  6:154173/135‑1 (MQ=255)
                                                                        cagcagCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATTAAGAAATGTATTGAATCAcgccgc                                                     >  6:249199/1‑135 (MQ=255)
                                                                        cagcagCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCGTGCCTGAATTTTTCGCAAAATGACGAAATGTATTGAATAAcgccg                                                      >  4:256188/1‑134 (MQ=255)
                                                                          gcagcGCGGTGTACGGCATGCCCAAGGCTGATGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCTTTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCa                                                   <  4:111348/135‑1 (MQ=255)
                                                                             gcgcGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGcc                                               >  8:262958/1‑136 (MQ=255)
                                                                                    gtagtgCATGCAAAAGGGTGCTGGGGCTATTGATGCGGCCGTGTAAATCCGCCGGCTGTAACGCTTTGCCGGGCGCTTGACTGAATTTTTCTAAAATTGACGAATTTTATTGAATCACGCCGCCAGGCCGGTAGt                                         <  4:115348/130‑1 (MQ=255)
                                                                                    gtACGGCATGCCCAAGGCTGCTGCGGCTATTGAGGCCGCCGTGTAAATCCGCCCGCTGGACCGCATTGCCGGGCGCTTGCCTGAATTTTTCGCAAAATCCCGAAATGTATTGAATCACGCCGCCAGGCCGGTAGt                                         <  3:173188/135‑1 (MQ=255)
                                                                                                cAAGGCTGCTGTGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAAc                             >  7:21094/1‑135 (MQ=255)
                                                                                                cAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATTACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGaa                              >  2:10015/1‑134 (MQ=255)
                                                                                                cAAGGCTGCTGCGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACGGAATTTTTGGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGaa                              >  5:122898/1‑134 (MQ=255)
                                                                                                cAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAAGTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTgca                              >  5:78159/1‑132 (MQ=255)
                                                                                                cAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCAATGCCGGGCGCGTGACTGAATTTTTCGCAAAATGACGAAAAGTATTGAATCACGCCGCCAGGCCGGGTGTGATCGGGTGaa                              >  3:115348/1‑134 (MQ=255)
                                                                                                     ctgctgCGGCTATTGATGCCGCCGTGTAAATACGCCCGCTGGAACGCATTGCCGGGCGCTTGCCTGAATTTTTCGCAAAATTCCGAAATGTATTGAATCACGCCGCCAGGGCGGTAGTGATCGGGTGTTCAatgat                       >  7:58570/1‑136 (MQ=255)
                                                                                                           cGGCTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATAACGCCGCCAGGCCGGTAGTGATCGGGTGTAAAATGatgaatga                 >  7:122739/1‑136 (MQ=255)
                                                                                                           cGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTAGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGatgaatg                  >  3:63833/1‑135 (MQ=255)
                                                                                                              cTATTGATGCCGCCGTGTAACACCGCCCGCGGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGAtt               <  2:249514/135‑1 (MQ=255)
                                                                                                              cTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGAtt               <  2:206673/135‑1 (MQ=255)
                                                                                                              cTATTGATGCCGCCGTGTAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGAtt               <  1:10015/135‑1 (MQ=255)
                                                                                                                         gccgTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGATTTACCGATCGaa    >  2:32218/1‑135 (MQ=255)
                                                                                                                           cgTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTTACTGAATTTTTAGCAAAATGACGAAATGTATTGAATCACCCCCCCCGGCCGGTAGTGATCGGGTGCACAATGATGAATGAATTACCGATCGAAgg  >  6:194969/1‑135 (MQ=255)
                                                                                                                               |                                                                                                                                  
GGGCCGCGACGGTGCCCAGGGCCTGAAGCTGATGCGTCAGCAGGGCTTTCTGACCATTGCCCAGGACCAGGCCAGCAGCGCGGTGTACGGCATGCCCAAGGCTGCTGCGGCTATTGATGCCGCCGTGGAAATCCGCCCGCTGGAACGCATTGCCGGGCGCTTGACTGAATTTTTCGCAAAATGACGAAATGTATTGAATCACGCCGCCAGGCCGGTAGTGATCGGGTGAACAATGATGAATGATTTACCGATCGAAGG  >  KT2440/1692139‑1692396

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: