Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 250,930 A→G I425T (ATC→ACC)  pepD ← cytosol nonspecific dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273250,9300AG100.0% 99.0 / NA 26I425T (ATC→ACC) pepDcytosol nonspecific dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/9);  total (17/9)

ACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAG  >  NZ_CP009273/250797‑251060
                                                                                                                                     |                                                                                                                                  
aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTc                                                                                                                               >  2:369579/1‑139 (MQ=255)
aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTc                                                                                                                               >  1:63017/1‑139 (MQ=255)
              tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGaa                                                                                                                 <  1:181409/139‑1 (MQ=255)
                        cTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGa                                                                                                       >  2:231532/1‑139 (MQ=255)
                           aTCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                                      <  2:283150/137‑1 (MQ=255)
                                     tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTa                                                                                          >  2:163810/1‑139 (MQ=255)
                                     tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTAAGAAAAAGTTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTa                                                                                          >  2:440613/1‑139 (MQ=255)
                                       tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                         >  2:382002/1‑138 (MQ=255)
                                       tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                         >  2:68837/1‑138 (MQ=255)
                                               tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAgg                                                                                                             <  1:428218/110‑1 (MQ=255)
                                               tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAgg                                                                                                             >  2:428218/1‑110 (MQ=255)
                                               tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                                >  1:302290/1‑139 (MQ=255)
                                               tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                                >  2:427017/1‑139 (MQ=255)
                                                    gTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCg                                                                           <  1:289184/139‑1 (MQ=255)
                                                                   aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                            >  2:20747/1‑107 (MQ=255)
                                                                   aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                            <  1:20747/107‑1 (MQ=255)
                                                                        gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTc                                                                                                                         >  1:51344/1‑73 (MQ=255)
                                                                        gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTc                                                                                                                         <  2:51344/73‑1 (MQ=255)
                                                                        gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                       >  1:53986/1‑139 (MQ=255)
                                                                                                            ttCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcgcg                   <  1:163810/139‑1 (MQ=255)
                                                                                                            tgccAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcgtg                   <  2:37155/137‑3 (MQ=255)
                                                                                                                     cgcgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGt          >  2:59969/1‑139 (MQ=255)
                                                                                                                       cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt         >  2:339174/1‑138 (MQ=255)
                                                                                                                       cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt         >  2:219999/1‑138 (MQ=255)
                                                                                                                       cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt         >  1:115003/1‑138 (MQ=255)
                                                                                                                             tAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAg  <  1:219999/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
ACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAG  >  NZ_CP009273/250797‑251060

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