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

AAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987282‑987445
                                                                                   |                                                                                
aaGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCg                                                                            <  2:649625/90‑1 (MQ=255)
    gATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGAtct                                                                        <  2:586262/90‑1 (MQ=255)
          gCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCa                                                                  >  2:179399/1‑90 (MQ=255)
                    gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        <  1:297150/90‑1 (MQ=255)
                    gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGt                                                        >  1:340546/1‑90 (MQ=255)
                    gAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCACTTCCGCCGt                                                        >  2:680548/1‑90 (MQ=255)
                        gTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGt                                                    <  1:94600/90‑1 (MQ=255)
                              atgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGcc                                              <  2:340546/90‑1 (MQ=255)
                                ggTAGTGCAGCGACCTGTGACGATTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAg                                            <  1:111245/90‑1 (MQ=255)
                                  tAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTc                                          >  1:329392/1‑90 (MQ=255)
                                  tAGGGCCGCGCCCTGTGCCGTTTTTGGCCGGGGGATGAAAGAGGGGTGGGATGTCGACCCCCCCCATCCCCGCCGTTGGTACAGCCAGTc                                          <  1:497596/90‑1 (MQ=255)
                                      gCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGt                                      >  2:424840/1‑90 (MQ=255)
                                           gACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                 >  2:554664/1‑90 (MQ=255)
                                                   acgacTTTGTGCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                         <  1:554664/90‑1 (MQ=255)
                                                          tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  >  1:664431/1‑90 (MQ=255)
                                                          tgtgCAGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  <  1:423577/90‑1 (MQ=255)
                                                               aGGCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacg             <  2:329392/90‑1 (MQ=255)
                                                                ggCGATGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga            >  2:562222/1‑90 (MQ=255)
                                                                     tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  <  2:26286/79‑1 (MQ=255)
                                                                     tGGAAGATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaa                  >  1:26286/1‑79 (MQ=255)
                                                                          gATGCGTCGGATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg  >  2:485795/1‑90 (MQ=255)
                                                                                   |                                                                                
AAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCG  >  NZ_CP009273/987282‑987445

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