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%	116.4 / NA	33	D86G (GAT→GGT)	dcp	dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/16);  total (17/16)

TCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATC  >  NC_000913/1626995‑1627221                                                                                                                                  |                                                                                                  tCAAGCCCCGGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATGGGTAt                                                                                          <  1:261531/139‑1 (MQ=255)    aGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTAt                                                                                          <  1:1157643/136‑1 (MQ=255)    aGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTAt                                                                                          >  2:1157643/1‑136 (MQ=255)    aGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATgcg                                                                                       >  2:907761/1‑139 (MQ=255)                        cTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAg                                                                   >  2:441438/1‑139 (MQ=255)                         tGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGa                                                                  <  1:1236705/139‑1 (MQ=255)                                     tCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTAGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcgc                                                      >  1:924661/1‑139 (MQ=255)                                                           cGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgtaag                                               <  2:1232508/124‑1 (MQ=255)                                                           cGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgtaag                                               >  1:1232508/1‑124 (MQ=255)                                                             ttCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                       <  1:536254/114‑1 (MQ=255)                                                             ttCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                       >  2:536254/1‑114 (MQ=255)                                                                    agataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                       <  1:333077/139‑1 (MQ=255)                                                                      atatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCa                     <  1:907761/139‑1 (MQ=255)                                                                       tatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgt                                                  <  1:913588/109‑1 (MQ=255)                                                                       tatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                       >  2:597724/1‑136 (MQ=255)                                                                       tatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                       <  1:597724/136‑1 (MQ=255)                                                                       tatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTAAGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgt                                                  >  2:913588/1‑109 (MQ=255)                                                                          cATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTcc                            <  2:562378/128‑1 (MQ=255)                                                                          cATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTcc                            >  1:562378/1‑128 (MQ=255)                                                                              aGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGt              <  1:844373/138‑1 (MQ=255)                                                                              aGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGt              >  2:844373/1‑138 (MQ=255)                                                                                    ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCgg                                                    <  1:1048447/94‑1 (MQ=255)                                                                                    ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCgg                                                    >  2:1048447/1‑94 (MQ=255)                                                                                    ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTAttgttg       >  2:347050/1‑139 (MQ=255)                                                                                         ccAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                       <  2:1020069/86‑1 (MQ=255)                                                                                         ccAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                       >  1:1020069/1‑86 (MQ=255)                                                                                         ccAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAAGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATc  >  2:1363475/1‑139 (MQ=255)                                                                                          cAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATg                                                                                         >  1:868085/1‑51 (MQ=255)                                                                                          cAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATg                                                                                         <  2:868085/51‑1 (MQ=255)                                                                                                 gCGGAAAACTGCTCGTCAAGACGCTGTAATtcacca                                                                                                 >  2:509891/1‑36 (MQ=255)                                                                                                 gCGGAAAACTGCTCGTCAAGACGCTGTAATtcacca                                                                                                 <  1:509891/36‑1 (MQ=255)                                                                                                       aaCTGCTCGTCAAGACGCTGTAATTCACCATTGGTAt                                                                                          <  2:987352/37‑1 (MQ=255)                                                                                                       aaCTGCTCGTCAAGACGCTGTAATTCACCATTGGTAt                                                                                          >  1:987352/1‑37 (MQ=255)                                                                                                                           tAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                       >  2:757061/1‑84 (MQ=255)                                                                                                                           tAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                       <  1:757061/84‑1 (MQ=255)                                                                                                                                  |                                                                                                  TCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATC  >  NC_000913/1626995‑1627221

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