Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	351,375	A→G	Y54C (TAC→TGC)	prpD →	2‑methylcitrate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	351,375	0	A	G	100.0%	104.1 / NA	30	Y54C (TAC→TGC)	prpD	2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/12);  total (18/12)

GAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTG  >  NC_000913/351245‑351507                                                                                                                                   |                                                                                                                                     gAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTg                                                                                                                             >  1:554010/1‑140 (MQ=255)            tGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTGGACTGCCCGGCCTGTAAAAAActgc                                                                                                                  >  1:280005/1‑140 (MQ=255)            tGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAActgc                                                                                                                  >  2:759030/1‑140 (MQ=255)            tGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAActgc                                                                                                                  >  1:701516/1‑140 (MQ=255)                  cGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGc                                                                                                            <  2:731640/140‑1 (MQ=255)                             tcgATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTccc                                                                                                 >  2:569352/1‑140 (MQ=255)                             tcgATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTccc                                                                                                 >  2:53065/1‑140 (MQ=255)                                 ttACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCa                                                                                             <  1:484100/140‑1 (MQ=255)                                      tCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACcgtc                                                                                        <  1:759030/140‑1 (MQ=255)                                          gAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTAc                                                                                    >  1:815679/1‑140 (MQ=255)                                                     tCAGCTCTAAAGTGGCCTACGACACCGCACAATACAGCCTGCTCGACATGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTg                                                                         <  2:280005/140‑1 (MQ=255)                                                               aGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCg                                                               <  1:635295/140‑1 (MQ=255)                                                               aGTGGCCTACGACACAGCACATTACGGCCTGCTCGACTCGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCg                                                               <  1:368306/140‑1 (MQ=255)                                                                    ccTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACt                                                          <  1:569352/140‑1 (MQ=255)                                                                               gCACATTACTGCCTGCCCGACATGCTCGGCTGCGGGCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAgtt                                               <  2:348963/140‑3 (MQ=255)                                                                                    ttACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGAcc                                          <  1:54293/140‑1 (MQ=255)                                                                                     tACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGACGTACCCACCGCTGTGAGCGTCCCCGGAACTCAGTTCCAGCTCCAccc                                         >  2:303453/1‑140 (MQ=255)                                                                                       ctgcctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCg                                       >  1:812018/1‑140 (MQ=255)                                                                                          cctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTcc                                    >  2:38770/1‑140 (MQ=255)                                                                                          cctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGCTCCAGCTCGACCCCGTcc                                    >  2:450380/1‑140 (MQ=255)                                                                                          cctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTACGCGTCCCCGGAACTCAGTTCAAGTTCGACCCCGTcc                                    >  1:348963/1‑140 (MQ=255)                                                                                             gcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGTGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAgg                                 <  2:618245/140‑1 (MQ=255)                                                                                             gcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAgg                                 <  1:38770/140‑1 (MQ=255)                                                                                               tCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAggcg                               >  1:644274/1‑140 (MQ=255)                                                                                                        tCGGCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCAAGCGCGACCCCGTCCAGGCGGCATTTAAc                      >  1:269179/1‑140 (MQ=255)                                                                                                           gCTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATc                   >  1:723557/1‑140 (MQ=255)                                                                                                            cTGCGGTCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTTCCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCg                  <  2:701516/140‑1 (MQ=255)                                                                                                                  tCTGGAAGCTCTCGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCga            >  1:531249/1‑140 (MQ=255)                                                                                                                            ctcGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCGGCTCGACCGCGTACAGGCGGCATTTAACATTGGCGCGATGATCCGCTg  >  2:315754/1‑140 (MQ=255)                                                                                                                            ctcGAATGCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTg  >  2:673506/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                     GAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTG  >  NC_000913/351245‑351507

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