Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,506,028 G→A A40T (GCC→ACC)  ydcS → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,506,0280GA100.0% 73.1 / NA 22A40T (GCC→ACC) ydcSABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/13);  total (9/13)

AGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGA  >  NZ_CP009273/1505950‑1506115
                                                                              |                                                                                       
aGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGAt                                                                             >  1:531003/1‑91 (MQ=255)
         ataatGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAAc                                                                    <  2:121173/91‑1 (MQ=255)
              gACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGa                                                               >  1:285336/1‑91 (MQ=255)
                ccGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACa                                                             >  1:347360/1‑91 (MQ=255)
                 cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACaa                                                            <  1:86928/91‑1 (MQ=255)
                 cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACaa                                                            <  2:143845/91‑1 (MQ=255)
                   cTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAAc                                                          <  1:363131/91‑1 (MQ=255)
                         ccgccgAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATaa                                                    <  1:495951/91‑1 (MQ=255)
                         ccgccgAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATaa                                                    <  1:557375/91‑1 (MQ=255)
                           gccgAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAAc                                                  >  1:695462/1‑91 (MQ=255)
                                   gCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGAc                                          >  1:605523/1‑91 (MQ=255)
                                   gCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGAc                                          <  2:548969/91‑1 (MQ=255)
                                   gCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGAc                                          <  2:214816/91‑1 (MQ=255)
                                             ttAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGc                                <  2:285336/91‑1 (MQ=255)
                                                    aaCCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGa                         >  2:269377/1‑91 (MQ=255)
                                                              gCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACa               >  2:762759/1‑91 (MQ=255)
                                                                gACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACaa                                                            <  2:208721/44‑1 (MQ=255)
                                                                gACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACaa                                                            >  1:208721/1‑44 (MQ=255)
                                                                gACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAgg             >  2:241935/1‑91 (MQ=255)
                                                                      atatTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcgc       <  2:234583/91‑1 (MQ=255)
                                                                      atatTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcgc       <  1:435279/91‑1 (MQ=255)
                                                                           atCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGa  <  2:32034/91‑1 (MQ=255)
                                                                              |                                                                                       
AGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGA  >  NZ_CP009273/1505950‑1506115

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