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

GGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987284‑987445
                                                                                 |                                                                                
gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAt                                                                          >  2:758748/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  2:1175738/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  2:386639/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:138744/1‑90 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:643427/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:6118/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:928633/90‑1 (MQ=255)
                  gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:443565/1‑90 (MQ=255)
                       tCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                   >  1:1051782/1‑90 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  2:492043/90‑1 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  1:930602/90‑1 (MQ=255)
                            atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  1:248274/90‑1 (MQ=255)
                              ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            <  1:347011/90‑1 (MQ=255)
                                tAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                          <  2:575988/90‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTcc                                                             >  1:665715/1‑70 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTcc                                                             <  2:665715/70‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  2:927208/75‑1 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:927208/1‑75 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         >  2:797156/1‑90 (MQ=255)
                                 aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         <  1:300272/90‑1 (MQ=255)
                                  gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  1:1082597/1‑90 (MQ=255)
                                  gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  2:1213787/1‑90 (MQ=255)
                                    gCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                      <  1:1241307/90‑1 (MQ=255)
                                        cGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTacac                                  <  1:386639/90‑1 (MQ=255)
                                        cGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTacac                                  >  2:350230/1‑90 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        >  1:1183634/1‑83 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCg                                        <  2:1183634/83‑1 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGCACAcc                                 >  2:1013418/1‑90 (MQ=255)
                                         gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAACTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                 >  2:943179/1‑90 (MQ=255)
                                            cTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAt                              >  2:824413/1‑90 (MQ=255)
                                                 acgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                         <  2:517878/90‑1 (MQ=255)
                                                   gacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGAcc                       <  1:874402/90‑1 (MQ=255)
                                                        tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  <  1:453799/90‑1 (MQ=255)
                                                          tgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAg                <  2:138744/90‑1 (MQ=255)
                                                             aGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacg             <  1:824413/90‑1 (MQ=255)
                                                              ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga            >  2:13141/1‑90 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            >  2:655478/1‑53 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            <  1:655478/53‑1 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         >  2:1157697/1‑56 (MQ=255)
                                                                   tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                         <  1:1157697/56‑1 (MQ=255)
                                                                        gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg  >  2:168124/1‑90 (MQ=255)
                                                                        gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg  <  2:1114960/90‑1 (MQ=255)
                                                                                 |                                                                                
GGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987284‑987445

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: