Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,936,855	T→G	A424A (GCT→GCG)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,936,855	0	T	G	96.5%	89.4 / ‑5.3	29	A424A (GCT→GCG)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (16/12);  minor base C (0/1);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCA  >  NC_000913/2936616‑2937098                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGTAAAGCTCTCAGCAAGAGGCCGACGCTTGCCt                                                                                                                                                                                                                                          <  1:66718/251‑1 (MQ=255)          aTGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTAc                                                                                                                                                                                                                                 >  2:16000/1‑251 (MQ=255)                            ggCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGAt                                                                                                                                                                                                               >  1:121737/1‑251 (MQ=255)                                              gccgagttgCGTTCCTACTGGGCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCTGCCATGTAATCAACTCCGGTTCTGCTGCGCTGGTCGGTTCCTGTAAACATCGCGACAGCGAAGGTAACCGGTCGATGAAGCCACACTGGGAAATCTTTCAGCATGAGGCGGTCGCGTGCCTCGCCGCGTCCTAATTGTGCCCGGCGATCCACGAATACTTCCGTgg                                                                                                                                                                                             <  1:106560/244‑1 (MQ=255)                                                 gATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGACCACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGCGCCCGGCGATCCACGAATACTTCCGTggcgg                                                                                                                                                                                          >  1:145186/1‑251 (MQ=255)                                                    gggCGTACCTACGGGTCACAAGAAGAAATTGAGCGTGTAACGGGGCATAACCTGGATGGACTGGCAGAACACGGCATAATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTACCCCGACGATGAAGCCACACGGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCGCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGcggcgg                                                                                                                                                                                       <  2:80406/249‑1 (MQ=255)                                                        gTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTAc                                                                                                                                                                                   >  2:51856/1‑251 (MQ=255)                                                                       cAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGCACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGCAGGTAACCCGCCGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGCTACTCTTTCCGccccctt                                                                                                                                                                    >  1:235017/1‑245 (MQ=255)                                                                            gCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCg                                                                                                                                                                >  1:193583/1‑250 (MQ=255)                                                                              aaTTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGGTACTTTTCCCGTTTCCtctccgatg                                                                                                                                                             >  2:7285/1‑243 (MQ=255)                                                                                                               aCTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAAc                                                                                                                             >  1:239179/1‑250 (MQ=255)                                                                                                                    aaGAACCCGGCTTCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTTCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACatcatc                                                                                                                       <  1:108775/250‑1 (MQ=255)                                                                                                                         cacGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAgg                                                                                                                  <  1:210311/251‑1 (MQ=255)                                                                                                                           gggCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTCCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAgaac                                                                                                                <  2:91741/250‑4 (MQ=255)                                                                                                                                tcatcCATATGATCATCTCAGGTTCTGCTCCGCTGGACGGATCCTGTCAACCACGCGACAGAGAAGGTAACCAGACGATGAAGACACAATGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGACTGGTGCCCAGCGATCCACGAATACTTCCGTGGCGGCGGTTACTGCTCCACCTTCCTTACCGAAGGCGGAGTCCCGTGCCCCCTGCCTCGCGTCCATATCATCAAAGGCCTgggg                                                                                                           >  2:137458/1‑250 (MQ=255)                                                                                                                                  acccATTTGACCAACCCCGGTTCTGCTGCGCTGGACGGTTCCTGAAAACAACGCGCCAGCGAAGGTACCCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGAcc                                                                                                         <  2:193583/249‑1 (MQ=255)                                                                                                                                     cATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCCAAGGCCTGGGACCGGt                                                                                                      >  2:58936/1‑251 (MQ=255)                                                                                                                                        atgtgcactcccGGTTCTGATGCGCGGGAGGTTTCCTGTAAACAGGGCGACAGGGAGGGTACCCCGAGGATGAAGCCACACTGGGAAATCTCTTCGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTAc                                                                                                    <  2:229605/240‑1 (MQ=255)                                                                                                                                                 ccgggTTCTGGTGCGCTGGACGTTTCCTGAAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAACCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTGACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATcgc                                                                                          <  2:151745/248‑1 (MQ=255)                                                                                                                                                              cgcggaaCGGTCCCTGAAAACAACGCGCGAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGc                                                                             <  2:66110/245‑1 (MQ=255)                                                                                                                                                                  ggccggttctgtaaaaaaagggaaCAGAGAAGGTACCCCGACGATGAAGCCACACGGGGAAATCTCTCAGCAGGAGGCGGACGCTGGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTgg                                                                         <  2:170040/228‑1 (MQ=255)                                                                                                                                                                                aacaacGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGAt                                                           >  2:59665/1‑251 (MQ=255)                                                                                                                                                                                      gcgACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCAt                                                     >  1:222696/1‑251 (MQ=255)                                                                                                                                                                                        gACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATg                                                    >  2:196073/1‑250 (MQ=255)                                                                                                                                                                                              gAAGGTAACCCGACGATGAAGCCACACTGGGAAATCGCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGGGGAATTGCCGAAGGATGTGCATGACATCCt                                             <  1:101397/251‑1 (MQ=255)                                                                                                                                                                                                             aTGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCCACATCATCAAAGGCCTGTGACCGGTACTGCCAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGGGCCTGACATCCCCAACAAACGCACCa                               >  2:93982/1‑250 (MQ=255)                                                                                                                                                                                                             aTGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGGTGTGCATGACATCCTCAACAAACGCACCa                               >  2:95834/1‑250 (MQ=255)                                                                                                                                                                                                                        cTGGGAAATCTCTCAGCAAGAGGCGGACGCTTGCCTAGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACAGCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAACCGCACCAACTCAACCTGGc                   >  2:61260/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaGAGGCGGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCa  <  1:188990/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GCAATGCTAATGGGTCACCAGCTCACCGGCACCGCTCAGGTATTTGCCGATGTGCGTACCTACTGGTCACCAGAAGCAATTGAGCGTGTAACGGGGCATAAACTGGATGGACTGGCAGAACACGGCATCATCCATTTGATCAACTCCGGTTCTGCTGCGCTGGACGGTTCCTGTAAACAACGCGACAGCGAAGGTAACCCGACGATGAAGCCACACTGGGAAATCTCTCAGCAAGAGGCTGACGCTTGCCTCGCCGCTACCGAATGGTGCCCGGCGATCCACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCA  >  NC_000913/2936616‑2937098

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