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

TATCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATCAGGCATTTGCGGGT  >  NC_000913/1626993‑1627235                                                                                                                                    |                                                                                                                tATCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGt                                                                                                          <  1:973810/139‑1 (MQ=255)           cAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCa                                                                                                >  1:1025900/1‑139 (MQ=255)             gggATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGt                                                                                              >  1:51205/1‑139 (MQ=255)                  tCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCg                                                                                         >  2:536253/1‑139 (MQ=255)                              cAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGc                                                                             >  1:139169/1‑139 (MQ=255)                                       tCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTGCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcgc                                                                    >  2:718755/1‑139 (MQ=255)                                       tCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcgc                                                                    >  1:154568/1‑139 (MQ=255)                                         tACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCgg                                                                  <  1:718755/139‑1 (MQ=255)                                               cgcgAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                                     >  2:879857/1‑130 (MQ=255)                                               cgcgAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAAcg                                                                     <  1:879857/130‑1 (MQ=255)                                               cgcgAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgtaagt                                                            >  2:893569/1‑139 (MQ=255)                                                 cgAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGgtaagtaa                                                          >  2:578098/1‑139 (MQ=255)                                                   aataatTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCaaaaaa                                                                                  <  1:454216/113‑1 (MQ=255)                                                   aataatTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCaaaaaa                                                                                  >  2:454216/1‑113 (MQ=255)                                                   aataatTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAAtt                                                        >  1:942696/1‑139 (MQ=255)                                                   aataatTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAAtt                                                        >  2:120492/1‑139 (MQ=255)                                                      aatTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATtctc                                                     <  1:816181/139‑1 (MQ=255)                                                       atTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTcc                                                    <  1:536253/139‑1 (MQ=255)                                                                tCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTc                                           <  1:1110778/139‑1 (MQ=255)                                                                tCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTc                                           >  2:1110778/1‑139 (MQ=255)                                                                 cAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTcc                                          >  1:1094311/1‑139 (MQ=255)                                                                     tagataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTg                                      <  2:942696/139‑1 (MQ=255)                                                                      agataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                                     <  2:1094311/139‑1 (MQ=255)                                                                      agataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                                     <  1:785525/139‑1 (MQ=255)                                                                      agataTCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGc                                     <  1:120492/139‑1 (MQ=255)                                                                        atatCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCa                                   <  1:893569/139‑1 (MQ=255)                                                                              ttAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAg                             >  1:588513/1‑139 (MQ=255)                                                                                      ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAg                                        >  2:570267/1‑120 (MQ=255)                                                                                      ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAg                                        <  1:570267/120‑1 (MQ=255)                                                                                      ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTAttgttg                     >  2:195475/1‑139 (MQ=255)                                                                                      ttCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTAtgttta                     >  2:1007247/1‑134 (MQ=255)                                                                                              gTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATCAgg             <  2:139169/139‑1 (MQ=255)                                                                                                         aaCTGCTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATCAGGCATTTGCGGGt  >  2:265475/1‑139 (MQ=255)                                                                                                              cTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGa                    >  2:401181/1‑116 (MQ=255)                                                                                                              cTCGTCAAGACGCTGTAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGa                    <  1:401181/116‑1 (MQ=255)                                                                                                                             tAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATtctc                                                     <  2:1062406/68‑1 (MQ=255)                                                                                                                             tAATTCACCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATtctc                                                     >  1:1062406/1‑68 (MQ=255)                                                                                                                                    |                                                                                                                TATCAAGCCCCAGGGATTCACGGCGCTGCCAGACAGCATCTACCCGCGCGAATAATTCACCGTTCAGATAGATATCATTAGCCAGTTCCGCCAGTTCAGCGGAAAACTGCTCGTCAAGACGCTGTAATTCATCATTGGTATGCGCCGCAGTCATCGCAAAAAAGACGCTGGTAACGCGGGTAAGTAATTCTCCGCTTTGTTCCAGTGCCAGAATAGTATTGTTGAAATCAGGCATTTGCGGGT  >  NC_000913/1626993‑1627235

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