Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,075,672 A→T H134L (CAC→CTC)  loiP → metalloprotease LoiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,075,6720AT100.0% 95.6 / NA 25H134L (CAC→CTC) loiPmetalloprotease LoiP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (10/15);  total (10/15)

CTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGCGAGTAGCTGCGGCCTCTGCGGGCGGG  >  NZ_CP009273/3075585‑3075757
                                                                                       |                                                                                     
cTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGt                                                                                    <  1:692263/91‑1 (MQ=255)
cTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGt                                                                                    <  2:389878/91‑1 (MQ=255)
cTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGt                                                                                    <  2:523234/91‑1 (MQ=255)
cTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGt                                                                                    <  1:505716/91‑1 (MQ=255)
      gCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGtt                                                                              <  2:634356/91‑1 (MQ=255)
                agTCTATAGCGGGCTGATGGATATGATGACGGATAAAGCAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGGTAGGCCAtgtg                                                                    >  2:235469/2‑91 (MQ=255)
                      taGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGaaaaaa                                                              <  1:28704/91‑1 (MQ=255)
                      taGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGaaaaaa                                                              <  1:793462/91‑1 (MQ=255)
                          gggCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGaa                                                          >  2:577812/1‑91 (MQ=255)
                               gATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAg                                                     >  1:195820/1‑91 (MQ=255)
                                         atgaCGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           <  1:707359/91‑1 (MQ=255)
                                         atgaCGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           <  2:237229/91‑1 (MQ=255)
                                         atgaCGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           >  1:113420/1‑91 (MQ=255)
                                         ataaCGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAACGCAGGTGGCACTTg                                           <  1:745916/88‑1 (MQ=255)
                                             cGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTAc                                       >  2:550668/1‑91 (MQ=255)
                                                 tAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCAtg                                                                      >  1:457347/1‑56 (MQ=255)
                                                 tAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCAtg                                                                      <  2:457347/56‑1 (MQ=255)
                                                  aaCGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATg                                  >  1:392364/1‑91 (MQ=255)
                                                        gTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           <  1:120417/76‑1 (MQ=255)
                                                        gTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           >  2:191419/1‑76 (MQ=255)
                                                        gTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           >  2:120417/1‑76 (MQ=255)
                                                        gTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTg                                           <  1:191419/76‑1 (MQ=255)
                                                        gTCGAAGCGGTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGc                            >  2:641011/1‑91 (MQ=255)
                                                                 gTGATCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGCGAGTAGCTg                   <  2:161476/91‑1 (MQ=255)
                                                                    aTCGGTCACGAAATGGGGCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGCGAGTAGCTGCgg                <  1:647857/91‑1 (MQ=255)
                                                                                  ggggCTCGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGCGAGTAGCTGCGGCCTCTGcgggcggg  >  2:763583/1‑91 (MQ=255)
                                                                                       |                                                                                     
CTAACGGCTGTATCCGCGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCGAAGCGGTGATCGGTCACGAAATGGGGCACGTGGCGTTAGGCCATGTGAAAAAAGGAATGCAGGTGGCACTTGGTACAAATGCCGTGCGAGTAGCTGCGGCCTCTGCGGGCGGG  >  NZ_CP009273/3075585‑3075757

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