Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	255,659	A→G	C20R (TGT→CGT)	pepD ←	aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	255,659	0	A	G	100.0%	209.6 / NA	58	C20R (TGT→CGT)	pepD	aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/25);  total (33/25)

CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGA  >  NC_000913/255530‑255793                                                                                                                                  |                                                                                                                                         cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                 >  2:353264/1‑139 (MQ=255)                     aGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCt                                                                                                           <  2:137173/139‑1 (MQ=255)                                 acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAAtt                                                                                               >  2:254605/1‑139 (MQ=255)                                      aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                >  1:1261858/1‑117 (MQ=255)                                      aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                <  2:1261858/117‑1 (MQ=255)                                      aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGAc                                                                                          >  2:1223819/1‑139 (MQ=255)                                            ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                >  2:800473/1‑111 (MQ=255)                                            ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                <  1:800473/111‑1 (MQ=255)                                              cGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                  >  1:721880/1‑139 (MQ=255)                                                     gAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                    >  1:1204876/1‑130 (MQ=255)                                                     gAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                    <  2:1204876/130‑1 (MQ=255)                                                     gAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTa                                                                           >  1:786249/1‑139 (MQ=255)                                                       aaCCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg                                                                         <  1:353264/139‑1 (MQ=255)                                                        aCCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                        <  2:721880/139‑1 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                 >  2:395286/1‑65 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                 <  1:395286/65‑1 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaa                                                                                                                               >  2:284036/1‑67 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaa                                                                                                                               <  1:284036/67‑1 (MQ=255)                                                                                     aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAt                                           >  2:371626/1‑139 (MQ=255)                                                                                          gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTgaga                                                                                           >  2:1334120/1‑86 (MQ=255)                                                                                          gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTgaga                                                                                           <  1:1334120/86‑1 (MQ=255)                                                                                             aGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                                   <  1:254605/139‑1 (MQ=255)                                                                                              gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                        <  1:378863/101‑1 (MQ=255)                                                                                              gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                        >  2:378863/1‑101 (MQ=255)                                                                                               ttGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     >  1:52750/1‑71 (MQ=255)                                                                                               ttGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     <  2:52750/71‑1 (MQ=255)                                                                                                tGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                  >  2:1492080/1‑89 (MQ=255)                                                                                                tGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                  <  1:1492080/89‑1 (MQ=255)                                                                                                          gATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                  <  1:1234465/79‑1 (MQ=255)                                                                                                          gATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                  >  2:1234465/1‑79 (MQ=255)                                                                                                            tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     >  2:869266/1‑58 (MQ=255)                                                                                                            tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     <  1:869266/58‑1 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                <  1:370934/44‑1 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAg                                                                                                                >  2:370934/1‑44 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     >  1:1358024/1‑55 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGgag                                                                                                     <  2:1358024/55‑1 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGc                        <  1:1340916/132‑1 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGc                        >  2:1340916/1‑132 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                 >  2:973837/1‑139 (MQ=255)                                                                                                                 cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt               >  1:790992/1‑139 (MQ=255)                                                                                                                 cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt               >  1:1087827/1‑139 (MQ=255)                                                                                                                   ggTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                   >  1:589757/1‑85 (MQ=255)                                                                                                                   ggTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                   <  2:589757/85‑1 (MQ=255)                                                                                                                   ggTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                 <  2:323387/119‑1 (MQ=255)                                                                                                                   ggTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                 >  1:323387/1‑119 (MQ=255)                                                                                                                   ggTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttg             >  1:60394/1‑139 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgtt              >  1:1282214/1‑139 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCgg                                                                                                             >  1:860485/1‑42 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCgg                                                                                                             <  2:860485/42‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcaca                                >  2:727266/1‑119 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcaca                                <  1:727266/119‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt            <  1:371626/139‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt            >  2:1245489/1‑139 (MQ=255)                                                                                                                         gAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                                   <  2:985225/111‑1 (MQ=255)                                                                                                                         gAATAGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGAAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                                   >  1:985225/1‑111 (MQ=255)                                                                                                                             aGAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAg   >  1:740089/1‑139 (MQ=255)                                                                                                                              gAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGa  <  2:740089/139‑1 (MQ=255)                                                                                                                              gAACGGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGa  >  2:537372/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                         CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCC‑‑CACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGA  >  NC_000913/255530‑255793

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: ‑