Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,949,404 T→C 12.9% F238F (TTT→TTC)  ilvA → L‑threonine dehydratase, biosynthetic; also known as threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,949,4040TC12.9% 68.2 / 4.8 31F238F (TTT→TTC) ilvAL‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base T (16/11);  new base C (2/2);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-03

CGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAA  >  CP009273/3949257‑3949527
                                                                                                                                                   |                                                                                                                           
cggcggTCTGGCTGCTGGCGTGCGGGTGCTGATCAAACAACTGATGCCGCAAATCAACGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCTGGTTGATCTGCCGCGCGTAGGGCTATTTg                                                                                                                            <  1:206149/149‑1 (MQ=255)
     gTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGaa                                                                                                                       <  1:109488/149‑1 (MQ=255)
      tgtGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAg                                                                                                                      <  2:138710/147‑1 (MQ=255)
            tgctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGAAAATCAAAGTGATCGCCATAGAAGCTGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAg                                                                                                                >  1:210963/1‑149 (MQ=255)
              ctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGTATGCAGCGCTGGTTGTGGGTTTTCTGGTTGGTCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCg                                                                                                              <  1:95516/149‑1 (MQ=255)
                 gcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCATGAAAGCAGAGCTGAATGCGGGTCATCCGGTTGATCTACCGCGCGTAGGGCTATTTGCTAAAGGGGTCGCGGtg                                                                                                           >  1:129621/1‑148 (MQ=255)
                 gcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTACGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAATGCGTAGCGGTa                                                                                                           >  1:327347/1‑149 (MQ=255)
                     ggcgGTGCTGATCAAACAACTGATGCGGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCCGTTTGATGAAGAAGTAGCGGTgaaag                                                                                                       >  2:302512/1‑148 (MQ=255)
                               aTCAAACAACTGATGCCGCACATCAAAGTGATCACCGTAAAAGCGAAAAACTCAGCCTGCCTAAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGGAAAACGCATCGGTGa                                                                                             >  1:681/1‑149 (MQ=255)
                                  aaaCAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGa                                                                                          >  1:242079/1‑149 (MQ=255)
                                  aaaCAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGAAATCCTGTAGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGa                                                                                          >  1:241951/1‑149 (MQ=255)
                                     caacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGTGCTGGATGCGGGTCATCCGGTTGATCTGCCGTGCGTAGGGCTATTTGCTGAAGGTGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                                       <  1:196604/149‑1 (MQ=255)
                                     caacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCTCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                                       <  1:195864/149‑1 (MQ=255)
                                            tGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGt                                                                                <  2:134125/149‑1 (MQ=255)
                                                        aaGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAg                                                                    <  2:288192/149‑1 (MQ=255)
                                                                gCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcga                                                            >  2:197459/1‑149 (MQ=255)
                                                                       aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTCGCTGAAGGCGTAGCGGTAAAACGCATCGGTGGCGAAACCTTCCGTTTATGCCAGGGGTATCTCGACGAcctc                                                     >  2:28098/1‑149 (MQ=255)
                                                                       aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTCGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAaatc                                                     >  2:28124/1‑149 (MQ=255)
                                                                               gACACCGCCTGCCTGAAACCAGCGCTGGATGCGGGTCATCCGGTTGAACTGCCGAGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTACGACGCATCGGGGACGAAACCTTCCGTTTATGCCAGGAGAACTACCACCACATCAGCACCGt                                             >  2:93680/1‑149 (MQ=255)
                                                                                aCTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTc                                            >  1:123231/1‑149 (MQ=255)
                                                                                aCTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCATAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTc                                            >  1:113132/1‑149 (MQ=255)
                                                                                aCTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTAATCTACCGCGCGAAGGGCTATTTGCTGAAGGAGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTc                                            >  1:81575/1‑149 (MQ=255)
                                                                                  tCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGATAAAGGCGTAACGGTAAAACGCATCGGTGCCGAAACCTTCCGTTTATGCCAGGAGTATCTCGCCGACATCATCACCGTCGa                                          >  1:79067/1‑149 (MQ=255)
                                                                                       ctgcctgAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcg                                     >  1:269011/1‑149 (MQ=255)
                                                                                         gcctgATGGCAGCGCTGGCTGCGGGTCATCCGTTTGATCTGCCGCGCGTAGGGCTATTTGCTGTAGGCGTAGCGGTACAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgat                                   <  2:83954/149‑1 (MQ=255)
                                                                                             gAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACGGTCGATAgcgatgcgg                               >  2:320117/1‑148 (MQ=255)
                                                                                             gAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcga                               <  2:26596/149‑1 (MQ=255)
                                                                                               aaGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCAATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCCCGACGACATCATCACCGTCGCAAGCGATGCGATc                             >  2:250132/1‑149 (MQ=255)
                                                                                                                         ctttaTCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGa   <  2:78731/145‑1 (MQ=255)
                                                                                                                          ttGATCTGCCGCGCGTAGGGCTATTCGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGTGTATCTCGATGACATCATCACTGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGaa  <  1:28124/149‑1 (MQ=255)
                                                                                                                          ttGATCTGCCGCGCGTAGGGCTATTCGCTGAAGGCGTAGAGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGaa  <  1:28098/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                           
CGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAA  >  CP009273/3949257‑3949527

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑