Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,931,836	A→G	D337G (GAT→GGT)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,931,836	0	A	G	100.0%	153.7 / NA	42	D337G (GAT→GGT)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/23);  total (19/23)

GGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCG  >  NZ_CP009273/2931750‑2931911                                                                                       |                                                                            ggAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGc                                                                         <  2:1051757/91‑1 (MQ=255) ggAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGc                                                                         <  1:992735/91‑1 (MQ=255) ggAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGc                                                                         <  1:811505/91‑1 (MQ=255) ggAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGc                                                                         <  1:1121490/91‑1 (MQ=255)      gggTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTaaa                                                                    >  2:127475/1‑91 (MQ=255)      gggTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTaaa                                                                    >  1:44913/1‑91 (MQ=255)           agaAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGc                                                               <  2:623557/91‑1 (MQ=255)               aaCTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAgtg                                                           <  2:767671/91‑1 (MQ=255)               aaCTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAgtg                                                           <  2:528679/91‑1 (MQ=255)                aCTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAgtgt                                                          <  1:384585/91‑1 (MQ=255)                aCTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAgtgt                                                          >  1:1087065/1‑91 (MQ=255)                 cTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGtgtg                                                         >  1:489710/1‑91 (MQ=255)                  tGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCTGGTGGCGTAAAAATGCAGTGTGa                                                        <  2:334187/91‑1 (MQ=255)                  tGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGa                                                        <  2:602801/91‑1 (MQ=255)                      cTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTAATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTa                                                    >  1:859068/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCGGTGTGATTTATTGTCgt                                             >  1:932346/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCGGTGTGATTTATTGTCgt                                             >  1:293040/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  1:708209/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  1:1101539/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             <  1:127475/91‑1 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  1:988098/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  2:779888/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  2:582143/1‑91 (MQ=255)                             gCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCgt                                             >  2:156384/1‑91 (MQ=255)                                gAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTc                                          <  2:397264/91‑1 (MQ=255)                                       ccTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACgc                                   >  1:51698/1‑91 (MQ=255)                                        cTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACgct                                  <  1:768663/91‑1 (MQ=255)                                                tGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCa                          >  2:201345/1‑91 (MQ=255)                                                    gattgaAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGa                      >  1:37687/1‑91 (MQ=255)                                                    gattgaAGAAGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGa                      >  2:794808/1‑91 (MQ=255)                                                        gaagaaGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCGGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  2:293040/91‑1 (MQ=255)                                                        gaagaaGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  2:552367/91‑1 (MQ=255)                                                        gaagaaGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  1:877183/91‑1 (MQ=255)                                                        gaagaaGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  1:956680/91‑1 (MQ=255)                                                        gaagaaGCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  1:1083239/91‑1 (MQ=255)                                                           gaagCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGc               <  1:1025527/91‑1 (MQ=255)                                                             agCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCtt             <  2:917176/91‑1 (MQ=255)                                                             agCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCtt             <  2:221921/91‑1 (MQ=255)                                                              gCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  >  2:1030247/1‑85 (MQ=255)                                                              gCTCGTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGa                  <  1:1030247/85‑1 (MQ=255)                                                                  gTCTGATCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATAc        >  1:396279/1‑91 (MQ=255)                                                                        tCGCGCCTGGCGCGGGTGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCAcgcg  <  2:446681/91‑1 (MQ=255)                                                                                       |                                                                            GGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCG  >  NZ_CP009273/2931750‑2931911

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