Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	2,708,478	T→G	32.0%	C46G (TGC→GGC)	CJ019_02558 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	2,708,478	0	T	G	32.0%	50.2 / 2.6	42	C46G (TGC→GGC)	CJ019_02558	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (8/5);  total (13/29)
Fisher's exact test for biased strand distribution p-value = 2.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCGTATCCTTGCCGGCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCACGCGGCTGTTCTAC  >  ACB122/2708348‑2708609                                                                                                                                   |                                                                                                                                    gCATCGTATCCTTGCCGGCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACCCCCGGGCCATTCAACCGATCTTCACCGAGAAGTGCCTggccggccc                                                                                                                                 >  4:16928/1‑134 (MQ=255)          ccTTGCCGGCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCCACCGATCTTCACCGAGAAGTGCGGGGGCGGCCCCGCCGGCac                                                                                                                        >  5:195777/1‑134 (MQ=255)             tGCCGGCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGGCGAGCCCGGCGATTTCCTACACCCGGGCCATTCCACCGATCTTCCCCGAGAAGTGCGGGGGCGGCCCCCCCGGCAccgc                                                                                                                     >  6:212055/1‑134 (MQ=255)                 ggCGCCTTCGCTCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGGCATTCAACCGATCTTCACCGAGAAGTGCGGGGGCGGCCCCGCCGGGAACgccgcc                                                                                                                 >  3:218599/1‑135 (MQ=255)                 ggCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTTCACCCGGGCCATTCAACCGATCTTCACCGCGAAGTGCGGGGCCCGGCACGCCGGCAacgccccc                                                                                                                 >  2:207172/1‑130 (MQ=255)                 ggCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACCCCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCGGCCCCGCCGGGCACgccgcc                                                                                                                 >  8:139038/1‑135 (MQ=255)                 ggCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGCTCTTCACCGAGAAGTGCGTGGCCAGCCCCCCCGGCCACgccgcc                                                                                                                 >  8:167723/1‑135 (MQ=255)                             cTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCATACACCCGGGACATTCAACCGATCTTAACCGAGAAGGGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTc                                                                                                     <  3:140410/135‑1 (MQ=255)                                   atcaGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTg                                                                                               <  4:249088/135‑1 (MQ=255)                                    tcaGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCt                                                                                                <  3:36342/133‑1 (MQ=255)                                      aGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCGGGCACGCCGGCACCCCCGCCGCCGGCCCGCTGAAGCggg                                                                                              >  7:177232/1‑133 (MQ=255)                                             cgGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCCCCGAGAAGTGCGTTGGCTGCCCCCCCAGCCACGCCGCCGCCAGCCAGCTCCAGCAGGGAAGcccc                                                                                      >  3:236477/1‑133 (MQ=255)                                              ggTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGGGGCGGGGCACCCCGGGAAAGAACCCCCCGGCCAGCCGCAGCTGGGAAGcccc                                                                                      >  5:238982/1‑132 (MQ=255)                                              ggTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGGGGCCGGCCACCCCCGGAACCCCCCCGCCAGCCCCCCCAACCGGGAAAGcccagg                                                                                    >  7:15595/1‑132 (MQ=255)                                                tATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAg                                                                                  <  8:190940/135‑1 (MQ=255)                                                 aTTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATGTTCACCGAGAAGTGCGGGGCCTGCCACGCCTGCAACGACGCCGCCTGTCAGCTCAAGCTGGAAAGCCCTGAAgg                                                                                 <  5:295652/135‑1 (MQ=255)                                                   ccGGGCCGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGCACATTCACCCGATCTTCACCGAGACGTGCGGGGCCTGCCACGCCGGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGgcg                                                                               <  5:227831/134‑1 (MQ=255)                                                     gggCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCCTCGCCGTCCCCCCCGGCCACCCCCCCCCCGGCCACCACAAACCGGCACGGCCGGAAGGcgcg                                                                             >  3:34871/1‑135 (MQ=255)                                                             cccccAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCGGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCgg                                                                     >  7:219800/1‑135 (MQ=255)                                                                 cAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGcc                                                                 <  2:242874/135‑1 (MQ=255)                                                                  aGTCGAGCCCGGCTATTTCCAACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGCCGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCa                                                                <  3:290626/135‑1 (MQ=255)                                                                  aGTCGAGCCCGGCTATATCCTACACCCGGGACATTCAACCGATCATCACCGAGAAGTGCGTGGCCGGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCa                                                                <  1:292131/135‑1 (MQ=255)                                                                      gAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCaa                                                            <  3:76206/135‑1 (MQ=255)                                                                      gAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGGGCGTGGCCGGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCaa                                                            <  7:119859/135‑1 (MQ=255)                                                                        gCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAgg                                                          <  1:183639/135‑1 (MQ=255)                                                                          ccGTCTATTTCCTCCCCCCGGGACATTCACCCGAGCTTCCCCGAGATGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTc                                                        <  3:101169/135‑1 (MQ=255)                                                                          ccGGCTATTTCCTACCCCCGGGACATTCAACCGCTCTTCCCCGCGACGTGCGGGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTc                                                        <  4:211921/135‑1 (MQ=255)                                                                          ccGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTc                                                        <  8:219800/135‑1 (MQ=255)                                                                                      tccacCCGGCCCATACACCCGAGCTGCCCCGCGCCGGGCGGGGCCGCCCCCGCCGGCAACGACGCCGCCTGCCAGCTCAAGCGGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTAcc                                              <  5:11006/131‑1 (MQ=255)                                                                                         acCCGGGCCATTCAACCGATCTTCACCGAGAAGGGCGGGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCCAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGc                                         <  5:154052/135‑1 (MQ=255)                                                                                         acCCGGGACATTCACCCGACCTTCACCGAGACGTGCGTGGCCGCCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGCCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGc                                         <  1:207172/135‑1 (MQ=255)                                                                                         acCCGGCACATGCAACCGATCTTCCCCGCGACGGCCGGCGCCCGCCACGCCTGCAACGCCGCCGCCTGCCAGCTCAAGCGGGCAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGc                                         <  1:250870/135‑1 (MQ=255)                                                                                           ccGGGAAATTCAACCGATCTTCACCCAGAAGGGCGGGGCCGGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGa                                       <  7:167723/135‑1 (MQ=255)                                                                                                aCATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCCCCCCCGCAACGACGCCGCCGGACCGCGCAAGCTGGAAAGCCCTGAACGCGCGGTACGCGGGGCCAGCAAGGGCCCCGGGTCCCAGGGCGCGCGGa                                  >  2:235219/1‑135 (MQ=255)                                                                                                   ttCAACCGCTCTTCACCGAGAAGTGCGGGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGACGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCa                               <  2:286407/135‑1 (MQ=255)                                                                                                                 aCCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTccc                     <  2:73619/131‑1 (MQ=255)                                                                                                                    gagaAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCAcg              <  5:156535/135‑1 (MQ=255)                                                                                                                    gagaAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCAcg              <  5:110519/135‑1 (MQ=255)                                                                                                                    gagaAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCAcg              <  7:139038/135‑1 (MQ=255)                                                                                                                    gagaAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCAcg              <  4:218599/135‑1 (MQ=255)                                                                                                                    gagaAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCAcg              <  3:16928/135‑1 (MQ=255)                                                                                                                           gCGTGGCCTGCCACGCCTGCCACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCACGCGGCTGt       <  6:238982/135‑1 (MQ=255)                                                                                                                                cccTGCCACGCGTGCAACGACGCCGCCTGCCAGATCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCGGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCACGCGGCTGTTCTac  <  1:49685/134‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    GCATCGTATCCTTGCCGGCGCCTTCGCCCTGCTCATCAGCGGCGCGGTATTCGGGCAGGCCCCCCAGTCGAGCCCGGCTATTTCCTACACCCGGGACATTCAACCGATCTTCACCGAGAAGTGCGTGGCCTGCCACGCCTGCAACGACGCCGCCTGCCAGCTCAAGCTGGAAAGCCCTGAAGGCGCGGTACGCGGGGCCAGCAAGGTCCCGGTGTACCAGGGCGAGCGGAGCAAGGCAGTCCCCACCACGCGGCTGTTCTAC  >  ACB122/2708348‑2708609

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