Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,626,996 C→T D129N (GAT→AAT)  dcp ← peptidyl‑dipeptidase Dcp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,626,9960CT100.0% 42.2 / NA 13D129N (GAT→AAT) dcppeptidyl‑dipeptidase Dcp
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/8);  total (5/8)

GACAAAACGTTGATGAATCACCTCCACCAGGCGGATGGATTCACTATCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCC  >  NC_000913/1626949‑1627037
                                               |                                         
gACAAAACGTTGATGAATCACCTCCACCAGGCGGATGGATTCACTATTaa                                         <  1:1632494/50‑1 (MQ=255)
           gatgaATCACCTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGAt                              >  1:382112/1‑50 (MQ=255)
             tgaATCACCTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTc                            <  1:1414153/50‑1 (MQ=255)
              gaATCACCTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTCa                           >  1:2642/1‑50 (MQ=255)
               aaTCACCTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTCa                           <  1:1934091/49‑1 (MQ=255)
               aaTCACCTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTCAc                          >  1:968830/1‑50 (MQ=255)
                     cTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTCACGGCGCt                    <  1:1204455/50‑1 (MQ=255)
                     cTCCACCAGGCGGATGGATTCACTATTAAGCCCCAGGGATTCACGGCGCt                    <  1:1778153/50‑1 (MQ=255)
                            aGGCGGATGGATTCACTATTAAGCCCCAGGGATTCACGGCGCTGCCAGa              >  1:1027416/1‑49 (MQ=255)
                            aGGCGGATGGATTCACTATTAAGCCCCAGGGATTCACGGCGCTGCCAGAc             <  1:1088814/50‑1 (MQ=255)
                              gCGGATGGATTCACTATTAAGCCCCAGGGATTCACGGCGCTGCCAGACAg           <  1:1868079/50‑1 (MQ=255)
                                 gatggatTCACTATTAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCAt        >  1:2196386/1‑50 (MQ=255)
                                       ttCACTATTAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTAccc  <  1:155803/50‑1 (MQ=255)
                                               |                                         
GACAAAACGTTGATGAATCACCTCCACCAGGCGGATGGATTCACTATCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCC  >  NC_000913/1626949‑1627037

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