Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	260,050	G→A	10.6%	R238H (CGC→CAC)	rssB_1 →	Regulator of RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	260,050	0	G	A	10.6%	107.6 / 2.9	43	R238H (CGC→CAC)	rssB_1	Regulator of RpoS
Reads supporting (aligned to +/- strand):  ref base G (31/6);  new base A (2/2);  total (33/10)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGCTGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAGGGGCGTTGTTATTGCG  >  ACB122/259920‑260181                                                                                                                                   |                                                                                                                                    gcggTGCTGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAagcgcagc                                                                                                                                 <  8:268171/135‑1 (MQ=255) gcggTGCTGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAagcgcagc                                                                                                                                 >  6:1650/1‑135 (MQ=255) gcggTGCTGATCGCTGCCGCCAATAGTGCGGAATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGAACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAcgcgcagc                                                                                                                                 <  4:393718/135‑1 (MQ=255) gcggTGCTGATCGCGGCCGCCAATAGTGCGGCTTTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAagcgcagc                                                                                                                                 <  4:374460/135‑1 (MQ=255)   ggTGCTGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAagcgcagcgc                                                                                                                               <  8:249839/135‑1 (MQ=255)        tGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGcc                                                                                                                          >  1:277021/1‑135 (MQ=255)        tGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGcc                                                                                                                          >  3:13994/1‑135 (MQ=255)        tGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCACAGCGCACGcc                                                                                                                          >  5:14803/1‑135 (MQ=255)        tGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCACAGCGCACGcc                                                                                                                          >  5:137296/1‑135 (MQ=255)             gcggccgccAATAGTGCGGCATTGCCCCGTGAAGCCCCGGGGCAGACGGTGCTGCAGGCCCGCAACACGCTTGGCAGTACCCAAAGCATGAACCTGATCCGGGGCATGCCCCTCAAGCGCCGCGCACGCCTCAGt                                                                                                                     <  7:94919/132‑1 (MQ=255)             caggccgccAATAGGGCGGCATTGCACCGTGACGCCCCGGTGCAGACGTGGCTGCCGGCCCTGAACAAGCTTGGCAGAACCCAAAGCATGAACGTGATCGTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGt                                                                                                                     <  6:359329/132‑1 (MQ=255)               tGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGa                                                                                                                   >  7:16917/1‑135 (MQ=255)                              ggCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGcctg                                                                                                    >  5:359329/1‑132 (MQ=255)                              ggCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGa                                                                                                    >  4:167173/1‑135 (MQ=255)                                cATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTGCCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCTCGCCTCAGTGACCCGCTGTTGGCGGAGc                                                                                                  >  7:204954/1‑135 (MQ=255)                                cATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGc                                                                                                  >  5:382131/1‑135 (MQ=255)                                cATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGc                                                                                                  >  3:147821/1‑135 (MQ=255)                                        cGTGAAGCCCCGTTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCCTGAACCTGATCCTGGGCATGACCCTCAAGAGCAGCGCACGCCTCAGTGACACGCTGTTGGCGGAGCACGCGAcc                                                                                          >  7:173410/1‑135 (MQ=255)                                        cGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGAcc                                                                                          >  5:117982/1‑135 (MQ=255)                                        cGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGAcc                                                                                          >  6:328156/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCCGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCACGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTa                                                                                     >  4:236616/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGCTCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTa                                                                                     >  7:297600/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGCCCCGTTa                                                                                     >  8:455139/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTa                                                                                     >  3:256523/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTa                                                                                     >  3:78325/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTa                                                                                     >  5:250195/1‑135 (MQ=255)                                             aGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGCCAGAACCCAAAGCATGAACATGATCCTGGGCATGAACCTCAAGCGCAGCGCACGCCTCACTGACCCGCTGTTGGCGGAGCAAGCGACCCGTTa                                                                                     >  7:266665/1‑135 (MQ=255)                                                  cGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTggg                                                                                >  7:386887/1‑135 (MQ=255)                                                  cGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTggg                                                                                >  3:403221/1‑135 (MQ=255)                                                              ggctgcCGGCCCTCAACACGCTTGGCAGTACCCAAAGCATCACCCTGACCCTGGCAATGCCCCTCACGCCCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGc                                                                    <  7:97925/134‑1 (MQ=255)                                                                  gcAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCTACCCGTTACTGGGACCTGTCACTGcaccc                                                                >  3:279777/1‑133 (MQ=255)                                                                  gcAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCacac                                                                >  1:350093/1‑135 (MQ=255)                                                                      gCCCTGAACAAGTTGGGCTGGACCCAAAGCATGAACCTGAGCGTGGGCATGCCCCTCACGCGCCGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGcc                                                            <  2:303117/135‑1 (MQ=255)                                                                        ccTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGa                                                          >  7:194493/1‑135 (MQ=255)                                                                         cTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCATCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGCCCCGTTACTGGGACCTCTCACTGCACACTGCCGa                                                          >  8:133079/1‑134 (MQ=255)                                                                                    tGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCAc                                              >  7:466875/1‑135 (MQ=255)                                                                                          cacccaacccaTCAACCTCACCCGGGGCATCACCCTCACCCCCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGgc                                        <  3:402133/134‑1 (MQ=255)                                                                                              cAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTCTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGgcgc                                      >  2:422466/1‑133 (MQ=255)                                                                                               aaaGCATCAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGgcgcg                                     >  7:11857/1‑133 (MQ=255)                                                                                                    aTGAACCTGATCCTGGGCTGGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTc                              <  8:261574/135‑1 (MQ=255)                                                                                                       aaCCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGaa                           >  6:25442/1‑135 (MQ=255)                                                                                                                cggggCAGGACCATCAACCACAGCGCCCGCCTCCGTTACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAg                  <  6:14803/133‑1 (MQ=255)                                                                                                                cTGGTCATGACCCTCAAGCGCCGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAg                  <  4:127594/135‑1 (MQ=255)                                                                                                                cTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGACGCACCTTGGCGCGCATGCTCGAACTGGACGAg                  >  1:168226/1‑135 (MQ=255)                                                                                                                 ggggCATCACCCTCAAGCACCGCGCACGCCGCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAg                  <  6:137296/133‑1 (MQ=255)                                                                                                                           ccTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCTTGCTCGAACTGGACGAGGGGCGTTGTTa       >  7:238664/1‑135 (MQ=255)                                                                                                                            cTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAGGGGCGTTGTTAt      >  3:446431/1‑135 (MQ=255)                                                                                                                             tCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCAGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAGGGGCGTTGTTAtt     <  8:16917/135‑1 (MQ=255)                                                                                                                                agcgcagcgcACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAGGGGCGTTGTTATTgcg  >  3:272731/1‑135 (MQ=255)                                                                                                                                   |                                                                                                                                    GCGGTGCTGATCGCTGCCGCCAATAGTGCGGCATTGCACCGTGAAGCCCCGGTGCAGACGTTGCTGCAGGCCCTGAACAAGCTTGGCAGTACCCAAAGCATGAACCTGATCCTGGGCATGACCCTCAAGCGCAGCGCACGCCTCAGTGACCCGCTGTTGGCGGAGCACGCGACCCGTTACTGGGACCTGTCACTGCACACTGCCGAATATGGCCGCACCTTGGCGCGCATGCTCGAACTGGACGAGGGGCGTTGTTATTGCG  >  ACB122/259920‑260181

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