Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,487,100	A→T	8.5%	intergenic (+116/+14)	yhdZ → / ← degS	amino acid ABC transporter ATP‑binding protein/serine endoprotease DegS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,487,100	0	A	T	8.5%	119.4 / 3.3	47	intergenic (+116/+14)	yhdZ/degS	amino acid ABC transporter ATP‑binding protein/serine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base A (27/16);  new base T (2/2);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCTG  >  NC_002947/1486960‑1487245                                                                                                                                             |                                                                                                                                                  gttCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGcc                                                                                                                                            <  1:46761/148‑1 (MQ=255)    ccTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCCCTACCCCAGGTTTAGCGCTGGGGTTGGAGCCTGCGTAATTACTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCggg                                                                                                                                        >  3:178971/1‑149 (MQ=39)        aGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGTTGCGGCCCttcgcgg                                                                                                                                    <  7:100055/149‑7 (MQ=38)                ccTCCATCGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCCCGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGCGCCTGCGTAATTCCTGGGGGCGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTctct                                                                                                                            <  4:18222/149‑1 (MQ=39)                 cTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTctct                                                                                                                            >  1:90200/1‑148 (MQ=255)                              gTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGTTTAGGGCTGGCGTTGGAGCCTGCGTTATTCCTGGGGGCGCGGGCGGGCCCGCGAACGGCCGCCACGCGGCCCTGGGTTTATTTGTCTTCCTGGCTGGCAGc                                                                                                              >  1:315230/1‑149 (MQ=38)                                   ggCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTTGCGCTGGCGTTGGAGCCTGCGTAATTGCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCGTGGGTTATTTGTCTTGCTGGCTGGCAGCCGgcg                                                                                                         >  1:200595/1‑149 (MQ=37)                                                            aCGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTc                                                                                >  7:339359/1‑149 (MQ=255)                                                              gCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCgg                                                                              >  3:301410/1‑149 (MQ=255)                                                              gCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCgg                                                                              >  4:109752/1‑149 (MQ=255)                                                                ccGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAg                                                                                         >  7:21598/1‑136 (MQ=255)                                                                ccGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAg                                                                                         <  8:21598/136‑1 (MQ=255)                                                                      ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCg                                                                                                      >  6:384647/1‑117 (MQ=255)                                                                      ccACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCg                                                                                                      <  5:384647/117‑1 (MQ=255)                                                                       cacaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                     >  2:200549/1‑149 (MQ=255)                                                                        acaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt                                                                    >  7:188883/1‑149 (MQ=255)                                                                        acaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt                                                                    >  3:46592/1‑149 (MQ=255)                                                                        acaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCt                                                                     >  5:231782/1‑148 (MQ=255)                                                                        acaGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTAGGCGATCAGCtt                                                                    >  4:26740/1‑149 (MQ=255)                                                                         caGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCtt                                                                    >  2:368154/1‑148 (MQ=255)                                                                             tttAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAAct                                                               <  4:301410/149‑1 (MQ=255)                                                                                aGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGACGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAActgct                                                            >  5:89631/1‑149 (MQ=255)                                                                                  cggtGGCTTTGTAGCCGGCGTAATTCTGGTGGGTGGGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACtgctgc                                                          <  3:58145/146‑1 (MQ=255)                                                                                   gcTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCg                                                                                            <  8:261936/114‑1 (MQ=255)                                                                                   gcTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCg                                                                                            >  7:261936/1‑114 (MQ=255)                                                                                      ggCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGt                                                      >  8:39686/1‑149 (MQ=255)                                                                                      ggCTTTGGAGCCTGCGTAATTCCTGTGGGAGCCGGCGTGCCCGCGAAGGGCCGCAACGCCGCCCTGGGTTTATTTGTCTTCCTTGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCCCCTCGGCGATGAGCATCAACTGCTGCCCGt                                                      >  1:29320/1‑149 (MQ=255)                                                                                       gCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCa                                                                                                                >  4:347941/1‑90 (MQ=255)                                                                                       gCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCa                                                                                                                <  3:347941/90‑1 (MQ=255)                                                                                            gccgCCGGCGTAATGCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTTGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTgcgc                                                 <  2:200595/145‑1 (MQ=255)                                                                                               gCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCa                                             <  1:368154/149‑1 (MQ=255)                                                                                                           ccTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGtt                                  <  2:315230/148‑1 (MQ=255)                                                                                                               tgGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTAtctt                             >  4:321440/1‑149 (MQ=255)                                                                                                                gggAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTAtcttc                            >  3:168184/1‑149 (MQ=255)                                                                                                                      gggCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTg                      <  4:178971/149‑1 (MQ=255)                                                                                                                      gggCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTg                      <  2:90200/149‑1 (MQ=255)                                                                                                                        gCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTggg                    <  1:200549/149‑1 (MQ=255)                                                                                                                            gCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCg                                                                               >  8:5539/1‑86 (MQ=255)                                                                                                                                acccaGGGCCGCGTTGCGGCCCTTGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTgg                     >  4:102383/5‑140 (MQ=255)                                                                                                                                acccaGGGCCGCGTTGCGGCCCTTGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTgg                     <  3:102383/136‑1 (MQ=255)                                                                                                                                acccaGGGCCGCGTTGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATcc            >  5:209098/5‑149 (MQ=255)                                                                                                                                  gAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTAGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCcgc          <  3:26740/149‑1 (MQ=255)                                                                                                                                        ccGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTAcc    >  8:128926/1‑149 (MQ=255)                                                                                                                                         cGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGc                                                                                        <  2:154055/64‑1 (MQ=255)                                                                                                                                         cGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGc                                                                                        >  1:154055/1‑64 (MQ=255)                                                                                                                                         cGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCt   >  1:224577/1‑149 (MQ=255)                                                                                                                                          gCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGCTCCGCGCTTCCTg  >  6:366341/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                  GTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCTGGCTTTGGAGCCTGCGTAATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCAACGCGGCCCTGGGTTTATTTCTCTTCCTGGCTGGCAGCCGGCGCCGGTGGCGGCCGCAGCCCCACCTCGGCGATCAGCTTCAACTGCTGCCCGTTGCGCATCACCTCGATGGTTATCTTCTCGTTGGGCTTGATCCGCGCTACCTG  >  NC_002947/1486960‑1487245

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