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%	106.5 / NA	34	G112S (GGC→AGC)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/17);  total (17/17)

GCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTGA  >  NC_000913/2935795‑2936050                                                                                                                           |                                                                                                                                      gCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                       <  2:530414/139‑1 (MQ=255) gCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                       <  1:826001/139‑1 (MQ=255) gCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                       <  1:432269/139‑1 (MQ=255)                 cagTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCa                                                                                                       >  1:352069/1‑139 (MQ=255)                              ggCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAg                                                                                          >  1:521197/1‑139 (MQ=255)                                    aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcaca                                                                                    >  2:583564/1‑139 (MQ=255)                                                       gCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGcc                                                                                 <  2:308468/123‑1 (MQ=255)                                                       gCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGcc                                                                                 >  1:308468/1‑123 (MQ=255)                                                            tGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCt                                                                                                                   >  2:474157/1‑84 (MQ=255)                                                            tGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCt                                                                                                                   <  1:474157/84‑1 (MQ=255)                                                             gCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATtct                                                            >  1:786816/1‑138 (MQ=255)                                                             gCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATtct                                                            <  2:786816/138‑1 (MQ=255)                                                                   gCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTa                                                     <  1:583564/139‑1 (MQ=255)                                                                      gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACgg                                                  >  1:109458/1‑139 (MQ=255)                                                                      gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACgg                                                  >  2:160894/1‑139 (MQ=255)                                                                        caaaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTc                                                >  2:524287/2‑139 (MQ=255)                                                                         aaaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                          >  1:38633/1‑112 (MQ=255)                                                                         aaaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                          <  2:38633/112‑1 (MQ=255)                                                                                    aTGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                    <  1:160894/139‑1 (MQ=255)                                                                                    aTGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                    <  2:521197/139‑1 (MQ=255)                                                                                              cccgcccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGa                          <  1:460774/139‑1 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGccc                                                                                                                          >  1:261142/1‑39 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGccc                                                                                                                          <  2:261142/39‑1 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                       >  2:277717/1‑42 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCgg                                                                                                                       <  1:277717/42‑1 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                                     >  2:171636/1‑44 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                                     <  1:171636/44‑1 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACGTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGAcaca                      >  1:566146/1‑139 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCAt                                                               <  2:56546/98‑1 (MQ=255)                                                                                                  cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCAt                                                               >  1:56546/1‑98 (MQ=255)                                                                                                       aGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGAt                                  >  2:102568/1‑122 (MQ=255)                                                                                                       aGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGAt                                  <  1:102568/122‑1 (MQ=255)                                                                                                                  ggCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATg      >  1:387550/1‑139 (MQ=255)                                                                                                                      ttAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTga  <  2:566146/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                      GCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTGA  >  NC_000913/2935795‑2936050

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