Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,035	C→T	G188D (GGC→GAC)	yjcS ←	linear primary‑alkylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,035	0	C	T	95.7%	63.0 / NA	23	G188D (GGC→GAC)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/14);  minor base G (0/1);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGAC  >  NC_000913/4305803‑4306271                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             cgtcgGTGCAATAATCCTCGGGTCCCCCTTTGCCAGCGGCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACTGTAAGAGTAGAGCGCCCGGCGGCGCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACTTCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGgtgt                                                                                                                                                                                                                            <  1:178151/251‑1 (MQ=255)                   cgtgtcgccgcttcacagcgccccGCCCAGGCCATTACCCCCATGGCCTTGCGCGTTGTCCGGCAGTAACAGCCCGTAAGGGTAGAGCGCACGGCGGCTCTGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGTGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATaa                                                                                                                                                                                                          <  1:148244/228‑1 (MQ=255)                   cGGGTCGCCCTTTGCCAGCGTCACGCCAAGGCCATTACCCATATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGGCACCATAGTGGTCGGTGGGTCTGTGAGTGTAGAt                                                                                                                                                                                                            <  1:203224/249‑1 (MQ=255)                        cGCCCGTTGCCAGCGTCACGCCAAGACCATTCCCCACATTGCCTGGCGCGTTGGGGGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGca                                                                                                                                                                                                     <  1:14864/251‑1 (MQ=255)                            cgtgcccAGCGCCACGCAAGGGCCTTTACCCCAATGGCCTTGCGCGTGGTGCGGCGGTACCAGACCGTAAGAGTAGAGCGCAGGGCGGCTATTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAAGCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaacaa                                                                                                                                                                                                 <  2:168718/246‑1 (MQ=255)                             gTTGCAAGCGTCACGCCAAGGCAATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTACCAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGGGATAATGCCTTTCACGTCACCTTAGTGGGAGGTGTGGCTGTGTGTGTCGATAACGGCAACAAt                                                                                                                                                                                                >  2:241271/1‑251 (MQ=255)                                     cGTCACGCCAAGGCCATTACCCCCTTTGCCTTGCGCGGTGTGCGGCATCAACAGCCCGAAAGAGAAGAGCGCACGGCGGCTCTTGATGTTACCCGCCAGCACGTTTTCGCTTATGGCTTCGTCCATAAAGCCTGCAGGGTCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTTAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCGGTGAGTGTAGATAACGGCAACAATCGGttttt                                                                                                                                                                                        <  1:21768/251‑1 (MQ=255)                                                                    ggcgcGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCCCGGCGGCCCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTg                                                                                                                                                          <  1:246389/249‑1 (MQ=255)                                                                         gTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTAATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCtt                                                                                                                                                    <  1:32112/251‑1 (MQ=255)                                                                           ggtggggcggTAACAAACCTTAAGAGTGGGGCGCAGGGGGTCTCATGATGTTGCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTTCCGGATTTAACATCGGCTTCAGAGATAATGCCTGTCACGTCACCATAGTGGTCGGTGGGGCTGTGAGTGTAGATAACGGCAACAAGCGGTTTTTGCGGACGATGCTGGAAGGAAAGGTCAAGTGCGGCtttt                                                                                                                                                  <  2:187933/242‑1 (MQ=255)                                                                                 gCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGGGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGc                                                                                                                                             >  1:102705/1‑250 (MQ=255)                                                                                   aGTAACCGACCGTACGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGAGGACGGTGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAgg                                                                                                                                          >  2:259172/1‑251 (MQ=255)                                                                                        cAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                                                     >  1:150026/1‑251 (MQ=255)                                                                                           aCCGTAAGAGTAGAGCGCACGACGGCTCATGATGTTACCCGCCAGCAAGTTTTCGCTGATGGCTTCGTCCATAAAACCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATAGGCTTCAGAGATAATGCCTTTCACGTCCCCATCGTGGTCGGGGTGGCTGTGAGTGTAGCTAACGGCAACAATTGGTTTTTGCGGCCGATGCTGGAATTAAAGGTCAAGAGCGGCTTTTGCGGccgggtgggcg                                                                                                                                   >  2:24650/1‑240 (MQ=255)                                                                                                     tAGAGCCCACCGCGACTAATGATGTTACCAGCCAGCACGTTTCCGCTGATGGCTTCGACCATAAAGACTGCAAGGGCAATCAACTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGACGGTGTGGCTGTGAGTGTAGATAACGGCGAAAATCGGTTTTTGCGGACGATGCTGGCAGTAAAGGTTCAGTGCAGCTTTTGCGACAggcggagcccccagcagggc                                                                                                                        >  2:124254/1‑236 (MQ=255)                                                                                                     aagagCGCACGGCGGCTCATGATGTTGCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGTATTTAACATCGGCTTCAGAGATAATGCCGTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                                                        <  1:111299/250‑1 (MQ=255)                                                                                                           gcccggcggCTCTGGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTGTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAAc                                                                                                                  <  1:259172/248‑1 (MQ=255)                                                                                                                                                cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTa                                                                             <  1:77926/251‑1 (MQ=255)                                                                                                                                                               ccATAAAGCCTGCAGGGGCAATCACCTGAACGTTGCCGGATTTAACATTGGCTTCAGAGATAATGCCGTTCACGTCACCATAGTGGTTGGTGGGGTTGTGAGTGTAGAGAACGGCAACAAGCGGTTTTTGCGGACGATGCTGGAAGTAAAGTTCAAGTGCGGCTTTTGCGGCAGTCGGCGTCACCAGCTGGTCGATAACAATAATGCCTTTCACGCCCTCAACGAACGTAATGTTAGAGATATCCTGGcc                                                               <  2:180124/250‑1 (MQ=255)                                                                                                                                                                               ggCAATCGCCTGAACTTTGCCTGATTTAACATCGGCTTCAGAGATAATTCCTTTCACGTCACCATAGTGGTCGGTTTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCTGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAttt                                              <  1:6799/251‑1 (MQ=255)                                                                                                                                                                               ggCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGGAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAttt                                              <  1:128660/251‑1 (MQ=255)                                                                                                                                                                                            aCTTTGCCGGATTTAACATCGGCTTCAGCGATAATGCCTTTCACGTCACCATAGAGCTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGCCGAAAACAAAAATGCCTTTCTCGCCCTCAACGAACGTAATGTTCGAGATATCCTGGCCGCCCCCCTGATACATTTTATCGGTGACttt                                 >  2:255833/1‑251 (MQ=255)                                                                                                                                                                                                                           tAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGCTACGTTTTATCGGTGACTTTGAACAGGCCCGAAATACCGTTGATTTGCGAc  >  1:204653/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             CGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGAC  >  NC_000913/4305803‑4306271

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