Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NODE_4	17,639	+A	intergenic (+1425/+305)	chaA → / ← kdsA	sodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NODE_4	17,638	1	.	A	95.8%	36.0 / ‑2.2	24	intergenic (+1424/+306)	chaA/kdsA	sodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (14/9);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17489‑17662                                                                                                                                                       |                         tCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑a                         <  1:944600/150‑1 (MQ=255) tCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑a                         <  2:3936228/150‑1 (MQ=255)  cGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aa                        >  1:1078749/1‑149 (MQ=255)   gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAGCAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:306855/150‑3 (MQ=255)   gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:2555976/150‑3 (MQ=255)   gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       >  1:3647451/1‑148 (MQ=255)   gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ac                        >  2:4025794/1‑150 (MQ=34)      gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  2:3486419/1‑145 (MQ=255)      gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  1:3280517/1‑145 (MQ=255)       taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   <  2:3195732/150‑7 (MQ=255)       taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   <  2:1703440/150‑7 (MQ=255)        aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  >  2:3527877/1‑143 (MQ=11)        aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  2:4035274/150‑8 (MQ=11)          ccACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCaa                 <  1:3457915/150‑1 (MQ=11)           cacaCCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAg                <  2:2728369/150‑1 (MQ=11)             cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              >  2:1358555/1‑150 (MQ=11)             cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              >  1:559896/1‑150 (MQ=255)              acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             <  1:1286041/150‑1 (MQ=11)               ccAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATa            <  2:1849679/150‑1 (MQ=11)                cAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAg           >  1:1446893/1‑150 (MQ=11)                 aacccGGCGTTATCCGGTGAGGCGCTATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          <  2:264901/147‑1 (MQ=11)                 aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  2:103813/1‑150 (MQ=11)                 aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  2:2062752/1‑150 (MQ=11)                  aGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGcc         >  1:3921703/1‑150 (MQ=11)                  aGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGcc         <  1:2585235/150‑1 (MQ=11)                   gCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCg        <  2:137627/150‑1 (MQ=11)                    ccGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGc       <  2:395294/150‑1 (MQ=11)                     cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTTAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCa      >  1:889873/1‑150 (MQ=11)                     cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGGTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑TACCAGAACAATTCAAGAATAACCGCa      >  1:1136621/1‑150 (MQ=11)                      ggCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑CACCAGAACAATTCAAGAATAGCCGCaa     >  1:617605/1‑150 (MQ=11)                       gCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaaa    >  2:2476387/1‑150 (MQ=11)                        cGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGGGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaaaa   >  1:2842371/1‑150 (MQ=11)                         gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:317786/1‑150 (MQ=11)                         gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:73119/1‑150 (MQ=11)                         gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  <  2:682607/150‑1 (MQ=11)                         gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAt  <  2:1042174/150‑1 (MQ=11)                         gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCCCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑TACCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:2447925/1‑150 (MQ=11)                                                                                                                                                       |                         TCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17489‑17662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑