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% 216.9 / NA 59N430D (AAT→GAT) pepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/24);  total (35/24)

AGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAA  >  NZ_CP009273/987279‑987449
                                                                                      |                                                                                    
aGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATg                                                                                   >  2:1015290/1‑90 (MQ=255)
  aaaGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTc                                                                                 <  2:2607510/90‑1 (MQ=255)
     gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCg                                                                                <  1:2864708/88‑1 (MQ=255)
     gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCg                                                                                >  2:2864708/1‑88 (MQ=255)
     gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAt                                                                              >  2:2005852/1‑90 (MQ=255)
      ggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAtc                                                                             <  2:1415005/90‑1 (MQ=255)
             gCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCa                                                                      >  2:1495391/1‑90 (MQ=255)
                       gAGCGTCATGCTGGTAGTGCAGCGACCTGTGCCGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAGCTCTCCCATTTCCGCCGt                                                            <  1:1675742/90‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTccgcc                                                              <  1:1509668/88‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTccgcc                                                              >  2:1509668/1‑88 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  1:1668462/90‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  2:2827417/1‑90 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  2:1388197/1‑90 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  1:1495391/90‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  1:1492935/90‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  1:796822/1‑90 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  2:452803/1‑90 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  2:667859/90‑1 (MQ=255)
                       gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGGCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  1:1776830/1‑90 (MQ=255)
                                 atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                                  <  1:1848328/90‑1 (MQ=255)
                                 atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                                  <  2:2757846/90‑1 (MQ=255)
                                 atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                                  <  1:1543680/90‑1 (MQ=255)
                                    gTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCAATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGt                                               >  1:2755286/1‑90 (MQ=255)
                                     tAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                              >  1:2545248/1‑90 (MQ=255)
                                      aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                             >  2:2989895/1‑90 (MQ=255)
                                       gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  2:833749/1‑74 (MQ=255)
                                       gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  1:833749/74‑1 (MQ=255)
                                       gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                            >  1:1885643/1‑90 (MQ=255)
                                       gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                            >  2:2129386/1‑90 (MQ=255)
                                         gCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                          >  2:432940/1‑90 (MQ=255)
                                         gCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                          >  1:2906206/1‑90 (MQ=255)
                                           aGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTac                                        <  1:8204/90‑1 (MQ=255)
                                            gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTaca                                       >  2:1414196/1‑90 (MQ=255)
                                            gCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTaca                                       >  1:2819620/1‑90 (MQ=255)
                                              gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                     >  1:1563184/1‑90 (MQ=255)
                                              gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                     >  1:1804703/1‑90 (MQ=255)
                                                 cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                       >  2:1917114/1‑69 (MQ=255)
                                                 cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                       <  1:1917114/69‑1 (MQ=255)
                                                 cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAt                                  >  2:2153231/1‑90 (MQ=255)
                                                 cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAt                                  >  2:101774/1‑90 (MQ=255)
                                                      acgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                             <  1:1377187/90‑1 (MQ=255)
                                                      acgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                             <  1:452803/90‑1 (MQ=255)
                                                       cgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGAc                            <  2:1827994/90‑1 (MQ=255)
                                                         acTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCg                          >  2:572397/1‑90 (MQ=255)
                                                             tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                      <  2:2435571/90‑1 (MQ=255)
                                                             tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                      <  2:404492/90‑1 (MQ=255)
                                                                   ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga                >  1:2997258/1‑90 (MQ=255)
                                                                   ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga                >  2:1166804/1‑90 (MQ=255)
                                                                   ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga                >  2:39672/1‑90 (MQ=255)
                                                                        tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACa           >  1:1278366/1‑90 (MQ=255)
                                                                        tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACa           <  1:1015290/90‑1 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            <  1:1886665/36‑1 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                            >  2:1886665/1‑36 (MQ=38)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                       >  2:2678720/1‑41 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                       <  1:2678720/41‑1 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg      <  1:1388197/90‑1 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg      >  2:1927830/1‑90 (MQ=255)
                                                                             gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg      >  1:880599/1‑90 (MQ=255)
                                                                                 cGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaaa  >  2:2634608/1‑90 (MQ=255)
                                                                                      |                                                                                    
AGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAA  >  NZ_CP009273/987279‑987449

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