Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,500,213	A→C	44.4%	L61R (CTG→CGG)	cysG_1 ←	Siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,500,213	0	A	C	44.4%	18.0 / 4.1	32	L61R (CTG→CGG)	cysG_1	Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (13/4);  new base C (12/3);  total (25/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGAC  >  ACB122/4500080‑4500343                                                                                                                                      |                                                                                                                                   cgcATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTGAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCCTCCGGATAGCCACGCACCAg                                                                                                                                   >  3:58695/1‑135 (MQ=255)        accacATTCACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCCCCTCAc                                                                                                                           >  8:232010/1‑136 (MQ=255)                  ccGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCt                                                                                                                   <  4:291019/134‑1 (MQ=255)                    ggCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCAGCCCCTCACCGCCACCTTcc                                                                                                                >  3:345858/1‑135 (MQ=255)                    ggCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGAACACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATATCCACCCACCAGCACCTCACCGCCAccgtac                                                                                                                >  3:81530/1‑131 (MQ=255)                       aGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCggg                                                                                                            >  8:317735/1‑136 (MQ=255)                         gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCCGCACCCGCCCAGGGCGGCTGCCGGGTAGCCCCCCACCCCCCCCGCCCCGCCCCCCTCCCcggg                                                                                                           >  5:95610/1‑134 (MQ=255)                          cTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGcc                                                                                                          <  4:246285/135‑1 (MQ=255)                                gCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCCGTCGCCGCGATCACCAGCCGGCCACCGACCAGGTCGGCTGCCTGATAGCCCCGCCCCCGCACCTCCCCGCCCCCTTACCCGGCCCAcgcg                                                                                                    >  5:176209/1‑135 (MQ=255)                                      gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCCCCAGCCGGCAACCGCCCAGGTCGGCGGCCTGATAGCCCCGCACCAGCCCCTCCCCGGCACCGTCCCGGGGCCACGCGGCCacc                                                                                              >  2:312043/1‑133 (MQ=255)                                      gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCTGATAGCCCCGCCCCCGCCCCTCCCCGCCCCCGTCCCGGGCCCACCCGGccc                                                                                                >  2:76033/1‑132 (MQ=255)                                      gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCCCCCGCACCTCCCCGCCACCTTCCCGGGCCAACGCGGCCAAc                                                                                              >  4:123682/1‑135 (MQ=255)                                       cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCCGTCGCCGCGCTCACCAGCCCGCAACCGCCCAGGTCGGCTGCCGGCAAGCCACGCACCCGCACCACACCGCCCCCTTCCCGGGGCCACGCCGCCAACt                                                                                             >  1:315206/1‑135 (MQ=255)                                       cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCCGCCCCGCCCCGCCCCCGTCCCGGGCCCACGCCGCCCAca                                                                                             >  4:230720/1‑134 (MQ=255)                                        aTCGGCCGACACCTGGGCATTTAGCCCAGGGTCGTCAGTCGCCGCGCTCCCCAGCCGGCAACCGCCCAGGTCGGCTGCCGGATAGCCACGCCCCCGCACCTCCCCGCCCCCTTCCCTGGCCCACCCGGCCCACTg                                                                                            >  2:100153/1‑135 (MQ=255)                                        aTCGGCCGACACCTGGGCATTCCGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCCCCGCCCAGGTCGGCGGCCGGCTAGCCCCGCCCCCGCCCCGCACCGCCCCCATCCCGGGCCCACGCGGccccctg                                                                                            >  5:164496/1‑130 (MQ=255)                                             ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGGCTGATAGGCACGCCCCCGCACCGCACCCGCCCCATCCCCGGCCACCGCCGCCCACAGGCCggc                                                                                      >  8:363951/1‑134 (MQ=255)                                             ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGGCGGGTAGCCCCGGCCCCGCCCCGCACCGCCCCCGTCCCGGGCCCCCCCGGCCCCCAGGCCgg                                                                                       >  5:309436/1‑134 (MQ=255)                                                  acCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCTGATAGCCACGCCCCCGCCCCGCACCGCCCCCTTCCCGGGCCAACGCGGCCCACTGGCCGTCGCCg                                                                                  >  4:293694/1‑135 (MQ=255)                                                          aTTCAGCCCAGGGTCGTCCGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCGGATAGCCACGCACCAGCCCCTCCCCGCCACCTGCCCGGGCCCACCCGGCCCACAGGCCGTCGACCTCCGGcccc                                                                         >  7:226622/1‑136 (MQ=255)                                                          aTTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCAGCACCTCCCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGGCGACGTCCGgcgc                                                                          >  6:58500/1‑135 (MQ=255)                                                           ttCAGCCCAGGGTCGTCTGGCGCGGCGATCCCCCGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACCCACCAGAACCTCACCGCCCCCTTCCCGGGCCAACGCGGCCAACTGCCCGTCGACGTCCGGCGcc                                                                         <  6:384289/135‑1 (MQ=255)                                                              aGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGGCACCCACCAGCACCTCCCCCCCACCTTCCCGGGCCACCCCGGCCAACAGGCCGTCGACGGCCGGCGcccccc                                                                     >  8:103595/1‑135 (MQ=255)                                                                ccccGGGTCGCCCGGCGCGGCTACCACCCCCCGCAACCCGACCCGGTCGGCGGCCTGATAGCCACGCACCCGCACCTCACCGCCACCTACCCGGGCCACCGCGCCCAACTGGCCGTCGACGTCCGGCGCcaccac                                                                    <  5:129487/135‑1 (MQ=255)                                                                 ccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCCGCCCCTCACCGCCCCCGTCCCGGGCCCACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCCcg                                                                   >  1:85088/1‑135 (MQ=255)                                                                  cAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCCCCCcgcg                                                                 >  8:192950/1‑136 (MQ=255)                                                                             gTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCAACGCCGCCCACAGGCCGTCGACGTCCGGCGCCCCCACGCGCACCAcgcccc                                                      >  7:187822/1‑134 (MQ=255)                                                                                  cgcgATCACCAGCCGGCAACCGCCCAGGTCGGCTGGCTGATAGCCACGCCCCAGCACCTCACCGCCACCGTCCCGGGCCAACGCCGCCAACTGGCCCTCGACGTCCCGCGCCACCACGCGCAGCACgccccccg                                                   >  2:234877/1‑132 (MQ=255)                                                                                    cgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCAACCCGGCCAACTGGCCGTCCACGTCCGGCGCCGCCACGCCCAGCACGCCGCCGGcct                                                >  2:327146/1‑133 (MQ=255)                                                                                             gccgccACCCGACCCGGGCGGCTGCCTGATACCCACGCACCACCCCCTCACCGCCACCTCCCCGGGCCACCGCGCCCACCGGGCCGGCGACGTCCGGCGCCCCCACGCGCAGCACGCCGCCGGCATCGGCcagca                                       <  3:266309/135‑1 (MQ=255)                                                                                                                ggcggcCGGATACCCACCCACCCCCACCTCACCCCCCCCTCCCCGGCCCACCGCGCCCAACTGGCGGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGc                    <  6:333912/135‑1 (MQ=255)                                                                                                                                  cACCCGAACCTCACCCCCCCGTTCCCGGGCCACCGCGGCCAAGTGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACcgac  <  1:234877/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   CGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGAC  >  ACB122/4500080‑4500343

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