Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,971,703	G→A	E148E (GAG→GAA)	wecC →	UDP‑N‑acetyl‑D‑mannosamine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,971,703	0	G	A	93.5%	78.6 / ‑0.8	32	E148E (GAG→GAA)	wecC	UDP‑N‑acetyl‑D‑mannosamine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (16/13);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGT  >  NC_000913/3971481‑3971947                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGCTCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGcc                                                                                                                                                                                                                          >  1:448274/1‑251 (MQ=255)  gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGccc                                                                                                                                                                                                                         <  1:131847/251‑1 (MQ=255)    aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGCGCGATCATGCGCCAGCTATGACCTCCGTCGCATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGCGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGACCAGGCGGACGTCAACATTGCTTACTGCCCgg                                                                                                                                                                                                                       >  1:249001/1‑251 (MQ=255)                          gTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATAGGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTGCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCGTACTGCCCGGAACGCGTGTTACCAGGACAgg                                                                                                                                                                                                  <  2:146511/250‑1 (MQ=255)                                  gCCGTTTAAGAGCGATCATGAGCCAGATCTGACCTCCGTTGCATCGGCTGCTCGCTCCATTGCGCCAGTGCTGCAAAAAAGCGAGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCCGCATGGTTAGCAGAGATGCGTCCGGCTCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGAGTTACAAGGACCGGTAATGGTCg                                                                                                                                                                                         >  2:374382/1‑251 (MQ=255)                                       ttAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTg                                                                                                                                                                                    >  1:32439/1‑251 (MQ=255)                                        tAAGGGCTATCATGAGCCAGATATGACCTACGTTGAATCTGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTg                                                                                                                                                                                    <  1:298969/250‑1 (MQ=255)                                                   aTGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGAAATGGTCGAGCTGATTAAAAACGAt                                                                                                                                                                        >  1:294931/1‑251 (MQ=255)                                                                      cGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGGGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGCTCTCACTTTACAGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGGTTACTGCCCGGAACGCGTGTTCCCAGGCCAGGTAATGGTCGGGCTGATTAAAAACGATCGCGTGATtggtggaagg                                                                                                                                                      >  2:436697/1‑246 (MQ=255)                                                                                   gctcgctcCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTc                                                                                                                                        <  1:248819/251‑1 (MQ=255)                                                                                     tcgctcCATTACGCCAGTGCTGAAAAAAGGCGCGCTGGTGCTCCTTGAATCCACCTCGCCGGTGGGGGCAACCGAGAAGCTGGCAGAATGGTTCGCAAAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTTCCCGGAACGCGGGTTACCCTGACAGGAAATGTGCGCGCTGCTTAAAAACGCTCGCGTGATTGGAAGTATGACGCAGGATTGGTCgg                                                                                                                                      >  2:314659/1‑251 (MQ=255)                                                                                        ctcCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGccc                                                                                                                                   >  1:241514/1‑251 (MQ=255)                                                                                                cgcCAGTGCTGAAAACAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGc                                                                                                                           >  1:158080/1‑251 (MQ=255)                                                                                                     gTGCTGAAAAAAGGCGCGCTGGAGAACCTTGAATCCACCTCGCCGGAGGGGTCAACCGAGAAGATGGCAGAACGGTAAGCAGAGCTGCGTCCGGATATCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGATTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGGTTGGTGGTATGACGCCGGGTTGTTCGGCCCGCGCCAGCGAACt                                                                                                                      >  2:426057/1‑251 (MQ=255)                                                                                                                 ggCGAGCTGGAGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCATAGAGGCGTCCGGATCTCACTTTCCCGCAGCATGTGGGCGAACAGTCGGACGACAACGTTGCGTACTGCCCGGAACGTGTGGTAGCAGGACAGGTAATGGTCGAGGTGATTAAAAACGATCGCGTGATTGGTGGTATGACGTCGGTTTGGTCGGCCCGCGCCAGCCAACGGTACAAAAtttt                                                                                                          <  2:151139/251‑1 (MQ=255)                                                                                                                   cgcgcTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTcc                                                                                                        <  1:110127/251‑1 (MQ=255)                                                                                                                                   gattCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTGCGGATCTCTCTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGGCAACATGGCTTACTGCCCGGAACGCGTGTTACCAGGACATGTAATGGTCGAGCTGATTAAAAACGAGCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCGCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAgtgtg                                                                                         <  2:170408/247‑1 (MQ=255)                                                                                                                                                              gagaAGATTGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCGGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTGTGTTCGGCGCGCGCCAGCGAACTGGACAAAAGTTGCCTCGAATGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAt                                                             <  1:147297/251‑1 (MQ=255)                                                                                                                                                              gagaAGAAGGCAGAATGGTTAGCCGAGTTGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGATTACTGCCCGGCACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGAGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGCACTGTACAAAATTTTCATCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAAt                                                             >  2:186274/1‑251 (MQ=255)                                                                                                                                                                        aaGAATGGTTAGCAGAGATACGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGCCGTCAACATTGCTTAATGCCCGGCACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGAGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCAAGCGAACTGTACAAAATTTTCCTCGAAGGGGAGTGTGTCGTCACTAACTCGCGGCCGGCGGACCTGTGTAAGCTc                                                   >  2:276598/2‑251 (MQ=255)                                                                                                                                                                            aTGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCg                                               >  1:180780/1‑251 (MQ=255)                                                                                                                                                                                  aGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGGCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGATCACCGAAAAc                                          >  2:328643/1‑250 (MQ=255)                                                                                                                                                                                     agagATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGTAACGCGTGTTACCAGGACCGGAAATGGTCGGGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCATCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCt                                      >  1:14314/1‑251 (MQ=255)                                                                                                                                                                                           gCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCcgcg                                >  1:355697/1‑251 (MQ=255)                                                                                                                                                                                           gCGTCCGGATCTCACTGTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCcgcg                                <  1:198810/251‑1 (MQ=255)                                                                                                                                                                                             gTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGAt                              >  1:41705/1‑251 (MQ=255)                                                                                                                                                                                               ccGGATCTCACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGAtg                             <  1:2978/250‑1 (MQ=255)                                                                                                                                                                                                       cACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTGCTGCCCGGAACGCGTGTGACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCTCGGACGGCGGAAATGTGTAAGCTAACCGAAAACAGATTCCGCGATGTGAATATc                     <  2:37143/250‑1 (MQ=255)                                                                                                                                                                                                       cACTTTCCCGCAGCAGGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGATTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAat                        <  1:97985/247‑1 (MQ=255)                                                                                                                                                                                                                   gcagGTGGGCGAACAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTTGTATGACGCCGGTTTGTGCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATg        <  2:158080/251‑1 (MQ=255)                                                                                                                                                                                                                     agGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGAGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGACTATCGCTTTTGCTATTGa       >  2:127798/1‑250 (MQ=255)                                                                                                                                                                                                                         gggCGAACCGGCGGACGTCAACATTGCTTACGTCCCGGAACGCGTGTTCCCAGGACAGGTAATGGTCTAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTTATGAATTGt  <  2:329167/251‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     TGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGCGCCAGTGCTGAAAAAAGGCGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGTGATTGGTGGTATGACGCCGGTTTGTTCGGCCCGCGCCAGCGAACTGTACAAAATTTTCCTCGAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGT  >  NC_000913/3971481‑3971947

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