Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,645,929	T→C	11.4%	F270L (TTC→CTC)	aceA →	Isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,645,929	0	T	C	11.4%	96.1 / 5.0	35	F270L (TTC→CTC)	aceA	Isocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (13/18);  new base C (2/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CCTGGCCGCTGACGTGTCGGGTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGA  >  ACB143/4645796‑4646053                                                                                                                                      |                                                                                                                             ccTGGCCGCTGACGTGTCGGGTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGATCGTGGTTGGCGAGCGCACCCGCGAAGGCTTc                                                                                                                            <  5:312473/136‑1 (MQ=255)  cTGGCCGCTGACGTGTCGGGTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCt                                                                                                                           <  3:25059/136‑1 (MQ=255)            aCGTGTCGGGTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCcg                                                                                                                 <  5:238653/136‑1 (MQ=255)                    ggTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGcgcccgcc                                                                                                         >  8:113044/1‑136 (MQ=255)                    ggTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGc                                                                                                          >  4:163827/1‑135 (MQ=255)                       gtACCGACCATCATCCTGGCGCGTACCGATGCCACCGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAACGTCCGCGCCGGCCTc                                                                                                       >  3:108385/1‑135 (MQ=255)                         aCCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGAc                                                                                                    >  5:124836/1‑136 (MQ=255)                              ccATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGc                                                                                               <  2:84991/136‑1 (MQ=255)                                atcatcCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGCCTGCGCCCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCa                                                                                             >  2:111805/1‑136 (MQ=255)                                  catcCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCAtt                                                                                           >  2:9506/1‑136 (MQ=255)                                   atcCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTg                                                                                          <  7:64822/136‑1 (MQ=255)                                         gcgcGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCcgc                                                                                     <  6:295201/135‑1 (MQ=255)                                         gcgcGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCcgc                                                                                     <  5:281233/135‑1 (MQ=255)                                         gcgcGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCcgc                                                                                     <  2:297291/135‑1 (MQ=255)                                         gcgcGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCcgc                                                                                     <  3:252176/135‑1 (MQ=255)                                         gcgcGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCcgc                                                                                     <  7:123267/135‑1 (MQ=255)                                          cgcgTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCgg                                                                                   >  6:294790/1‑136 (MQ=255)                                          cgcgTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCCTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCgg                                                                                   >  5:167696/1‑136 (MQ=255)                                          cgcgTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCCTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCgg                                                                                   >  5:163933/1‑136 (MQ=255)                                            cgTACCGATGCCAACGCCGCCGACCTGCTGCCCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCCCCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCggcc                                                                                 >  2:174842/1‑136 (MQ=255)                                            cgTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCAATGCCCGCggcc                                                                                 >  3:131428/1‑136 (MQ=255)                                                  gATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCggcctggcct                                                                           <  6:47627/136‑1 (MQ=255)                                                  gATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCggcctggcct                                                                           <  7:36891/136‑1 (MQ=255)                                                      ccAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGc                                                                       <  3:163827/136‑1 (MQ=255)                                                      ccAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGc                                                                       >  3:135125/1‑136 (MQ=255)                                                                  aCCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCt                                                                                                                           <  4:268080/72‑1 (MQ=255)                                                                  aCCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCt                                                                                                                           >  3:268080/1‑72 (MQ=255)                                                                                  cTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGa                                            >  8:121778/1‑135 (MQ=255)                                                                                         ccGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCa                                    <  6:156927/136‑1 (MQ=255)                                                                                         ccGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCa                                    <  3:87345/136‑1 (MQ=255)                                                                                                   aGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCt                          <  2:227234/136‑1 (MQ=255)                                                                                                       gTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGAc                      >  2:285160/1‑136 (MQ=255)                                                                                                                          aCCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCg   <  3:309432/136‑1 (MQ=255)                                                                                                                           cccGCGAAGGCCTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGa  <  6:167696/136‑1 (MQ=255)                                                                                                                           cccGCGAAGGCCTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGa  <  6:163933/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                             CCTGGCCGCTGACGTGTCGGGTGTACCGACCATCATCCTGGCGCGTACCGATGCCAACGCCGCCGACCTGCTGACCAGCGACTGCGACCCGTACGACCAGCCGTTCGTGGTTGGCGAGCGCACCCGCGAAGGCTTCTACAAGGTCCGCGCCGGCCTCGACCAGGCCATTGCCCGCGGCCTGGCCTACGCCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGA  >  ACB143/4645796‑4646053

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