Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,861,266	A→C	32.4%	A108A (GCA→GCC)	cbiA →	Cobyrinate a,c‑diamide synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,861,266	0	A	C	32.4%	52.9 / 2.6	42	A108A (GCA→GCC)	cbiA	Cobyrinate a,c‑diamide synthase
Reads supporting (aligned to +/- strand):  ref base A (18/9);  new base C (12/3);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTGGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTGGCCAATCGGGT  >  ACB122/1861140‑1861397                                                                                                                                |                                                                                                                                    ccTGTGGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACtt                                                                                                                             >  5:164265/1‑135 (MQ=255)   tgtgGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCCCTTTg                                                                                                                           >  4:44267/1‑135 (MQ=255)   tgtgGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCACTTTgg                                                                                                                          >  7:48446/1‑136 (MQ=255)   tgtgGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCACTTTg                                                                                                                           >  4:11477/1‑135 (MQ=255)   tgtgGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTgg                                                                                                                          >  8:313588/1‑136 (MQ=255)        atgatCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCACTTTGGCGTg                                                                                                                      >  6:58972/1‑135 (MQ=255)        atgatCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTg                                                                                                                      >  1:83107/1‑135 (MQ=255)        atgatCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGc                                                                                                                     >  7:330996/1‑136 (MQ=255)               gTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGc                                                                                                               >  1:213325/1‑135 (MQ=255)                        aaaGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACCCCACTTTGGCGTGCCGGTGCTGGCGGTgct                                                                                                     >  7:236331/1‑134 (MQ=255)                             cgtcgtCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCCTCAAGCGCCGACCTGGCCCCCCCCATTTGCCTGCCGGTGCGGGCCGTGATCGAc                                                                                                 >  3:49520/1‑135 (MQ=255)                             cggcggctgcggTGGGACGCTGCGGGCGAGGCTGACCTCACCCTCATCTAAGGGGTGAGGGGGCGGTT‑CGACGGTACCCCTTCAAGCGCCGACCGGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACg                                                                                                <  8:251036/130‑1 (MQ=255)                              gtcgtgTGCTGTGGGAGGCTGCGGGCCTGGCGGACCTGACCATCATCGAGGGGGTGAGGGGGGTGTT‑CCCCGTTACCCCGTAAAGCCCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACg                                                                                                <  6:201231/135‑1 (MQ=255)                                  tcTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGcac                                                                                            >  4:83522/1‑135 (MQ=255)                                   ctgctgTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAg                                                                                          <  8:330996/136‑1 (MQ=255)                                    tgctgTGGTAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGACCCCCGTAAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAg                                                                                          <  3:197829/135‑1 (MQ=255)                                       tgtgGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCACTTTGGCGTGCCGGGGCTGGCGGTGATCGACGGCCCAGc                                                                                         >  3:295201/1‑133 (MQ=255)                                              aGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACATTGGCGTGCCGGTGCTGGCGGTGATCGCCCGCACCGCCAAGGCCCa                                                                                >  6:318291/1‑135 (MQ=255)                                                 tgctgGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCCCTTTGGCGTGCCGGTGCTGGCGGTGATCGCCGGCACAGCCATGGCCCAGAcc                                                                            >  7:204513/1‑136 (MQ=255)                                                   ctgGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCCCCCCTTTTGCGGGGCGGGGCGGGCGGGGGTCGCCGGCCCCGCCCAGGCCCCGACCt                                                                           >  5:239365/1‑135 (MQ=255)                                                    gggCGAGGCTGACCTGAGCCTCCTCGACGGGGGGAGGGGGCTGGC‑CGCCGGCCCCCCGTCAAGCGCCGACCGGGCACGCCCCTTGGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCtt                                                                          <  1:242788/134‑1 (MQ=255)                                                         cggCTGACCTCAGCCTCATC‑CAGGGGGGATGGGGCGGTTCCTCCGGAACCCGGACATGCGCCTACGTGCCACGCCACTTTGGGGGGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGgcg                                                                     <  6:389487/134‑1 (MQ=255)                                                                       cTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTTGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGcc                                                      >  7:234906/1‑136 (MQ=255)                                                                       cTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGc                                                       >  5:152524/1‑135 (MQ=255)                                                                       cTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGc                                                       >  2:210630/1‑135 (MQ=255)                                                                         cATCGACGGGGGCAGGGGGCGGGG‑CCACGGAACCCGGTCAAGCGCCGCCCGGGCCCCCCACTTGGGCGGCCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTg                                                    <  7:56264/136‑1 (MQ=255)                                                                         cATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACATTGGCGTGCCGGTGCTGGCGGTGTTAGCCGGCCCACCCAAGGCCCAGACCTTCGGCCCCCCGGCCCCGGGccg                                                     >  5:317077/1‑134 (MQ=255)                                                                          agcGACGGGGGGAGGGGGCGGTT‑CCACGGACCCCCGCCAAGCGCCCACCGGGCCCGCCACTTTGGGGGGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTg                                                    <  3:182104/133‑1 (MQ=255)                                                                             gAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCCCTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGgcgc                                                >  8:238209/1‑136 (MQ=255)                                                                             gAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGCCCTTCGGCGCCCTGGCCCTGGGCCTGgcgc                                                >  7:77791/1‑136 (MQ=255)                                                                                       tGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAg                                                                                          <  6:202509/84‑1 (MQ=255)                                                                                       tGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAg                                                                                          >  5:202509/1‑84 (MQ=255)                                                                                               tt‑CGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCCCTTTGGCCTGCCGCGGCGGGGGGTGCTCGACGGCCCCCCCAAGGCCCCGCCCGTCGGCGCCCCGGCCCTGGGCCGGGCGCGGGTCCCGGCCgcgatg                               >  1:331927/1‑130 (MQ=255)                                                                                                  cGACGGTACCCCGTCAAGCGCCGACCTGGCCCGCCCCTTTGGCGTGGCGGTGCTGGCGGTGATCGACGGCCCAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGGACCAGGCCGCTCTGCCg                            >  8:70020/1‑136 (MQ=255)                                                                                                    aCGGTACCCCGTCAAGCGCCGACCTGGCCCGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGCTCTGCCGtt                          >  8:251428/1‑136 (MQ=255)                                                                                                                cAAGCGCCGACCTGGCCCGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCCCAGCCATGGCCCAGCCCTTCGGCGCCCGGGCCCGGGGCCTGGCGCGGGTCCAGGCCGCTCTGCCGTTCGGCGGCGGGc               >  3:302025/1‑135 (MQ=255)                                                                                                                cAAGCGCCGACCTGGCACGCCACATTGGCGTGCCGGTGGTGGCGGTGATCGCCGGCACCGCCCAGGCCCAGCCCTTCGGCGCCCGGGCCCGTGGCCTGGCGCGGGACCAGGCCGCTCGGCCCTTCGCCGGCGGGcg              >  7:358417/1‑135 (MQ=255)                                                                                                                  aGCGCCGCCCGGGCCCGCCACTTGGGCGGGCCGGTACTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTgg            <  8:264523/136‑1 (MQ=255)                                                                                                                  aGCGCCGACCTGGCCCGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCCCAGCCATGGCCCAGCCCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCCGGCCGCGCTGCCGGTCGCCGGCggggtg             >  2:24054/1‑130 (MQ=255)                                                                                                                   gcgcCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTgg            <  3:83522/135‑1 (MQ=255)                                                                                                                            tGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTGGCCAATCGGGt  >  7:140491/1‑136 (MQ=255)                                                                                                                             ggCACGCCACTTTGGCGGGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTGGCCAATCGGGt  <  6:239365/135‑1 (MQ=255)                                                                                                                                |                                                                                                                                    CCTGTGGATGATCGGTGCCGAGGAAAGCCGTCGTCTGCTGTGGGAAGCTGCTGGCGAGGCTGACCTGATCCTCATCGAAGGGGTGATGGGGCTGTT‑CGACGGTACCCCGTCAAGCGCCGACCTGGCACGCCACTTTGGCGTGCCGGTGCTGGCGGTGATCGACGGCACAGCCATGGCCCAGACCTTCGGCGCCCTGGCCCTGGGCCTGGCGCGGTACCAGGCCGATCTGCCGTTCGCCGGCGTGCTGGCCAATCGGGT  >  ACB122/1861140‑1861397

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