Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,935,917	G→A	G112S (GGC→AGC)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,935,917	0	G	A	100.0%	126.8 / NA	40	G112S (GGC→AGC)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/20);  total (20/20)

GCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGT  >  NC_000913/2935785‑2936047                                                                                                                                     |                                                                                                                                   gCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCa                                                                                                                                 <  1:686591/136‑1 (MQ=255) gCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCa                                                                                                                                 >  2:686591/1‑136 (MQ=255)     ctgctTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACg                                                                                                                          >  2:1073332/1‑139 (MQ=255)                                          ccTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCt                                                                                     >  2:438210/1‑139 (MQ=255)                                          ccTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCt                                                                                     <  1:1345606/139‑1 (MQ=255)                                            tCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                                                                                   >  1:60656/1‑139 (MQ=255)                                             cACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                  <  1:1073332/139‑1 (MQ=255)                                                        tAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                       >  1:1047399/1‑139 (MQ=255)                                                          aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                     >  1:269774/1‑139 (MQ=255)                                                                tGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCg                                                                                                                     >  1:1124341/1‑85 (MQ=255)                                                                tGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCg                                                                                                                     <  2:1124341/85‑1 (MQ=255)                                                                tGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAg                                                               <  1:620505/139‑1 (MQ=255)                                                                          aTGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                            >  2:804985/1‑68 (MQ=255)                                                                          aTGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                            <  1:804985/68‑1 (MQ=255)                                                                                gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                            >  2:723692/1‑62 (MQ=255)                                                                                gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                            <  1:723692/62‑1 (MQ=255)                                                                                gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                  >  2:632508/1‑104 (MQ=255)                                                                                gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                  <  1:632508/104‑1 (MQ=255)                                                                                  gAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTc                                             >  1:759148/1‑139 (MQ=255)                                                                                   aaaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCa                                            <  2:1047399/139‑1 (MQ=255)                                                                                        aTCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAATAGCACTGaa                                                                                                                            >  1:970451/1‑54 (MQ=255)                                                                                        aTCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                            <  2:970451/54‑1 (MQ=255)                                                                                         tCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                    >  2:155936/1‑61 (MQ=255)                                                                                         tCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                    <  1:155936/61‑1 (MQ=255)                                                                                         tCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                     >  2:483366/1‑108 (MQ=255)                                                                                         tCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                     <  1:483366/108‑1 (MQ=255)                                                                                              aTGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                 <  2:759148/139‑1 (MQ=255)                                                                                                       aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCAt                                                     <  1:845091/110‑1 (MQ=255)                                                                                                       aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCAt                                                     >  2:845091/1‑110 (MQ=255)                                                                                                        cccgcccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGTCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGa                       <  2:60656/139‑1 (MQ=255)                                                                                                            cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACg                                                >  1:227790/1‑110 (MQ=255)                                                                                                            cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACg                                                <  2:227790/110‑1 (MQ=255)                                                                                                            cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                 <  2:1073116/125‑1 (MQ=255)                                                                                                            cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                 >  1:1073116/1‑125 (MQ=255)                                                                                                            cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGAcaca                   >  2:748774/1‑139 (MQ=255)                                                                                                                             gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAg                                                                                                   >  1:12938/1‑42 (MQ=255)                                                                                                                             gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAg                                                                                                   <  2:12938/42‑1 (MQ=255)                                                                                                                             gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGt  <  2:269774/139‑1 (MQ=255)                                                                                                                                 tAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGtt                                     >  1:1165699/1‑100 (MQ=255)                                                                                                                                 tAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGtt                                     <  2:1165699/100‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   GCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGT  >  NC_000913/2935785‑2936047

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