Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,824,767	G→A	100%	C65C (TGC→TGT)	kguE ←	epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,824,767	0	G	A	100.0%	147.1 / NA	45	C65C (TGC→TGT)	kguE	epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCACATCGCAGTGCTCATCGAAGTGCCCCAAGGAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGGCATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTGCTCCTGGCCACGTTGGCGCA  >  NC_002947/3824624‑3824909                                                                                                                                                |                                                                                                                                               gTCACATCGCAGTGCTCATCGAAGTGCCCCAAGGAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTc                                                                                                                                           >  5:236965/1‑149 (MQ=255)                   tCGAAGTGCCCCAAGGAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAAt                                                                                                                         <  3:53420/149‑1 (MQ=255)                       aGTGCCCCAAGGAAACCTTCAGCGCCACTGCGCCACGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGAcc                                                                                                                     >  6:199648/1‑149 (MQ=255)                            cccAAGGAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGcc                                                                                                                >  4:264122/1‑149 (MQ=255)                                 ggAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAgcg                                                                                                           <  3:264122/149‑1 (MQ=255)                                  gAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCTTGCCCa                                                                                                              <  6:240880/145‑1 (MQ=255)                                  gAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCTTGCCCa                                                                                                              >  5:240880/1‑145 (MQ=255)                                  gAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCa                                                                                                              <  4:321844/145‑1 (MQ=255)                                  gAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCa                                                                                                              >  3:321844/1‑145 (MQ=255)                                   aaaCCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCa                                                                                                         <  1:164331/149‑1 (MQ=255)                                     aCCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGc                                                                                                       <  6:236965/149‑1 (MQ=255)                                      ccTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCaccg                                                                                                      >  6:267559/1‑146 (MQ=255)                                      ccTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCt                                                                                                      <  3:313164/149‑1 (MQ=255)                                       cTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTg                                                                                                     <  2:293365/149‑1 (MQ=255)                                        ttCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCCCGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGc                                                                                                    >  3:250216/1‑149 (MQ=255)                                        ttCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGc                                                                                                    >  5:316522/1‑149 (MQ=255)                                        ttCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGc                                                                                                    >  2:146865/1‑149 (MQ=255)                                         tCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGcc                                                                                                                >  7:40992/1‑136 (MQ=255)                                         tCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGcc                                                                                                                <  8:40992/136‑1 (MQ=255)                                                 aCTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTgc                                                                                           >  6:188292/1‑149 (MQ=255)                                                 aCTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTgc                                                                                           >  3:179477/1‑149 (MQ=255)                                                      gcCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTgg                                                                                      <  1:167265/149‑1 (MQ=255)                                                        cAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCgg                                                                                                                             >  1:325208/1‑108 (MQ=255)                                                        cAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCgg                                                                                                                             <  2:325208/108‑1 (MQ=255)                                                             gcCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATg                                                                                                                        <  8:147017/108‑1 (MQ=255)                                                             gcCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATg                                                                                                                        >  7:147017/1‑108 (MQ=255)                                                             gcCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACaa                                                                               >  6:12161/1‑149 (MQ=255)                                                                   gCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAAttcttc                                                                         <  5:199648/149‑1 (MQ=255)                                                                       tttCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTcgc                                                                       <  6:149817/147‑1 (MQ=255)                                                                        ttCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCgcgc                                                                     >  4:147147/1‑148 (MQ=255)                                                                          ccAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAAt                                                                  >  7:57122/1‑149 (MQ=255)                                                                           cAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAAt                                                                  >  4:248287/1‑148 (MQ=255)                                                                                 cTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGAc                                                             <  3:147147/147‑1 (MQ=255)                                                                                    cACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCgg                                                         <  2:79694/148‑1 (MQ=255)                                                                                         gcagcGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTg                                                                                                                  >  4:80363/1‑86 (MQ=255)                                                                                         gcagcGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTg                                                                                                                  <  3:80363/86‑1 (MQ=255)                                                                                          cagcGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGc                                                  <  5:188292/149‑1 (MQ=255)                                                                                                     gCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGc                                        >  2:262369/1‑148 (MQ=255)                                                                                                         gCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCcagca                                   <  1:262369/149‑1 (MQ=255)                                                                                                                 cgGTCCACAGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGa                           <  6:316522/149‑1 (MQ=255)                                                                                                                      cacaGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTg                      >  7:196134/1‑149 (MQ=255)                                                                                                                      cacaGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTg                      <  8:141460/149‑1 (MQ=255)                                                                                                                         aGTTCAAGGGGGGTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTGCTc                   <  1:146865/149‑1 (MQ=255)                                                                                                                                ggggggTGGAGTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTGCTCCTGGCCa            <  5:12161/149‑1 (MQ=255)                                                                                                                                          gTACAGACATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTGCTCCTGGCCACGTTGGCGca  >  5:169185/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                               GTCACATCGCAGTGCTCATCGAAGTGCCCCAAGGAAACCTTCAGCGCCACTGCGCCAAGCGCCCGGGCTATTTCCAGCTTCTGCACCAGCAGCGGGTCGGGCACGCCTTGCGCGGTCCACAGTTCAAGGGGGGTGGAGTACAGGCATTCCAGGCGCAGCGCGGCAATGACCCTGCCCAGCGCAGCTGCGTCCGGTGCGCTGGCGAACAATTCTTCGCGCAATTCGACGCGGCTGACGCCAGCGGCTGCCAGCAGGTCCAGAAACTGCTCCTGGCCACGTTGGCGCA  >  NC_002947/3824624‑3824909

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