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%	165.7 / NA	56	V470D (GTT→GAT)	prkC	non‑specific serine/threonine protein kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (32/24);  total (32/24)

ATACTGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTATGTAGAA  >  CP000730/1233363‑1233646                                                                                                                                              |                                                                                                                                               aTACTGGTGAACGTGTTGAACGTGGGGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAt                                                                                                                                          <  2:867183/148‑1 (MQ=255)     tGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaa                                                                                                                                      >  1:729917/1‑148 (MQ=255)     tGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaa                                                                                                                                      >  1:1019452/1‑148 (MQ=255)     tGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGaa                                                                                                                                      >  1:448550/1‑148 (MQ=255)        tGAACATGTTGAACGTGGTGAGAGTGTTGATGTTGTTATATCAAAGGGCCATGAAAAGGTTAAAATGCCAAATGTCACTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGAGACGATTGaaaaa                                                                                                                                   <  2:175816/148‑1 (MQ=255)         gAACGTGTTAAACGTGGTGAAAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTCCCTAAGGAGGAAGCCGTGCAGAAATTAAAATCGTTAGTTTTTAAAGATGATACGATTGAAAAAg                                                                                                                                  <  2:1239423/148‑1 (MQ=255)            cGTTTTGAACGTGGTGAAAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtat                                                                                                                               <  2:342131/148‑1 (MQ=255)            cGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtat                                                                                                                               <  2:48491/148‑1 (MQ=255)            cGTGTTGAACGTGGGAACAGTGTTGATGTTGTTAGATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACATAAGGAGGAAGCCCTGCAGAAATTAAAATCGTTAGGTTTTAAAGATGATACGATTGAAAAAGtat                                                                                                                               <  1:756390/148‑1 (MQ=255)              ggtTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGtatat                                                                                                                             <  2:815305/147‑1 (MQ=255)                                tgtTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCaaaa                                                                                                           >  1:364362/1‑148 (MQ=255)                                    taagtttttATATCAAAGGGCCCTTAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGAt                                                                                                       <  2:631477/145‑1 (MQ=255)                                    gATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAAGAt                                                                                                       >  1:245506/1‑148 (MQ=255)                                         tgttATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACAtt                                                                                                  >  1:814376/1‑148 (MQ=255)                                            tacatCAAAGGGCCCGTAAAAGGTTAAAACGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCa                                                                                               <  1:1209407/145‑1 (MQ=255)                                              tatCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCaaa                                                                                             >  1:218662/1‑148 (MQ=255)                                                       gCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAgtgt                                                                                    >  2:1253018/1‑148 (MQ=255)                                                        cccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGAATTGAAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTa                                                                                   <  1:1227740/148‑1 (MQ=255)                                                         ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGGGGAAGCCTTGGAGGAGTTAGAATCGTTAGGTATTAAAGATGATACGATTGAAAAAGTAGATAATAATCAAGCGCCGAAAGAGTAAATTGAAAATCAAAGTGaaa                                                                                  >  2:882142/1‑148 (MQ=255)                                                         ccTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATGAAAGTGTaa                                                                                  >  1:59899/1‑148 (MQ=255)                                                              aacgccatcaatgtcaaatgtcaTTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                                                                             >  2:946082/10‑148 (MQ=255)                                                              aaaGGTTAAAATGCCAAATGTCATTGGTTTACCCAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                                                                             >  2:774992/1‑148 (MQ=255)                                                              aaaGGTTAAAATGACAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATCATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCa                                                                             <  1:880020/148‑1 (MQ=255)                                                                aGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATCTAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCaaa                                                                           >  2:827468/1‑148 (MQ=255)                                                                  gTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATCCTAATCAAGCGCCAAAAGGATACATTGCAAATCAACGTGTAGCCGCAAATa                                                                         <  2:729917/148‑1 (MQ=255)                                                                     aaaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCGAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                                                      >  2:374567/1‑148 (MQ=255)                                                                     aaaaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTg                                                                      >  2:181952/1‑148 (MQ=255)                                                                       aaTGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGaa                                                                    >  2:125035/1‑148 (MQ=255)                                                                       aaTGCCAAATGACCTGGGTTTACCTAAGGAGGAAGCCTTGCAGAACTTAAAGTCGTAAGGTATGTGAGAGGAAACGATTGAAAAAGTAAAGAATAATAAAGCGCCAAAAGGAACAATTGCAACTCAAAGTGTAACCGCAACAACTGaa                                                                    >  2:139644/1‑148 (MQ=255)                                                                                   cATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCAt                                                        <  1:356345/148‑1 (MQ=255)                                                                                   cATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCAt                                                        <  1:491197/148‑1 (MQ=255)                                                                                          ttACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTa                                                 >  1:677006/1‑148 (MQ=255)                                                                                          ttACCTAAGGAGGAAGCCTTGCAAAAATTAAAATCGTTAGGTCTTAAAGATGATACGAGTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTa                                                 >  2:729593/1‑148 (MQ=255)                                                                                           tACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTaa                                                <  1:1050619/148‑1 (MQ=255)                                                                                                aaGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTa                                           >  1:1128615/1‑148 (MQ=255)                                                                                                aaGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTa                                           >  1:1128557/1‑148 (MQ=255)                                                                                                 aggaggAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTaa                                          >  2:539416/1‑148 (MQ=255)                                                                                                 aggaggAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAAGAATAAAGCGCCAAAAGGATACATTGCTAGTTAAATTGTAACCGCAAATACTAAAATCGCTATTCATCATGCTACTATTaa                                          >  1:411475/1‑148 (MQ=255)                                                                                                   gaggaAGCGTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAAc                                        <  2:1065240/148‑1 (MQ=255)                                                                                                   gaggaAGCCTTGCAGAAATTAAAATCTTTAGGTCTGAAAGATGATACGATTGAAAAAGTATATGATGATCGAGCGCCAAAAGGCTACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAAc                                        <  2:380045/148‑1 (MQ=255)                                                                                                   gaggaAGCCTTGCAAAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAAc                                        <  2:427080/148‑1 (MQ=255)                                                                                                    aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACt                                       >  2:412392/1‑148 (MQ=255)                                                                                                    aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACt                                       >  1:620830/1‑148 (MQ=255)                                                                                                    aggaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACt                                       >  2:662017/1‑148 (MQ=255)                                                                                                      gaAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGAATGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTAATGATTCTAATATTAAACaat                                     >  2:320760/1‑148 (MQ=255)                                                                                                          ccTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGa                                 >  1:90190/1‑148 (MQ=255)                                                                                                                 gAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAACTACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTa                          >  1:871995/1‑148 (MQ=255)                                                                                                                 gAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTa                          >  1:227778/1‑148 (MQ=255)                                                                                                                 gAAATTAAAATCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTa                          >  1:936147/1‑148 (MQ=255)                                                                                                                         aaTCGTTAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTaa                  <  2:218662/148‑1 (MQ=255)                                                                                                                              ttAGGTCTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATCCTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAg             <  1:68636/148‑1 (MQ=255)                                                                                                                              gtAGGTCTTAAAGATGAAAAGATTGAACAAGGAAATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAg             <  2:411475/147‑1 (MQ=255)                                                                                                                                   gtttAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTAt        <  1:413617/146‑1 (MQ=255)                                                                                                                                   gcTTAAAGATGATACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTAt        <  1:1001194/147‑1 (MQ=255)                                                                                                                                       aaaGAAGATAAGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGCGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTATGTag    <  2:52107/148‑1 (MQ=255)                                                                                                                                         aGATGATGCGGTTGAAAAAGTGTATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTGTAATATTAAACTATATGAATCTTTAGGCATTAAGGAAGTTTATGTagag  >  2:299228/1‑147 (MQ=255)                                                                                                                                              |                                                                                                                                               ATACTGGTGAACGTGTTGAACGTGGTGACAGTGTTGATGTTGTTATATCAAAGGGCCCTGAAAAGGTTAAAATGCCAAATGTCATTGGTTTACCTAAGGAGGAAGCCTTGCAGAAATTAAAATCGTTAGGTCTTAAAGATGTTACGATTGAAAAAGTATATAATAATCAAGCGCCAAAAGGATACATTGCAAATCAAAGTGTAACCGCAAATACTGAAATCGCTATTCATGATTCTAATATTAAACTATATGAATCTTTAGGCATTAAGCAAGTTTATGTAGAA  >  CP000730/1233363‑1233646

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