Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,793,849	A→G	12.9%	E51G (GAA→GGA)	rbsB →	Ribose import binding protein RbsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,793,849	0	A	G	12.9%	80.4 / 4.0	31	E51G (GAA→GGA)	rbsB	Ribose import binding protein RbsB
Reads supporting (aligned to +/- strand):  ref base A (12/15);  new base G (2/2);  total (14/17)
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

CTGGCCCTCGCTGTCATTTCCTCGCTTAGCCTCGCCTTGCCATTGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAGCGCT  >  ACB143/2793725‑2793984                                                                                                                             |                                                                                                                                        ctgGCCCTCGCTGTCATTTCCTCGCTTAGCCTCGCCTTGCCATTGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGgcg                                                                                                                                 <  6:133597/133‑1 (MQ=255)       cTCGCTGTCATTTCCTCGCTTAGCCTCGCCATGCCAATGTCAGCCGCCCCCGCTGAAGACAAGCCCAACCTCGCCCTTGGCATGAAGTCGCTGGCCAACCACTTCATACCCACCAACGAAAACCGCGCCCAGGaa                                                                                                                         >  4:215987/1‑134 (MQ=255)           cTGTCATTTCCTCGCTTAGCCTCGCCTTGCCATTCTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCACCGAGTTCTTCCGCACCATGGAAGACGGCGACAAGGACTAc                                                                                                                      >  8:297680/1‑134 (MQ=255)                      tagCTTAGCCTAGCCTTGCCATTGTCTGCCCCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTATTACGCACCAAGGAAGACGGCGCCAAGGACTACCAGAAAACCCAtt                                                                                                         >  8:33485/3‑135 (MQ=255)                        gCTTAGCCTCGCCTTGCCATTGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGTCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGc                                                                                                        <  7:61574/135‑1 (MQ=255)                                           ttGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCg                                                                                    <  2:245271/136‑1 (MQ=255)                                           ttGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCg                                                                                    <  4:164357/136‑1 (MQ=255)                                             gTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAATACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGc                                                                                   <  1:234752/135‑1 (MQ=255)                                                     cccACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAACCCCCATGCCGACGAGTTCGAGCTGATCGCCCACGGCAt                                                                          >  7:312158/1‑136 (MQ=255)                                                     cccACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCAt                                                                          >  7:352226/1‑136 (MQ=255)                                                           cTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGaa                                                                    >  7:183487/1‑136 (MQ=255)                                                               aGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACga                                                                 <  7:267312/135‑1 (MQ=255)                                                                      cccAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGaca                                                         >  1:245583/1‑136 (MQ=255)                                                                      cccAAGGGCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGaca                                                         <  5:125966/136‑1 (MQ=255)                                                                       ccAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGacac                                                        >  8:320719/1‑136 (MQ=255)                                                                                  cTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGGAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAg                                                 >  6:5328/1‑132 (MQ=255)                                                                                  cTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGGAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAg                                                 >  6:360/1‑132 (MQ=255)                                                                                      tCATGACGTCGGTGGCCAACGAGTTCTTCCGCACCAGGGACGACGGCGCCAAGGACTACCAGAAAACCCATCCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATcc                                          <  3:215987/135‑1 (MQ=255)                                                                                      tCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATcc                                          <  8:183487/135‑1 (MQ=255)                                                                                        aTGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCTAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCa                                       >  6:318407/1‑136 (MQ=255)                                                                                             gTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCAtcgtc                                  <  8:91582/136‑1 (MQ=255)                                                                                                gCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAg                               >  7:244388/1‑136 (MQ=255)                                                                                                gCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAg                               <  1:405000/136‑1 (MQ=255)                                                                                                    gCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAg                               >  6:367593/1‑132 (MQ=255)                                                                                                      cacccaGTTCTTCCGCCCCATGCAAGACGGCGCCAAGGACTACCATAAAACCCATGCCTACTAGTTCGATCTGATCGCCAACTGCATCAAGACCGATACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATg                         <  2:386200/131‑1 (MQ=255)                                                                                                      caccGAGTTCTTCCGCACCACCGACGACGGCGCCAAGGACTACCAGAAATCCCATGCAGCCGATTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTCGCGAGCAGATCCGCATCGTCGAGCAGATg                         <  2:210945/133‑1 (MQ=255)                                                                                                            gTTCTTCCGCCCCATGGGAGACGGCGCCAAGGACTACCAGAAACCCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAAt                   <  5:5328/136‑1 (MQ=255)                                                                                                            gTTCTTCCGCACCATGGGAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAAt                   <  5:360/136‑1 (MQ=255)                                                                                                               cttcCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGc                 <  2:289961/135‑1 (MQ=255)                                                                                                                    gCACCATGGAAGACGGCGCCAAGGACTACCAGAAACCCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGgcgc           <  6:22106/136‑1 (MQ=255)                                                                                                                            gAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAgcgc   >  1:90229/1‑136 (MQ=255)                                                                                                                             aaGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCTCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAGCGCt  >  7:294549/1‑136 (MQ=255)                                                                                                                             aaGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAGCGCt  >  6:326718/1‑136 (MQ=255)                                                                                                                             aaGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAGCGCt  >  5:232599/1‑136 (MQ=255)                                                                                                                             |                                                                                                                                        CTGGCCCTCGCTGTCATTTCCTCGCTTAGCCTCGCCTTGCCATTGTCTGCCGCCCACGCTGAAGACAAGCCCAAGGTCGCCCTTGTCATGAAGTCGCTGGCCAACGAGTTCTTCCGCACCATGGAAGACGGCGCCAAGGACTACCAGAAAACCCATGCCGACGAGTTCGAGCTGATCGCCAACGGCATCAAGAACGAGACCGACACTGGCGAGCAGATCCGCATCGTCGAGCAGATGGTCAATGCCGGCGCCAAAGCGCT  >  ACB143/2793725‑2793984

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