Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,626,825	C→T	V186I (GTT→ATT)	dcp ←	dipeptidyl carboxypeptidase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,626,825	0	C	T	100.0%	131.3 / NA	41	V186I (GTT→ATT)	dcp	dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/19);  total (22/19)

GCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCT  >  NC_000913/1626695‑1626940                                                                                                                                   |                                                                                                                    gcCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgc                                                                                                             <  2:1296814/139‑1 (MQ=255)           gggTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGatttatt                                                                                                   >  2:135090/1‑139 (MQ=255)                           cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGc                                                                                   <  1:903911/139‑1 (MQ=255)                             gAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTg                                                                                 <  2:1034190/139‑1 (MQ=255)                                     cATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAAct                                                                         >  2:267317/1‑139 (MQ=255)                                        ttGTTATCCAGACCTTTCTCGCGAACCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActggc                                                                      >  1:1001651/1‑139 (MQ=255)                                              tCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCa                                                                >  1:1336416/1‑139 (MQ=255)                                                  gACCTTTCTCGCGAACCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGt                                                            <  2:1001651/139‑1 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                    <  2:1351859/114‑1 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                    >  1:203163/1‑114 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                    >  1:1351859/1‑114 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                    <  2:203163/114‑1 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCa                                                                >  2:140311/1‑134 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCa                                                                <  1:140311/134‑1 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTc                                                           >  1:91147/1‑139 (MQ=255)                                                   aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTc                                                           >  1:975695/1‑139 (MQ=255)                                                            gcgAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCt                                                  <  2:1336416/139‑1 (MQ=255)                                                               aGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGatttattt                                                                                                  <  2:1146230/88‑1 (MQ=255)                                                               aGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGatttattt                                                                                                  >  1:1146230/1‑88 (MQ=255)                                                                  cgccTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTc                                            <  2:975695/139‑1 (MQ=255)                                                                         gccgccAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACtt                                     >  2:852338/1‑139 (MQ=255)                                                                          ccgccAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGACTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACttt                                    >  1:479120/1‑139 (MQ=255)                                                                              cAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAAt                                <  1:135090/139‑1 (MQ=255)                                                                               aGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAAtt                               >  2:1337874/1‑139 (MQ=255)                                                                                                  tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                       >  1:317005/1‑96 (MQ=255)                                                                                                  tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                       <  2:317005/96‑1 (MQ=255)                                                                                                  tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTgcgc            >  2:709090/1‑139 (MQ=255)                                                                                                    aTTCCTGCCAGCTGCGCGATAACGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCaa          >  1:59781/1‑139 (MQ=255)                                                                                                        cTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGc                                                         >  2:418381/1‑88 (MQ=255)                                                                                                        cTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGc                                                         <  1:418381/88‑1 (MQ=255)                                                                                                          gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTaaa                                                                           <  1:835872/68‑1 (MQ=255)                                                                                                          gCCAGCTACGCGATATCGTTCACAGTCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTaaa                                                                           >  2:835872/1‑68 (MQ=255)                                                                                                            cAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCt  <  1:1337874/139‑1 (MQ=255)                                                                                                                tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                       >  2:920374/1‑66 (MQ=255)                                                                                                                tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                       <  1:920374/66‑1 (MQ=255)                                                                                                                            ttCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                       >  2:881375/1‑70 (MQ=255)                                                                                                                            ttCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                       <  1:881375/70‑1 (MQ=255)                                                                                                                            ttCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTAt                      <  2:569347/103‑1 (MQ=255)                                                                                                                            ttCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTAt                      >  1:569347/1‑103 (MQ=255)                                                                                                                              cacaATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                  <  2:107771/57‑1 (MQ=255)                                                                                                                              cacaATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                  >  1:107771/1‑57 (MQ=255)                                                                                                                                   |                                                                                                                    GCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCT  >  NC_000913/1626695‑1626940

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