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

CGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCA  >  NC_000913/1626694‑1626946                                                                                                                                    |                                                                                                                          cgcCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccg                                                                                                                    <  2:730614/139‑1 (MQ=255)                 gtATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTtgctg                                                                                                    >  2:334492/1‑139 (MQ=255)                  tATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTtgctgc                                                                                                   >  1:60155/1‑139 (MQ=255)                       agcagcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTGCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTa                                                                                              >  2:777970/1‑139 (MQ=255)                           gcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                          <  2:60155/139‑1 (MQ=255)                           gcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                          <  2:517027/139‑1 (MQ=255)                            cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgcc                                                                                                                  >  2:60015/1‑114 (MQ=255)                            cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgcc                                                                                                                  <  1:60015/114‑1 (MQ=255)                            cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGc                                                                                         >  1:1067861/1‑139 (MQ=255)                              gAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTg                                                                                       <  1:377622/139‑1 (MQ=255)                              gAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTg                                                                                       <  1:777970/139‑1 (MQ=255)                                 tCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGtt                                                                                    >  2:1003174/1‑139 (MQ=255)                                               tCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCa                                                                      <  1:1003174/139‑1 (MQ=255)                                                 cAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAgg                                                                    <  2:1067861/139‑1 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTc                                                                 >  1:1160029/1‑139 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                             <  1:501591/139‑1 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAg                                                             <  2:22022/139‑1 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGt                                                                                               <  1:777911/104‑1 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACAGCCGGATTTATTTGCTGCCAGt                                                                                               >  2:777911/1‑104 (MQ=255)                                                             gcgAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCt                                                        <  2:74450/139‑1 (MQ=255)                                                                            cgccAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCa                                                 >  2:1163184/1‑131 (MQ=255)                                                                            cgccAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCa                                                 <  1:1163184/131‑1 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccg                                                                                                                    >  2:30023/1‑57 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccg                                                                                                                    <  1:30023/57‑1 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTg                                  <  2:736866/139‑1 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTg                                  <  1:334492/139‑1 (MQ=255)                                                                                                   tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                           >  2:659735/1‑66 (MQ=255)                                                                                                   tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                           <  2:24393/66‑1 (MQ=255)                                                                                                   tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                           >  1:24393/1‑66 (MQ=255)                                                                                                   tcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                           <  1:659735/66‑1 (MQ=255)                                                                                                               gCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                        >  2:396980/1‑73 (MQ=255)                                                                                                               gCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                        <  1:396980/73‑1 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTgcgc                  >  2:1144571/1‑125 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTgcgc                  <  1:1144571/125‑1 (MQ=255)                                                                                                                   cgcgATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCa  >  1:611598/1‑139 (MQ=255)                                                                                                                          tCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGc                                                               <  2:304185/71‑1 (MQ=255)                                                                                                                          tCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGc                                                               >  1:304185/1‑71 (MQ=255)                                                                                                                                    |                                                                                                                          CGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCA  >  NC_000913/1626694‑1626946

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