Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	5,515,237	C→T	G157S (GGC→AGC)	CJ019_05187 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	5,515,237	0	C	T	96.1%	172.2 / ‑1.8	52	G157S (GGC→AGC)	CJ019_05187	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (22/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGCGCGAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCCGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGCAGGCAGGTAGCCCTCGA  >  ACB143/5515107‑5515365                                                                                                                                   |                                                                                                                                 tGCGCGAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGAt                                                                                                                               >  7:78420/1‑134 (MQ=255) tGCGCGAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGAt                                                                                                                               >  5:265784/1‑134 (MQ=255)  gcgcGAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCtgatga                                                                                                                             <  5:336146/135‑1 (MQ=255)      gAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATc                                                                                                                        >  8:39685/1‑136 (MQ=255)      gAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATc                                                                                                                        <  7:47290/136‑1 (MQ=255)      gAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATc                                                                                                                        >  7:197266/1‑136 (MQ=255)      gAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATc                                                                                                                        >  2:291995/1‑136 (MQ=255)       aaCACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCg                                                                                                                       <  1:52166/136‑1 (MQ=255)         cacGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTc                                                                                                                     <  1:149947/136‑1 (MQ=255)              tGTTGCGAAACGCTTTCACCAGCTGGCTAAGACGTTCGAACAGCGCCGGCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGGTGGCACCGAGGATCATGCTGATGAGCATCGTCACGaa                                                                                                                <  1:189611/136‑1 (MQ=255)                ttGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAAc                                                                                                               <  7:39685/135‑1 (MQ=255)                gtGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAAc                                                                                                               <  7:319419/134‑1 (MQ=255)                  gCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACAt                                                                                                             <  4:312994/135‑1 (MQ=255)                     aaaCGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTgc                                                                                                         >  1:125699/1‑136 (MQ=255)                     aaaCGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTgc                                                                                                         >  6:217525/1‑136 (MQ=255)                         gCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCa                                                                                                     >  6:9685/1‑136 (MQ=255)                                   tcagGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGcc                                                                                           <  5:146714/135‑1 (MQ=255)                                   gcagGCTAAGACGTTCGAACATCGCCGCCGCCAGGGGCTTGCGACGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGcc                                                                                           <  4:252378/136‑1 (MQ=255)                                   gcagGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGcc                                                                                           <  1:362629/136‑1 (MQ=255)                                       gCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCCGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCccac                                                                                        >  1:47193/1‑135 (MQ=255)                                       gCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCCGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCccac                                                                                        >  2:90894/1‑135 (MQ=255)                                          aaGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAg                                                                                                                            <  8:263121/96‑1 (MQ=255)                                          aaGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAg                                                                                                                            >  7:263121/1‑96 (MQ=255)                                           aGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCt                                                                                   <  6:273649/136‑1 (MQ=255)                                            gACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTg                                                                                  <  6:321640/136‑1 (MQ=255)                                                    aaCAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCCGCTGCAAATCGc                                                                          >  4:106965/1‑136 (MQ=255)                                                    aaCAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGc                                                                          >  7:288083/1‑136 (MQ=255)                                                    aaCAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGc                                                                          >  7:286567/1‑136 (MQ=255)                                                    aaCAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGc                                                                          >  3:37809/1‑136 (MQ=255)                                                       aGCGCCGCCGCCAGGGGCTTGCGGCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCa                                                                       <  7:124474/136‑1 (MQ=255)                                                       aGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCAGCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCa                                                                       <  4:110255/136‑1 (MQ=255)                                                                 ccAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCtt                                                             >  6:168623/1‑136 (MQ=255)                                                                            cgcCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGc                                                  <  2:10168/136‑1 (MQ=255)                                                                             gcCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGATCATCGTCACGATCATGTGCGCCATCCCCTTGCCCACCAGCTGCAACTCGATCAGGTGGCTCTTGATCCAGTCGcc                                                 >  8:317335/1‑136 (MQ=255)                                                                             gcCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGcc                                                 >  5:42453/1‑136 (MQ=255)                                                                             gcCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGcc                                                 >  3:116167/1‑136 (MQ=255)                                                                               cGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGa                                               <  3:362685/136‑1 (MQ=255)                                                                               cGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGa                                               <  5:160070/136‑1 (MQ=255)                                                                               cGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGa                                               <  5:9685/136‑1 (MQ=255)                                                                               cGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGa                                               <  5:191245/136‑1 (MQ=255)                                                                                          gggAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCGGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCtt                                     <  5:320385/135‑1 (MQ=255)                                                                                              aTGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAAct                                <  1:344014/136‑1 (MQ=255)                                                                                                      cAAAGCGATGATGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCgg                          <  3:159428/134‑1 (MQ=255)                                                                                                            gatgatGGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTggca                  <  5:168623/136‑1 (MQ=255)                                                                                                                 tgGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTc                                                                        >  3:330043/1‑77 (MQ=255)                                                                                                                 tgGCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTc                                                                        <  4:330043/77‑1 (MQ=255)                                                                                                                   gCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGcaggcagg            >  7:5030/1‑135 (MQ=255)                                                                                                                   gCACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGcaggcagg            >  8:111386/1‑135 (MQ=255)                                                                                                                    cACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGCAGGCAGGTa          <  7:242026/136‑1 (MQ=255)                                                                                                                    cACCGAGGATCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGCAGGCAGGTa          <  1:346703/136‑1 (MQ=255)                                                                                                                            aTCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGCAGGCAGGTAGCCCTCGa  >  6:211856/1‑136 (MQ=255)                                                                                                                             tCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATg                                             <  1:336235/92‑1 (MQ=255)                                                                                                                             tCATGCTGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATg                                             >  2:336235/1‑92 (MQ=255)                                                                                                                                   |                                                                                                                                 TGCGCGAACACGATGTTGCGAAACGCTTTCACCAGCAGGCTAAGACGTTCGAACAGCGCCGCCGCCAGGGGCTTGCGCCGGGAAATATCGGGAATGCGCTGCAAAGCGATGATGGCACCGAGGATCATGCCGATGAGCATCGTCACGAACATGTGCGCCATACCCTTGCCCACCAGCTGCAAATCGCTCAGGTGGCTCTTGATCCAGTCGCCGATGGCCACCTTGAACTCTGCGGCACTGGCAGGCAGGTAGCCCTCGA  >  ACB143/5515107‑5515365

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