Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,335,202 A→C H193P (CAC→CCC)  kdsD → arabinose‑5‑phosphate isomerase KdsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,335,2020AC100.0% 122.1 / NA 34H193P (CAC→CCC) kdsDarabinose‑5‑phosphate isomerase KdsD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/13);  total (21/13)

TGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATCCCGCATGTTAAGAA  >  NZ_CP009273/3335127‑3335286
                                                                           |                                                                                    
tGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTCACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCAc                                                                        <  1:75971/90‑1 (MQ=255)
          cGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaa                                                              <  1:234430/90‑1 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  2:304864/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  1:213158/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  2:283848/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  2:394724/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  2:140869/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  2:405191/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  1:440670/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaaa                                                             >  1:71229/1‑90 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaa                                                              <  1:172545/89‑1 (MQ=255)
           gATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTaa                                                              >  2:172545/1‑89 (MQ=255)
                  tCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTctgc                                                      >  2:351610/1‑90 (MQ=255)
                        tCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGcgcg                                                <  1:304864/90‑1 (MQ=255)
                        tCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGcgcg                                                >  1:215532/1‑90 (MQ=255)
                           cgcTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTaa                                             <  2:300933/90‑1 (MQ=255)
                             cTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAAc                                           <  1:86003/90‑1 (MQ=255)
                             cTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAAc                                           <  1:3480/90‑1 (MQ=255)
                             cTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAAc                                           <  2:215532/90‑1 (MQ=255)
                               gTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGa                                         >  1:499585/1‑90 (MQ=255)
                                    aaGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGAtatta                                    <  2:71229/90‑1 (MQ=255)
                                         cgcgGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTg                                                                      <  1:80648/51‑1 (MQ=255)
                                         cgcgGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTg                                                                      >  2:80648/1‑51 (MQ=255)
                                         cgcgGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCAt                               >  1:208451/1‑90 (MQ=255)
                                         cgcgGCTTTACTGCTGAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCAt                               >  1:327124/1‑90 (MQ=255)
                                                    tgctgAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATga                    <  2:440670/90‑1 (MQ=255)
                                                     gctgAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAAGTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATgag                   >  1:394234/1‑90 (MQ=255)
                                                      ctgAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATgaga                  <  1:405191/90‑1 (MQ=255)
                                                       tgAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGAt                 >  2:347477/1‑90 (MQ=255)
                                                       tgAAGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGAt                 >  2:253658/1‑90 (MQ=255)
                                                          aGATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATccc              >  2:503952/1‑90 (MQ=255)
                                                           gATTTTGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATCCCg             <  2:98254/90‑1 (MQ=255)
                                                               ttGCGCTCTCACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATCCCGCATg         >  2:356697/1‑90 (MQ=255)
                                                                      ctcACCCCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATCCCGCATGTTAAGaa  >  2:375895/1‑90 (MQ=255)
                                                                           |                                                                                    
TGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTAAACTTCTGCTGCGCGTAAACGATATTATGCATACGGGCGATGAGATCCCGCATGTTAAGAA  >  NZ_CP009273/3335127‑3335286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: