Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,955,371	G→A	23.5%	P128L (CCG→CTG)	atoC_3 ←	Regulatory protein AtoC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,955,371	0	G	A	23.5%	66.4 / 15.0	34	P128L (CCG→CTG)	atoC_3	Regulatory protein AtoC
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base A (5/3);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGAATGTAGCGCGCCAGCACTTCCTTGCCAGTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCACGCTCGACGGCGGCAT  >  ACB143/4955245‑4955506                                                                                                                               |                                                                                                                                        tGAATGTAGCGCGCCAGCACTTCCTTGCCAGTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCcttctt                                                                                                                                <  6:56829/136‑1 (MQ=255)             gcCAGCACTTCCTTGCCAGTACCGGACTCGCCGGAAATCAGTCCGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCacac                                                                                                                    <  2:148027/136‑1 (MQ=255)                          tGCCAGTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCATGCCACCGGCCCTTCATCGCCGGTCCCACGCCCGGCGGCATg                                                                                                       >  8:260969/1‑136 (MQ=255)                           gCCAGTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATg                                                                                                       <  2:66249/135‑1 (MQ=255)                               gTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCt                                                                                                  <  8:172787/136‑1 (MQ=255)                               gTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCt                                                                                                  <  8:165472/136‑1 (MQ=255)                               gTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCt                                                                                                  <  3:220383/136‑1 (MQ=255)                               gTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCt                                                                                                  <  8:183742/136‑1 (MQ=255)                               gTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCt                                                                                                  >  4:102464/1‑136 (MQ=255)                                  ccGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCa                                                                                               >  2:42685/1‑136 (MQ=255)                                          gCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGtt                                                                                        >  2:134067/1‑135 (MQ=255)                                           ccGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGa                                                                                      <  3:269421/136‑1 (MQ=255)                                                      tACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCGCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTg                                                                           <  1:289437/136‑1 (MQ=255)                                                         ggTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGc                                                                        >  7:166883/1‑136 (MQ=255)                                                          gTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCt                                                                       <  3:160333/136‑1 (MQ=255)                                                          gTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCt                                                                       <  5:28480/136‑1 (MQ=255)                                                              aGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCATGGGCTCGaa                                                                   >  5:237663/1‑136 (MQ=255)                                                               gTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAAc                                                                  <  1:140490/136‑1 (MQ=255)                                                                tCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACg                                                                 <  8:116803/136‑1 (MQ=255)                                                                 cACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCTGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACgg                                                                >  4:2203/1‑136 (MQ=255)                                                                       gcgcCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGAc                                                          >  2:171308/1‑136 (MQ=255)                                                                        cgcCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGAc                                                          >  7:217709/1‑135 (MQ=255)                                                                                gggCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTc                                                 >  2:220833/1‑136 (MQ=255)                                                                                 ggCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGACCGGTTTGACCAGGTAGTcg                                                >  5:313565/1‑136 (MQ=255)                                                                                     gccAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGCTTGAGCAGCGCCTTGTGCTCGAACGGTTTGACCAGGTAGTcggcgg                                            >  4:260567/1‑136 (MQ=255)                                                                                       cAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGgcggcg                                          <  6:239028/136‑1 (MQ=255)                                                                                             cAGCACCTGGCGACTGGCCGGCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTgg                                    <  4:4430/136‑1 (MQ=255)                                                                                             cAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTgg                                    <  7:37571/136‑1 (MQ=255)                                                                                                  aCTGGCGACT‑GCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATg                              >  3:233715/1‑136 (MQ=255)                                                                                                      gCGACTGGCCGGCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATgg                             <  1:7768/134‑1 (MQ=255)                                                                                                        gACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGACGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTc                         >  4:282515/1‑136 (MQ=255)                                                                                                                ggCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCAc                 >  2:54394/1‑136 (MQ=255)                                                                                                                ggCTCGCAGGCCACCAGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTGGGGCTCGAACGGTTTGACCAGGTAGTCCGCGGCGCCTTGGCGCCTGGCCTCGACCGCAc                 >  4:211618/1‑136 (MQ=255)                                                                                                                       aGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCACGCTCGAc          <  7:101433/136‑1 (MQ=255)                                                                                                                               ggCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCACGCTCGACGGCGGCAt  <  8:292508/136‑1 (MQ=255)                                                                                                                               ggCCCTTCTACGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCACGCTCGACGGCGGCAt  <  3:28075/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                        TGAATGTAGCGCGCCAGCACTTCCTTGCCAGTACCGGACTCGCCGGAAATCAGTACGGTCGAGTCACTGCGCGCCACACGGGCCGCCAGCTCCAGCAACTGGCGACTGGCCGGCTCGCAGGCCACCGGCCCTTCTTCGCCGGTCACACGCCCGGCGGCATGCCGCTCCACGAGGTTGAGCAGCGCCTTGGGCTCGAACGGTTTGACCAGGTAGTCGGCGGCGCCTTGGCGCATGGCCTCGACCGCACGCTCGACGGCGGCAT  >  ACB143/4955245‑4955506

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