Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,514,291	A→G	18.2%	K336K (AAA→AAG)	curA →	curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,514,291	0	A	G	18.2%	45.2 / 4.0	22	K336K (AAA→AAG)	curA	curcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base G (2/2);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTT  >  CP009273/1514049‑1514491                                                                                                                                                                                                                                                   |                                                                                                                                                                                                         aCAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAg                                                                                                                                                                          >  1:66686/1‑275 (MQ=255) aCAGAGCTACCACCCAGTCCGGATCGTTCACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACNGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtggt                                                                                                                                                                                                     >  1:499371/1‑248 (MQ=255)                                                    cTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAGGTGGTGATCCGCGTGGCGGGTGATGATTAAAGGTTCAACGGCGGCGTAAGCCGCCATAAAAATGCAGGCAAAANGCTTaacat                                                                                                                       >  2:273439/1‑271 (MQ=255)                                                     tGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGAt                                                                                                                                                                                                  >  1:263626/1‑199 (MQ=255)                                                         aaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCcgc                                                                                                                                                                                              <  2:263621/199‑1 (MQ=255)                                                                 ttCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAGGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaat                                                                                                             <  1:273444/272‑1 (MQ=255)                                                                                                          cGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaacta                                                                                                       >  1:271481/1‑237 (MQ=255)                                                                                                              tCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACca                                                                                                   <  2:271476/237‑1 (MQ=255)                                                                                                                                        cAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCtgctgc                                                                                  >  2:719855/1‑228 (MQ=255)                                                                                                                                            gggTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                                                                              <  1:719865/228‑1 (MQ=255)                                                                                                                                                       gATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAGAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGatta                                                  >  1:627738/1‑245 (MQ=255)                                                                                                                                                           aaaTCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCa                                              <  2:627729/245‑1 (MQ=255)                                                                                                                                                                        cgAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTATAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                                                                                  >  2:403502/1‑180 (MQ=255)                                                                                                                                                                            gagaTTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                                                                              <  1:403509/177‑1 (MQ=255)                                                                                                                                                                            gaaATTACTGTCGTGTTAGTGAATTCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGGTCCGCATGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGtt  <  1:397716/272‑1 (MQ=255)                                                                                                                                                                                                    gcgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCAGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTCATTAAATAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGatta                                                  >  1:176168/1‑200 (MQ=255)                                                                                                                                                                                                                  tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaacta                                                                                                       >  1:257256/1‑133 (MQ=255)                                                                                                                                                                                                                      ggCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACca                                                                                                   <  2:257251/133‑1 (MQ=255)                                                                                                                                                                                                                                      aaaaCTTCGGCAAGGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGat                                                    >  2:486281/1‑164 (MQ=255)                                                                                                                                                                                                                                          cTTCGGCAAGGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattatt                                                <  1:486288/164‑1 (MQ=255)                                                                                                                                                                                                                                            tCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                                                                              >  2:372317/1‑132 (MQ=255)                                                                                                                                                                                                                                                cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaaca                                                                          <  1:372323/132‑1 (MQ=255)                                                                                                                                                                                                                                                   aGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTaa                                                                                                                                                         >  1:92135/1‑50 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                         ACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTT  >  CP009273/1514049‑1514491

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