Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	645,191	A→G	D252D (GAT→GAC)	citE ←	citrate (pro‑3S)‑lyase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	645,191	0	A	G	100.0%	71.2 / NA	22	D252D (GAT→GAC)	citE	citrate (pro‑3S)‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/12);  total (10/12)

CGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTG  >  NZ_CP009273/645055‑645301                                                                                                                                         |                                                                                                               cGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTc                                                                                                              <  1:385164/139‑1 (MQ=255)        ttCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCActtct                                                                                                       <  1:273945/139‑1 (MQ=255)                gggAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGt                                                                                               <  1:375127/139‑1 (MQ=255)                  gagaGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGGCCACCTCTTTTTGGGGCg                                                                                             >  2:629187/1‑139 (MQ=255)                  gagaGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCg                                                                                             >  1:147530/1‑139 (MQ=255)                   agagCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCgg                                                                                            <  1:695454/139‑1 (MQ=255)                   agagCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCgg                                                                                            <  2:705327/139‑1 (MQ=255)                   agagCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCgg                                                                                            <  1:231010/139‑1 (MQ=255)                     agCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGGTCACCTCTTTCTTGGTCgggt                                                                                          >  2:695454/1‑137 (MQ=255)                       cACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCg                                                                                        <  1:644593/139‑1 (MQ=255)                                              ccGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAAt                                                                 <  1:259559/139‑1 (MQ=255)                                               cGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGGCCACCTCTTTTTGGGGCGGGGCGTAGAGGGTGTGGAGCAGGTCAATc                                                                >  2:259559/1‑139 (MQ=255)                                                ggACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCt                                                               <  2:147530/139‑1 (MQ=255)                                                ggACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGGGCCACTTCTTTCTGGGTCGGGGGGTAGAGGGTGTGCAGGAGGTCAATCt                                                               >  2:137609/1‑139 (MQ=255)                                                        caccaTCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCggg                                                       <  1:137609/139‑1 (MQ=255)                                                                tGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGGGCGTAGAGGGTGTGCAGCAGATCAATCTGACGCGGGTTGATCAg                                               >  1:590716/1‑139 (MQ=255)                                                                tGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGGCCACTTCTTTCTGGGTCGGGGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAg                                               >  1:427118/1‑139 (MQ=255)                                                                               ccACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCa                                  >  1:735688/1‑137 (MQ=255)                                                                                        ggCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGGCCCCTTCTTTCTGGGGCGGGGCGTAGAGGGTGTGCAGCAGGTCAATCTGACGCGGGGTGGTCAGCGCTTTGGCGTCCAAGCCCAGCTg                       >  2:518663/1‑139 (MQ=255)                                                                                              cgcgAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGGCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGAt                 >  1:282169/1‑139 (MQ=255)                                                                                                          cTTCAGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGcttc     <  1:697352/139‑1 (MQ=255)                                                                                                              aGCGGCTTCTACGACGCGGCGGGCGTGGTCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTg  <  2:735688/138‑1 (MQ=255)                                                                                                                                         |                                                                                                               CGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTG  >  NZ_CP009273/645055‑645301

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