Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,532,409	T→G	intergenic (‑758/+28)	sgcX ← / ← yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,532,409	0	T	G	96.2%	82.3 / ‑3.7	26	intergenic (‑758/+28)	sgcX/yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (14/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

GATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGC  >  NC_000913/4532174‑4532651                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   gATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGACATa                                                                                                                                                                                                                                       <  2:63139/249‑1 (MQ=255)    aCCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCggg                                                                                                                                                                                                                                  <  2:166767/251‑1 (MQ=255)               ccAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATAACAGGCACTGAGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCATGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTCATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGACAGAAATACCGGGAACCTgttgtt                                                                                                                                                                                                                       >  2:238222/1‑251 (MQ=255)                 aGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTAc                                                                                                                                                                                                                     <  1:1243/251‑1 (MQ=255)                     tttGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCggg                                                                                                                                                                                                                 >  1:200757/1‑251 (MQ=255)                                    gCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAac                                                                                                                                                                                                   <  1:21592/250‑1 (MQ=255)                                                      acaATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTTTGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCAtt                                                                                                                                                                                <  1:61509/251‑1 (MQ=255)                                                        aaTATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCACAGCATTcc                                                                                                                                                                              >  2:84413/1‑251 (MQ=255)                                                                                              aTCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAg                                                                                                                                        <  1:236263/251‑1 (MQ=255)                                                                                                              ccTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGGTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcg                                                                                                                        <  2:25431/251‑1 (MQ=255)                                                                                                                                 ttgtGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGCCGCGACGCTTAGCTCGGCCCTGACTTCGGCGCCGAAGTCATCGTc                                                                                                       >  2:47566/1‑249 (MQ=255)                                                                                                                                                   cTTATGGTGACCTCTTGATTAATAAACATACTGACTATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGCACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTCGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGCCTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCcagccagc                                                                                   >  2:50227/1‑251 (MQ=255)                                                                                                                                                             cATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGACCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGATTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGt                                                                         >  1:238953/1‑251 (MQ=255)                                                                                                                                                                 ttGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGc                                                                     >  1:183573/1‑251 (MQ=255)                                                                                                                                                                                   gAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATa                                                    <  1:174030/250‑1 (MQ=255)                                                                                                                                                                                             tgtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAg                                          <  1:46359/250‑1 (MQ=255)                                                                                                                                                                                              gtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                                         <  2:215705/250‑1 (MQ=255)                                                                                                                                                                                                    gTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTg                                  >  1:25797/1‑251 (MQ=255)                                                                                                                                                                                                     ggATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTgg                                 <  2:77755/250‑1 (MQ=255)                                                                                                                                                                                                        tACATCAGATCAACGCAAAGGTAAAGCCTCGATATTGCCCGCCAGAACTACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGCAGCCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                              >  2:251712/1‑251 (MQ=255)                                                                                                                                                                                                            tCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc                           >  1:126192/1‑250 (MQ=255)                                                                                                                                                                                                              aGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt                        >  2:29790/1‑251 (MQ=255)                                                                                                                                                                                                                              aaaGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCcagccagc                                                                                   >  1:160958/1‑176 (MQ=255)                                                                                                                                                                                                                              aaaGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCcagccagc                                                                                   <  2:160958/176‑1 (MQ=255)                                                                                                                                                                                                                                 gCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCATAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAt     >  2:64612/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATacc  >  2:243644/1‑249 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   GATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGC  >  NC_000913/4532174‑4532651

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