Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 987,365 A→G N430D (AAT→GAT)  pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273987,3650AG100.0% 101.2 / NA 32N430D (AAT→GAT) pepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/17);  total (15/17)

GGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAAC  >  NZ_CP009273/987285‑987450
                                                                                |                                                                                     
ggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAtc                                                                              <  2:1451357/90‑1 (MQ=255)
         tttATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCAtt                                                                     >  1:713085/1‑90 (MQ=255)
             tttCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTccg                                                                 <  1:283280/90‑1 (MQ=255)
                 gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                             >  2:358325/1‑90 (MQ=255)
                      tCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                        >  2:1237266/1‑90 (MQ=255)
                         tgatgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAg                                                     <  1:220415/90‑1 (MQ=255)
                             ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                                 <  2:970915/90‑1 (MQ=255)
                             ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                                 >  2:339451/1‑90 (MQ=255)
                             ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                                 <  1:358325/90‑1 (MQ=255)
                              gTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGt                                                >  2:965933/1‑90 (MQ=255)
                                aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                              >  2:944541/1‑90 (MQ=255)
                                   gCAGCGCCCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                           <  2:1147742/90‑1 (MQ=255)
                                     aGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTac                                         <  1:1149316/90‑1 (MQ=255)
                                      gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                               <  2:119966/83‑1 (MQ=255)
                                      gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                               >  1:119966/1‑83 (MQ=255)
                                      gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTaca                                        >  1:634661/1‑90 (MQ=255)
                                      gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTaca                                        >  1:528847/1‑90 (MQ=255)
                                       cGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTacac                                       <  2:1279597/90‑1 (MQ=255)
                                        gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                      >  1:646597/1‑90 (MQ=255)
                                                 cgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGAc                             <  2:229399/90‑1 (MQ=255)
                                                       tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                       >  1:51314/1‑90 (MQ=255)
                                                       tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                       <  1:1434613/90‑1 (MQ=255)
                                                       tgtgCAGGCGATGGAAGATGCGTCGGATGTCGACCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                       <  2:1562983/90‑1 (MQ=255)
                                                             ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                           >  2:1540072/1‑64 (MQ=255)
                                                             ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                           <  1:1540072/64‑1 (MQ=255)
                                                                  tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACa            <  1:346742/90‑1 (MQ=255)
                                                                       gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg       <  2:1392677/90‑1 (MQ=255)
                                                                       gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg       <  2:416423/90‑1 (MQ=255)
                                                                       gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg       >  2:606403/1‑90 (MQ=255)
                                                                           cGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaaa   >  2:831837/1‑90 (MQ=255)
                                                                           cGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaaa   >  1:82632/1‑90 (MQ=255)
                                                                            gTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAAc  <  1:170557/90‑1 (MQ=255)
                                                                                |                                                                                     
GGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAAC  >  NZ_CP009273/987285‑987450

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: