Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,611	G→T	G156V (GGT→GTT)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,611	0	G	T	100.0%	149.4 / NA	39	G156V (GGT→GTT)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/23);  total (16/23)

TCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGG  >  NC_000913/4033534‑4033688                                                                              |                                                                              tCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCgtggttgg                                                                            <  1:743998/81‑1 (MQ=255) tCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCgtggttgg                                                                            <  2:743999/81‑1 (MQ=255)         tgctgcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGc                                                                    <  1:744000/81‑1 (MQ=255)          gctgcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGct                                                                   >  1:744055/1‑81 (MQ=255)           ctgcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGctc                                                                  >  2:744061/1‑81 (MQ=255)           ctgcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGctc                                                                  <  2:744002/81‑1 (MQ=255)           ctgcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGct                                                                   <  1:744001/80‑1 (MQ=255)             gcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTa                                                                >  2:744063/1‑81 (MQ=255)             gcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTa                                                                >  1:744057/1‑81 (MQ=255)             gcAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTa                                                                >  1:744049/1‑81 (MQ=255)              cAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGctct                                                                 <  1:744003/79‑1 (MQ=255)                 tGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATcgc                                                            <  1:744004/81‑1 (MQ=255)                     ttGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATcgcg                                                           >  2:744053/1‑78 (MQ=255)                        gCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGaa                                                      <  2:744005/80‑1 (MQ=255)                         cGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAAt                                                    <  2:744006/81‑1 (MQ=255)                          gggATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATg                                                   >  2:744062/1‑81 (MQ=255)                           ggATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATgc                                                  >  1:744067/1‑81 (MQ=255)                           ggATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATgc                                                  <  1:744007/81‑1 (MQ=255)                                 ggATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATgcccggcc                                            <  1:744008/81‑1 (MQ=255)                                  gATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATgcccggccc                                           <  1:744009/81‑1 (MQ=255)                                  gATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATgcccggccc                                           <  1:744010/81‑1 (MQ=255)                                    tcatcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGc                                         <  1:744011/81‑1 (MQ=255)                                      atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCt                                        >  1:744070/1‑80 (MQ=255)                                      atcGTGTTGGCGGTTGCGATACTGCCTATCCACGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTg                                       <  2:744012/81‑1 (MQ=255)                                        cGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGTTCTATCGCGCAGAAATGCCCGGCCCGCTGa                                      >  2:744075/1‑80 (MQ=255)                                         gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGa                                      >  1:744066/1‑79 (MQ=255)                                         gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGa                                      >  1:744058/1‑79 (MQ=255)                                                        aTACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATgcgc                     >  2:744064/1‑81 (MQ=255)                                                         tACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGcc                    >  1:744078/1‑81 (MQ=255)                                                         tACTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGcc                    <  2:744013/81‑1 (MQ=255)                                                           cTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGccc                   <  2:744014/80‑1 (MQ=255)                                                           cTGCCTATCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCcg                  <  2:744015/81‑1 (MQ=255)                                                                  tCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAtt            <  2:744017/80‑1 (MQ=255)                                                                  tCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAtt            <  1:744019/80‑1 (MQ=255)                                                                  tCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTg           <  2:744016/81‑1 (MQ=255)                                                                  tCCTCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTg           <  1:744018/81‑1 (MQ=255)                                                                     tCGGCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCg         >  2:744065/1‑80 (MQ=255)                                                                       ggCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCgg        <  2:744020/79‑1 (MQ=255)                                                                       ggCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGa       <  2:744021/80‑1 (MQ=255)                                                                        gCGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGaaa     >  2:744069/1‑81 (MQ=255)                                                                         cGTGGTTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAAc    <  2:744022/81‑1 (MQ=255)                                                                          gtggttggGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACg   <  2:744023/81‑1 (MQ=255)                                                                          gtggttggGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACg   <  1:744024/81‑1 (MQ=255)                                                                           tggttgggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACgg  <  2:744025/81‑1 (MQ=255)                                                                           tggttgggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACgg  >  1:744081/1‑81 (MQ=255)                                                                              |                                                                              TCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGG  >  NC_000913/4033534‑4033688

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