Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,868,981	T→C	9.5%	G68G (GGT→GGC)	adeP_1 →	Adenine permease AdeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,868,981	0	T	C	9.5%	113.9 / 4.9	42	G68G (GGT→GGC)	adeP_1	Adenine permease AdeP
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

GGATTTTCAAACTCAAGCAACACGGCAGCACCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCA  >  ACB122/4868848‑4869111                                                                                                                                      |                                                                                                                                   ggATTTTCAAACTCAAGCAACACGGCAGCACCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTg                                                                                                                                   >  4:276686/1‑135 (MQ=255)      ttCAAACTCAAGCAACACGGCAGCACCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCAACATGGCCTACATCATCTACGTCAACCCCAACATCATGGCCGACTCCGGCATCGACAATGGTGCCGCt                                                                                                                              >  2:184843/1‑135 (MQ=255)                         gcagcaCCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATcg                                                                                                           <  6:143005/135‑1 (MQ=255)                         gcagcaCCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATcg                                                                                                           <  7:241521/135‑1 (MQ=255)                         gcagcaCCGTCAGAACCGAAATGATCGCTGGGGGTAACACCTTCATCACGATGGCCTACATCATCTTCGTCAACCCCAACATCTTTGCCGACGCCGGCATAGACCATTGTGCCGCTTTTGTCGCCACCTGAATcg                                                                                                           <  7:208931/135‑1 (MQ=255)                          cagcaCCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATcg                                                                                                           <  5:215537/134‑1 (MQ=255)                          cagcaCCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATcg                                                                                                           >  6:215537/1‑134 (MQ=255)                                  tCAGAACCGAAATGATCGCTGGGGTGACCACCGTCATCCCCATGCCCTCCTTCATCTTCGTCACCCCCACCATCATGGCCGCCGCCGGCATCGCCCATGGTGGCGCTTAGGTCGCCAGCTGCATCCCCGCCgcgc                                                                                                  >  6:84825/1‑135 (MQ=255)                                     gAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTgg                                                                                               >  3:73413/1‑135 (MQ=255)                                        ccGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGc                                                                                             <  3:276686/134‑1 (MQ=255)                                              tGATCGCTGGGGTGACCCCCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGCCGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGctgcctgc                                                                                      >  4:202635/1‑135 (MQ=255)                                              tGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGctgccggc                                                                                      >  2:260194/1‑135 (MQ=255)                                                    cTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATgg                                                                                <  5:261458/135‑1 (MQ=255)                                                       gggTGACCACCTTCATCACCAAGGCCTACATAATCTTATTCAACCCCAACATCATGGCCGACGCCGGCATAGACCATGGTGCCCCTTTTGTCGCCACCTGCATCGCACCCGCGCTGGGCTGACTGCGTATGCGcc                                                                             >  8:1614/1‑135 (MQ=255)                                                               accTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGCGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGc                                                                      >  8:294348/1‑134 (MQ=255)                                                               accTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGCGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGc                                                                      >  8:290282/1‑134 (MQ=255)                                                                      tcaCCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTTTCGCCACCTGCATCCCCCCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACAGGc                                                              >  5:52777/1‑135 (MQ=255)                                                                       caCCATGGCATACATCATCTTCGGCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGcc                                                             <  5:19020/135‑1 (MQ=255)                                                                           aTGGCCTACATCATCTTCGGCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTg                                                         <  2:1420/135‑1 (MQ=255)                                                                              gCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCTGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGc                                                      >  7:142240/1‑135 (MQ=255)                                                                               ccTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGcc                                                     >  2:249916/1‑135 (MQ=255)                                                                                   catcatCTTCGCCAACCCCAACAGCATGGCCGACGCCGGCATCGACCATGGTGCCGCTGTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGgc                                                 <  1:100305/135‑1 (MQ=255)                                                                                    ctcatcTTCGTCAACCCCAAACTCATGGCCGCCGCCGGCATCGACCCTGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGgc                                                 <  5:223871/133‑1 (MQ=255)                                                                                           tCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCa                                         >  6:122383/1‑135 (MQ=255)                                                                                           tCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACAGGCCGGTGGGCCTGGCGCCGGGCa                                         >  4:78384/1‑135 (MQ=255)                                                                                           tCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCAGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCa                                         >  5:10143/1‑135 (MQ=255)                                                                                                     cAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCaa                               >  4:249238/1‑135 (MQ=255)                                                                                                          tcagGGCCGACGCCGGCATCGACCAAGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCct                          <  7:137507/135‑1 (MQ=255)                                                                                                           catGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCctt                         >  8:268774/1‑135 (MQ=255)                                                                                                           catGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCctt                         <  3:184261/135‑1 (MQ=255)                                                                                                                     gCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAGTGCCTTCGTTACCTac               >  3:36606/1‑135 (MQ=255)                                                                                                                     gCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTac               <  2:174426/135‑1 (MQ=255)                                                                                                                     gCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTac               <  7:268774/135‑1 (MQ=255)                                                                                                                             cGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGt       >  5:252273/1‑135 (MQ=255)                                                                                                                             cGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGt       <  6:52777/135‑1 (MQ=255)                                                                                                                                 cATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGc   >  8:162110/1‑135 (MQ=255)                                                                                                                                 cATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCAGCCTGCTGATGGGCCTGTCCGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCATTGCCTTCTTTACCTACACCGTGGTCGGc   >  8:252643/1‑135 (MQ=255)                                                                                                                                  aTGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCCTGGGGCTCACTGCCGTCTTTACCTACACCGTGGTCGGCa  >  4:63054/1‑135 (MQ=255)                                                                                                                                  aTGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCa  <  3:111970/135‑1 (MQ=255)                                                                                                                                  aTGGTGCCGCTTTGGGCGCCACCTGCATCGCCGCCGCGCGGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTGCTTTACCTACACCGTGGTCGGCa  <  1:260194/135‑1 (MQ=255)                                                                                                                                  aTGGCGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCa  <  7:290282/135‑1 (MQ=255)                                                                                                                                  aTGGCGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCa  <  7:294348/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   GGATTTTCAAACTCAAGCAACACGGCAGCACCGTCAGAACCGAAATGATCGCTGGGGTGACCACCTTCATCACCATGGCCTACATCATCTTCGTCAACCCCAACATCATGGCCGACGCCGGCATCGACCATGGTGCCGCTTTTGTCGCCACCTGCATCGCCGCCGCGCTGGGCTGCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCA  >  ACB122/4868848‑4869111

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