Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	351,307	A→G	K31K (AAA→AAG)	prpD →	2‑methylcitrate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	351,307	0	A	G	100.0%	188.7 / NA	54	K31K (AAA→AAG)	prpD	2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/23);  total (31/23)

CAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTG  >  NC_000913/351229‑351384                                                                               |                                                                              caacaaCATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGt                                                                             >  2:79758/1‑81 (MQ=255) caacaaCATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATCACGTCATGAACTACGAAATCAGCTCTAAGGt                                                                             >  2:79757/1‑81 (MQ=255)  aacaacATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTg                                                                            <  2:79698/81‑1 (MQ=255)  aacaacATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTg                                                                            <  2:79697/81‑1 (MQ=255)   acaacaTCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTgg                                                                           >  2:79735/1‑81 (MQ=255)   acaacaTCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTgg                                                                           >  1:79751/1‑81 (MQ=255)         tCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACg                                                                     <  2:79700/81‑1 (MQ=255)         tCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACg                                                                     <  2:79699/81‑1 (MQ=255)          ccGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGa                                                                    >  1:79760/1‑81 (MQ=255)               cGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCg                                                               <  1:79702/81‑1 (MQ=255)               cGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCg                                                               >  2:79762/1‑81 (MQ=255)               cGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCg                                                               <  1:79701/81‑1 (MQ=255)                ggAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCCCGAACTACGAAATCAGCTCTACGGTGGCCTACGTCACCg                                                               >  1:973396/1‑80 (MQ=255)                  aaTTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACTAAATCAGCTCTAAGGTGGCCTACGACACCGcac                                                            <  1:79703/81‑1 (MQ=255)                  aaTTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGca                                                             >  1:79741/1‑80 (MQ=255)                  aaTTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGca                                                             >  1:79738/1‑80 (MQ=255)                   aTTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGcaca                                                           >  1:79749/1‑81 (MQ=255)                     ttGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACAt                                                          >  2:79733/1‑80 (MQ=255)                      tGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACAtt                                                         <  2:79704/80‑1 (MQ=255)                      tGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTa                                                        >  1:79759/1‑81 (MQ=255)                       gATCGTGAAATCGTTGATATCGTCGATTACATCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTa                                                        <  2:79705/80‑1 (MQ=255)                        aTCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTAct                                                      >  1:79744/1‑81 (MQ=255)                           gTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTActgcc                                                   <  2:79706/81‑1 (MQ=255)                             gAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTActgcctg                                                 >  1:79755/1‑81 (MQ=255)                               aaTCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCt                                               <  1:79707/81‑1 (MQ=255)                               aaTCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCt                                               >  1:79776/1‑81 (MQ=255)                                 tCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCg                                             <  2:79708/81‑1 (MQ=255)                                  cGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGa                                            >  2:79748/1‑81 (MQ=255)                                     tGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGacac                                         >  1:79789/1‑81 (MQ=255)                                      gATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGacac                                         <  2:79709/80‑1 (MQ=255)                                      gATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACg                                        >  2:79754/1‑81 (MQ=255)                                        tatCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCt                                      <  1:79711/81‑1 (MQ=255)                                        tatCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCt                                      <  2:79710/81‑1 (MQ=255)                                            gtcgATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGc                                  >  2:79761/1‑81 (MQ=255)                                              cgATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTg                                >  2:79782/1‑81 (MQ=255)                                              cgATTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTg                                >  2:79772/1‑81 (MQ=255)                                                aTTACGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCg                              <  2:79712/81‑1 (MQ=255)                                                    cGTCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCt                          >  1:79785/1‑81 (MQ=255)                                                      tCATGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTgg                        <  1:79713/81‑1 (MQ=255)                                                        aTGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGTCTGCTCGACACGCTCGGTTGCGGTCTgg                        >  1:79771/1‑79 (MQ=255)                                                        aTGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTgg                        >  2:79770/1‑79 (MQ=255)                                                        aTGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGa                       <  2:79714/80‑1 (MQ=255)                                                         tGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAg                     <  1:79715/81‑1 (MQ=255)                                                         tGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAg                     >  2:79764/1‑81 (MQ=255)                                                         tGAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAg                     >  2:79765/1‑81 (MQ=255)                                                          gAACTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGc                    >  2:79756/1‑81 (MQ=255)                                                             cTACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGctct                 <  2:79716/81‑1 (MQ=255)                                                              tACGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCtctc                >  1:79774/1‑81 (MQ=255)                                                                cGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCtctc                <  1:79717/79‑1 (MQ=255)                                                                cGAAATCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCtctc                <  1:79718/79‑1 (MQ=255)                                                                    aTCAGCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATAc          >  1:79763/1‑81 (MQ=255)                                                                        gCTCTAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCgg      <  2:79719/81‑1 (MQ=255)                                                                           ctAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCt   >  1:79781/1‑81 (MQ=255)                                                                            tAAGGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTg  >  1:79767/1‑81 (MQ=255)                                                                               |                                                                              CAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTG  >  NC_000913/351229‑351384

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