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%	204.4 / NA	57	C20R (TGT→CGT)	pepD	aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/25);  total (32/25)

GAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACA  >  NC_000913/255545‑255788                                                                                                                   |                                                                                                                                  gAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                           <  2:276059/139‑1 (MQ=255)    tCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                        <  2:734815/139‑1 (MQ=255)            ttACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                <  1:791443/139‑1 (MQ=255)                  acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                          >  2:230619/1‑139 (MQ=255)                      gATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      <  1:230619/139‑1 (MQ=255)                       aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      >  2:586558/1‑138 (MQ=255)                       aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      <  1:586558/138‑1 (MQ=255)                       aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                     >  2:751373/1‑139 (MQ=255)                       aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                     >  2:560105/1‑139 (MQ=255)                       aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                     >  1:641660/1‑139 (MQ=255)                         cgcgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                   >  1:765530/1‑139 (MQ=255)                         cgcgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                   >  1:321871/1‑139 (MQ=255)                           cgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                 <  1:138650/139‑1 (MQ=255)                           cgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                 >  1:454393/1‑139 (MQ=255)                             ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                 >  2:776839/1‑137 (MQ=255)                             ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                 <  1:776839/137‑1 (MQ=255)                             ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                               >  1:14310/1‑139 (MQ=255)                             ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                               >  2:571271/1‑139 (MQ=255)                               cGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                             <  2:333931/139‑1 (MQ=255)                                gACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                      <  2:388538/113‑1 (MQ=255)                                gACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                      >  1:388538/1‑113 (MQ=255)                                     ggAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGGACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt                                                                       >  2:484824/1‑139 (MQ=255)                                         aCCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                   <  1:457702/139‑1 (MQ=255)                                           cTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGtct                                                                 <  2:14310/139‑1 (MQ=255)                                                           accaacAATGTATTCAGCGAGTTGCTCTTCATGATCGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                 <  2:485563/139‑1 (MQ=255)                                                           accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                 <  1:371058/139‑1 (MQ=255)                                                           accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                 >  1:480935/1‑139 (MQ=255)                                                           accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                 <  2:321871/139‑1 (MQ=255)                                                                      aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAt                                      >  2:755037/1‑139 (MQ=255)                                                                       ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCg                                                                                                         <  1:13378/71‑1 (MQ=255)                                                                       ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCg                                                                                                         >  2:13378/1‑71 (MQ=255)                                                                       ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg                                                                    >  1:400439/1‑108 (MQ=255)                                                                       ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg                                                                    <  2:400439/108‑1 (MQ=255)                                                                       ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAtt                                     >  2:648516/1‑139 (MQ=255)                                                                          aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                  >  2:198525/1‑139 (MQ=255)                                                                           gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      <  2:694552/86‑1 (MQ=255)                                                                           gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      >  1:694552/1‑86 (MQ=255)                                                                             gAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGt                               >  1:256083/1‑139 (MQ=255)                                                                              aGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                              >  1:606407/1‑139 (MQ=255)                                                                              aGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                              >  1:20389/1‑139 (MQ=255)                                                                               gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTca                             <  1:198525/139‑1 (MQ=255)                                                                                             tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGaa               >  1:523743/1‑139 (MQ=255)                                                                                                gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt            <  2:606407/139‑1 (MQ=255)                                                                                                gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt            >  1:265069/1‑139 (MQ=255)                                                                                                  cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt          >  1:609920/1‑139 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      <  2:260454/60‑1 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                      >  1:260454/1‑60 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                             <  2:518799/69‑1 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                             >  1:518799/1‑69 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                   <  2:566427/111‑1 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                   >  1:566427/1‑111 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt       >  2:595674/1‑139 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt       >  1:206914/1‑139 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt       <  2:765530/139‑1 (MQ=255)                                                                                                     gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt       <  1:648516/139‑1 (MQ=255)                                                                                                        aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTAAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCaa    >  2:592080/1‑139 (MQ=255)                                                                                                          gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca  <  2:480935/139‑1 (MQ=255)                                                                                                                   |                                                                                                                                  GAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACA  >  NC_000913/255545‑255788

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