Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	3,236,935	C→T	G389D (GGT→GAT)	uxaC ←	uronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,236,935	0	C	T	100.0%	147.7 / NA	44	G389D (GGT→GAT)	uxaC	uronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/20);  total (24/20)

TCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATACCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCTTCGCC  >  CP009273/3236870‑3236999                                                                  |                                                                 tcagCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCg                                                                >  1:477091/1‑68 (MQ=255)  cagCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGt                                                               <  1:454120/68‑1 (MQ=255)  cagCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGt                                                               <  1:3929937/68‑1 (MQ=255)  cagCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGt                                                               <  1:2601210/68‑1 (MQ=255)  cagCATACCGACAAAGCGGCGCAGCAGACCGCGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGt                                                               <  1:420500/68‑1 (MQ=255)   agCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATTTCGTc                                                              >  1:3249267/1‑68 (MQ=255)   agCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTc                                                              >  1:4029129/1‑68 (MQ=255)   agCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTc                                                              >  1:2244898/1‑68 (MQ=255)   agCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTc                                                              <  1:5991187/68‑1 (MQ=255)     cATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCtt                                                            <  1:864987/68‑1 (MQ=255)           gACAACGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGCCGTTCCATATCGTCTTTCTGaa                                                      >  1:3205414/1‑67 (MQ=255)           gACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTTGGTCATCTGACGTTCCATATCGTCTTTCTGAt                                                      >  1:5409816/1‑68 (MQ=255)           gACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGAt                                                      >  1:1279765/1‑68 (MQ=255)           gACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGAt                                                      >  1:4858667/1‑68 (MQ=255)            aCAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGGCATCTGACGTTCCATATCGTCTTTCTGATc                                                     <  1:1661812/68‑1 (MQ=255)             cAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCg                                                    <  1:579070/68‑1 (MQ=255)             cAAAGCGGCTCAGCAGACCGAGCAGCGCCAGATGGGTCATCTCACGTTCCATATCGTCTTTCTGATCg                                                    <  1:284481/68‑1 (MQ=255)              aaaGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGt                                                   >  1:3513326/1‑68 (MQ=255)               aaGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGtt                                                  <  1:3413290/68‑1 (MQ=255)               aaGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGtt                                                  <  1:2843865/68‑1 (MQ=255)               aaGCGGCGCAGCAGACCGAGCTGCGCCAGTGGGGGCATCTGACGTTCCATATCGTCTTTCTGATCGtt                                                  <  1:2659431/68‑1 (MQ=255)                aGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTa                                                 >  1:3019706/1‑68 (MQ=255)                aGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTa                                                 >  1:4357539/1‑68 (MQ=255)                      tCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAaccac                                           >  1:1261575/1‑68 (MQ=255)                              ccGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGa                                   >  1:3772444/1‑68 (MQ=255)                                gAGCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAAc                                 >  1:3403421/1‑68 (MQ=255)                                  gCTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCg                               >  1:3808967/1‑68 (MQ=255)                                   cTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGa                              <  1:6001820/68‑1 (MQ=255)                                   cTGCGCCAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGa                              <  1:893506/68‑1 (MQ=255)                                         cAGTTGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCa                        <  1:466802/68‑1 (MQ=255)                                            ttGGGTCTTCTGCCGTTCCATATCGTCGTTCTGATCGTTAACCCACCAGCCGGAACCGAACTGCAttt                     >  1:4154448/1‑68 (MQ=255)                                            ttGGGTCATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCAttt                     >  1:5452531/1‑68 (MQ=255)                                               ggTCATCTGACGATCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGc                  <  1:1931384/68‑1 (MQ=255)                                                  cATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCg               >  1:4654791/1‑68 (MQ=255)                                                  cATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCg               >  1:3871032/1‑68 (MQ=255)                                                  cATCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCg               >  1:727134/1‑68 (MQ=255)                                                    tcagACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGc             <  1:1682267/65‑1 (MQ=255)                                                    tCTTACGTTCCATATCGTCTTTCTGATCGTTAAACCACCTGCGGGAACCGAACTGCATTTTGCGCGGc             <  1:499003/68‑1 (MQ=255)                                                    tCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGc             <  1:5522954/68‑1 (MQ=255)                                                    tCTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGc             <  1:2764657/68‑1 (MQ=255)                                                     cTGACGTTCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCa            >  1:3608760/1‑68 (MQ=255)                                                           ttCCATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCtt      >  1:5895613/1‑68 (MQ=255)                                                           ttCCATATCGTCGTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCtt      >  1:1751023/1‑68 (MQ=255)                                                              cATATCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCTTCGc   >  1:2850201/1‑68 (MQ=255)                                                               atatCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGACCTGCATTTTGCCCGGCATACCTTCGcc  <  1:4824408/68‑1 (MQ=255)                                                               atatCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCTTCGcc  <  1:790700/68‑1 (MQ=255)                                                               atatCGTCTTTCTGATCGTTAAACCACCAGCCGCAACCGAACTGCATTTTGCCCGGCATACCTTCGcc  <  1:3287233/68‑1 (MQ=255)                                                                  |                                                                 TCAGCATACCGACAAAGCGGCTCAGCAGACCGAGCTGCGCCAGTTGGGTCATCTGACGTTCCATACCGTCTTTCTGATCGTTAAACCACCAGCCGGAACCGAACTGCATTTTGCCCGGCATACCTTCGCC  >  CP009273/3236870‑3236999

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