Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,627,124	T→C	D86G (GAT→GGT)	dcp ←	dipeptidyl carboxypeptidase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,627,124	0	T	C	100.0%	102.2 / NA	29	D86G (GAT→GGT)	dcp	dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/14);  total (15/14)

GGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTT  >  NC_000913/1627005‑1627215                                                                                                                        |                                                                                            gggATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTGTGCGCCGCAGt                                                                          <  1:353920/139‑1 (MQ=255) gggATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                          <  1:238525/139‑1 (MQ=255) gggATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                          >  2:639714/1‑139 (MQ=255) gggATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                          >  2:29144/1‑139 (MQ=255)   gATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCa                                                                        >  1:651131/1‑139 (MQ=255)                         cATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAc                                                  <  1:639714/139‑1 (MQ=255)                          aTCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                 <  1:29144/139‑1 (MQ=255)                           tCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcgc                                                <  2:638020/139‑1 (MQ=255)                           tCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcgc                                                >  1:638020/1‑139 (MQ=255)                                                  gTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGt                         >  2:160504/1‑139 (MQ=255)                                                     cAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTcc                      <  2:651131/139‑1 (MQ=255)                                                        atagataTCATTAGCCAGTGCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGt                   <  2:115566/139‑1 (MQ=255)                                                          agataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                 <  1:160504/139‑1 (MQ=255)                                                           gataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGcc                >  2:720424/1‑139 (MQ=255)                                                             tatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAg              >  1:675571/1‑139 (MQ=255)                                                              atCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGCATGCGCCGCAGTCATCGCAAAAAAGACGCTGATAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGa             >  1:72382/1‑139 (MQ=255)                                                                    aGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTa       >  1:38098/1‑139 (MQ=255)                                                                                   ttCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAc                                                  >  2:424596/1‑81 (MQ=255)                                                                                   ttCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAc                                                  <  1:424596/81‑1 (MQ=255)                                                                                   ttCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgt                                            >  1:88261/1‑87 (MQ=255)                                                                                   ttCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgt                                            <  2:88261/87‑1 (MQ=255)                                                                                       gCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                          >  2:385330/1‑53 (MQ=255)                                                                                       gCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                          <  1:385330/53‑1 (MQ=255)                                                                                                  gtCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGaa            >  1:515999/2‑104 (MQ=255)                                                                                                  cTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGaa            <  2:515999/104‑1 (MQ=255)                                                                                                       cAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTAttgtt  <  1:527044/109‑1 (MQ=255)                                                                                                       cAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTAttgtt  >  2:527044/1‑109 (MQ=255)                                                                                                               tGTAATTCACCATTGGTATGCGCCGCAGTCATCGCaaaa                                                                >  2:265724/1‑39 (MQ=255)                                                                                                               tGTAATTCACCATTGGTATGCGCCGCAGTCATCGCaaaa                                                                <  1:265724/39‑1 (MQ=255)                                                                                                                        |                                                                                            GGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTT  >  NC_000913/1627005‑1627215

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑