Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,805,508	C→T	50.0%	G535S (GGT→AGT)	pykA ←	pyruvate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,805,508	0	C	T	50.0%	‑6.4 / 21.5	20	G535S (GGT→AGT)	pykA	pyruvate kinase
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (6/4);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTTACCTTGAGCAGCATCAATCGTAACTAACATGTTATTGCTTATGTTTTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGTTTCAGCAACTAACGTAGTACC  >  CP000730/1805383‑1805646                                                                                                                              |                                                                                                                                           tttACCTTGAGCAGCATCAATCGTAACTAACATGTTATTGCATATGTATTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTTTTAACGAACAATTTCACTTGGTAGTGTAATACTATGttcttgtatgat                                                                                                                             >  2:25188/1‑134 (MQ=255)     ccTTGAGCAGCATCAATCGTAACTAACATGTTATTGCTTATGTTTTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAg                                                                                                                         >  1:15618/1‑141 (MQ=255)              gcaTCAATCGTAACTAACATGTTATTGCTTATGTTTTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTACTTAAGCCTAACGCt                                                                                                                <  2:284330/141‑1 (MQ=255)                                         cTTATGTTTTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCa                                                                                     <  2:161440/141‑1 (MQ=255)                                                          tttttCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAAc                                                                    <  1:89699/141‑1 (MQ=255)                                                                tACACCTACAACTGTTTGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAAc                                                              <  1:107271/141‑1 (MQ=255)                                                                    ccTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTg                                                          <  2:82311/141‑1 (MQ=255)                                                                                 gAATACCTTTTTCTAAACCAACAATTGGACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGACTAAAGCTTTTTCTACATAAGGTGCAAAGATGTCATCGATGAAGATAGTAACGATAAATTTGTGAGGTTAAACtt                                             >  2:193944/1‑141 (MQ=255)                                                                                     aCCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTAcc                                         <  1:127947/141‑1 (MQ=255)                                                                                      ccTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACACAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCt                                        <  2:244432/141‑1 (MQ=255)                                                                                                aaCCGACAATTGCACTTGGTGATGTAATACTAGTTTCTTCAGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTAACCTTCTAAATCtt                              >  1:204979/1‑141 (MQ=255)                                                                                                aaCCAACAATTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCtt                              >  1:242938/1‑141 (MQ=255)                                                                                                   caacaaTTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCCAAAGCTCTTTCTACATAAGGTACAAACGTTTCATCCATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTGCCTTCTAAATCTTTaa                           <  1:170387/141‑1 (MQ=255)                                                                                                    aacaaTTGCACTTGGTGGTGTAATACTACACTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAAc                          >  2:262872/1‑141 (MQ=255)                                                                                                    aacaaTTGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTATATCTTTaaa                          >  2:95051/1‑140 (MQ=255)                                                                                                         ttGCACTTGGTGATGTAATACTATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGttt                     >  2:281087/1‑141 (MQ=255)                                                                                                         ttGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTAAGATAAATCTTTACCTTCTAAAGCTTTAACAGttt                     >  1:227778/1‑141 (MQ=255)                                                                                                           gCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGTTTCa                   <  1:33243/141‑1 (MQ=255)                                                                                                                         aaTACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGTTTCAGCAACTAACgtagt     <  1:20649/141‑1 (MQ=255)                                                                                                                            aCCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGTTTCAGCAACTAACGTAGTAcc  >  1:50767/1‑141 (MQ=255)                                                                                                                              |                                                                                                                                           TTTACCTTGAGCAGCATCAATCGTAACTAACATGTTATTGCTTATGTTTTTAACAGCTTTTTCTACACCTACAACTGTTGGAATACCTTTTTCTAAACCAACAATTGCACTTGGTGATGTAATACCATTTTCTTCTGTAATTAAGCCTAAAGCTTTTTCTACATAAGGTACAAACGTTTCATCGATTGAGTTAGTAACGATAACTTTGTCAGATAAATCTTTACCTTCTAAATCTTTAACAGTTTCAGCAACTAACGTAGTACC  >  CP000730/1805383‑1805646

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: ‑