Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,937,210	C→T	L543F (CTC→TTC)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,210	0	C	T	90.0%	72.2 / 1.0	30	L543F (CTC→TTC)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (11/16);  total (12/18)
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.08e-01

CGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATAT  >  NC_000913/2936967‑2937454                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTc                                                                                                                                                                                                                                               >  1:189023/1‑251 (MQ=255)    gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCa                                                                                                                                                                                                                                             <  1:105256/250‑1 (MQ=255)         catcatCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCAtgctgc                                                                                                                                                                                                                                       >  1:233674/1‑251 (MQ=255)                 aaGGGCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCTGTTACTGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCg                                                                                                                                                                                                                               <  2:6139/251‑1 (MQ=255)                  aGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGCATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCgg                                                                                                                                                                                                                              >  1:250380/1‑251 (MQ=255)                        gggACCGGTACTGAAAATCTCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGgtatgta                                                                                                                                                                                                                        <  2:139649/251‑1 (MQ=255)                                         gcgcgGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGAGGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAagag                                                                                                                                                                                                       <  2:43652/250‑1 (MQ=255)                                         ccgcgGAAGGCTGGATCGTCGAATTGCCGAAGGATGTGAATGACATCACCATCAAACGCTCCAACTCAACCTGGCCAACAATCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATTGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATTCACAACGTTGAagag                                                                                                                                                                                                       <  2:233674/250‑1 (MQ=255)                                          cgcgGAAGGCTGGAGCGAGGCCTTGCCGAAGGATGTGCATGACATCCACAACATACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGCCCATCGGCCCCCTTGGCGACGACTTTATCACTTTCGCCTCCATGCTGCGTCTCCCAGTATGTATGCACCACGTTgcaacg                                                                                                                                                                                                       >  2:271514/1‑245 (MQ=255)                                                        gCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGt                                                                                                                                                                                        >  1:200683/1‑251 (MQ=255)                                                               aTTGCCGAAAGATGTTCATGACAGCCTCAACAAACGCACCAACACAACCGGGCCAACCACCGGGTTTGCACCGCGCCTCGCCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGTGCGCTAACCATGGGGGTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCGGCGTATCCCGGTATGTATGCACAACGTGTAAGAGACCAAAGTGTATCGTCCTTCTg                                                                                                                                                                                 <  2:148025/251‑1 (MQ=255)                                                                                     aaccacaacaaACCCCCCAAACCAACCTGGCCCACCCCCGGGTTTGCACCGCTCCGCACCGGTAAGGGGTCGTTTACGGAAGTGTACGTTGTAAAGGCATACTGGGGCTCTAACCATGGGGTTCTGACCATCGGCCACGTCTGCGCAGACTTTATCAGTTTCGCCTCCATGCTGCGTGTGCCGGTATGTGTGCACAACGTTGAATAGACCAAATTGTATCGGCTTTCTGCCTGGGCAGCGCGCTGCATGGa                                                                                                                                                           <  2:272389/246‑1 (MQ=255)                                                                                               aaCGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGc                                                                                                                                                 >  1:162466/1‑251 (MQ=255)                                                                                                       aaCTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTa                                                                                                                                         <  1:261265/251‑1 (MQ=255)                                                                                                                   ccAACCGCCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCa                                                                                                                             <  1:240402/251‑1 (MQ=255)                                                                                                                     aaCCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGcc                                                                                                                              >  1:51512/1‑248 (MQ=255)                                                                                                                     aaCCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCGCAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGcc                                                                                                                              >  2:32876/1‑248 (MQ=255)                                                                                                                        caccTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGCGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCGGGATTACCGCGCTTGCCAGAACt                                                                                                                        >  2:51310/1‑251 (MQ=255)                                                                                                                                             cACCGGTAAAGGGCCTTTGACGGATGTGTCCTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACACCGTTGAAGAGACCACAGTGTCTCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATCTTGCAGGCCAGGATTACCGTGCTTGCCAGAACTCCGGTCCGtgtgaacaactat                                                                                                   >  2:179064/1‑239 (MQ=255)                                                                                                                                               ccGGTCAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTa                                                                                                  >  1:245129/1‑250 (MQ=255)                                                                                                                                                       gggCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAgg                                                                                                                                             <  1:78661/199‑1 (MQ=255)                                                                                                                                                         gctgTTTACGGATGTGTACCCGGTAATGGCGAACTGGGGCGCTAACCACGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTGTAGCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTAGCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCAAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTc                                                                                       <  2:268244/248‑1 (MQ=255)                                                                                                                                                                   gATGTGTCCTCGTTAATTGCGAACTGGGGCGCGTACCATGGGGTTCTGACCATCGCCCTCTTTGGTGCAGACTTTATCACTTTCGCCTCCAGGCTGCGTATCCCGTTATGTATGCACAACGTTGAAGAGACCAATGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTTCGTTGTCCAAGCTTTAATAAAGTACTCCCtgatgtga                                                                             <  2:287627/251‑1 (MQ=255)                                                                                                                                                                             cGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCg                                                                   <  1:51310/251‑1 (MQ=255)                                                                                                                                                                                                       cATGGGGTTCTGTCCATCGGCCACGTTGGCGCAGTCTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATCCACAACGTTGAAGAGACCAAAGTCTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATAGTGAAGGCCAGGATTTCCGCTCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAAcaagca                                         <  1:159627/251‑1 (MQ=255)                                                                                                                                                                                                       cATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAAcaagca                                         <  1:123439/251‑1 (MQ=255)                                                                                                                                                                                                          ggggTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAAcaagcaag                                       <  1:102182/250‑1 (MQ=255)                                                                                                                                                                                                                  gACCATCGGCCACTGTGGCTCAGACTTTATCACTTGCGCCTCCAGGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGGCCTTCTGCCTGGGCTGCGCACGGCAGGGATATTGAAGGCCATGTTTACCTCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGGTAATACATTACACCCTGATGTGATGCCCGGCCGCTCCGTCTTCCGGGCCTGAACACTCAAGAGTGGTTa                               <  2:200683/250‑1 (MQ=255)                                                                                                                                                                                                                      aTCGGCCACGTTGGCGCAGACTTTATCACTTTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGAACGGCATGGATATTGAAGGCCAGGATTACCGAGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGAAGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCgg                          >  1:257030/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCGCAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTatat  <  1:32876/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGGCGCAGACTTTATCACTCTCGCCTCCATGCTGCGTATCCCGGTATGTATGCACAACGTTGAAGAGACCAAAGTGTATCGTCCTTCTGCCTGGGCTGCGCACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATAT  >  NC_000913/2936967‑2937454

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