Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,233,504	T→A	V470D (GTT→GAT)	prkC →	non‑specific serine/threonine protein kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,233,504	0	T	A	100.0%	52.1 / NA	18	V470D (GTT→GAT)	prkC	non‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/9);  total (9/9)

AATACTGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCA  >  CP000730/1233362‑1233633                                                                                                                                               |                                                                                                                                  aaTACTGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGa                                                                                                                              >  1:44189/1‑148 (MQ=255)        gTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaaaa                                                                                                                       >  1:296667/1‑148 (MQ=255)           aaCGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGt                                                                                                                    <  2:233456/148‑1 (MQ=255)             cGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtat                                                                                                                  <  1:125188/148‑1 (MQ=255)                  tGAACGTGGTGACAGTGTTGAGGTTGTTATATCAAAGGGCCCTGAAAAGGGTAAACTGCCAAATGTCATTGGTTTACCTTAGGAGGAGGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATataat                                                                                                             >  2:180436/1‑148 (MQ=255)                      cGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATc                                                                                                         >  2:209445/1‑148 (MQ=255)                                     gATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                                                                          >  1:49903/1‑148 (MQ=255)                                                       ggCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGGTGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAgtg                                                                        >  2:238830/1‑148 (MQ=255)                                                               aaaGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                                                                <  2:296667/148‑1 (MQ=255)                                                                      aaaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                                         <  2:19788/148‑1 (MQ=255)                                                                      aaaaTGCCAAATGTCATAGGTTTACCTAAGGAGGAACCCTAGCAGAAATTAAAATCGTTAGGTCTTAAAGAAAAAAAAAAAAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                                         <  2:279231/148‑1 (MQ=255)                                                                           gCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATc                                                    >  1:34977/1‑148 (MQ=255)                                                                            ccAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCg                                                   <  2:72129/148‑1 (MQ=255)                                                                                aTGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAAAAATTAAAATCGTTAGGTCTTAAAGAAAAAAAAAAAAAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTAt                                               <  2:72218/148‑1 (MQ=255)                                                                                                   ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGCTCCATTGCAACTAAACGTGTAGACGACACTACTGAACTCGCTATTCATGATTCTAATATTaaa                            >  1:258750/1‑148 (MQ=255)                                                                                                   ggaggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaaa                            >  1:237666/1‑148 (MQ=255)                                                                                                                  gAAATTAAAATCGTTAGGGCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTa             <  1:19833/148‑1 (MQ=255)                                                                                                                             cGTTAGGTCTTAAAGATAATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCa  <  2:271988/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                  AATACTGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCA  >  CP000730/1233362‑1233633

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