Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	130,120	T→C	D673D (GAT→GAC)	acnB →	bifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	130,120	0	T	C	100.0%	110.3 / NA	31	D673D (GAT→GAC)	acnB	bifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/20);  total (11/20)

TCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACA  >  NZ_CP009273/130036‑130194                                                                                     |                                                                           tCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCg                                                                      <  2:103815/91‑1 (MQ=255) tCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCg                                                                      <  1:749112/91‑1 (MQ=255)    gggCATGGAAAAATGGCTGGGGAATCCTGAGCTGCTGGAAGCCCATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGat                                                                   >  2:181728/1‑91 (MQ=255)         tGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGCCCTGGCGGATATTaa                                                              >  2:635092/1‑91 (MQ=255)         tGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaa                                                              >  1:786521/1‑91 (MQ=255)          ggAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaaa                                                             <  1:280635/91‑1 (MQ=255)          ggAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaaa                                                             <  1:311511/91‑1 (MQ=255)          ggAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaaa                                                             <  2:619116/91‑1 (MQ=255)          gcaaaaaTGGCTGGCGAATCCTGAGCTGGTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTaaa                                                             <  2:190097/89‑1 (MQ=255)            aaaaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAaga                                                           <  2:613821/91‑1 (MQ=255)              aaaTGGCTGGCGAATCCTGAGCTGCTGGTAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGc                                                         <  1:843273/91‑1 (MQ=255)              aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGc                                                         <  2:65938/91‑1 (MQ=255)              aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGc                                                         >  2:128887/1‑91 (MQ=255)              aaaTGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGc                                                         >  1:782011/1‑91 (MQ=255)                   gctggcGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATc                                                    <  2:138692/91‑1 (MQ=255)                                gAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGa                                       >  2:290757/1‑91 (MQ=255)                                gAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGa                                       <  1:273781/91‑1 (MQ=255)                                  gctgctGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAAc                                     <  1:634050/91‑1 (MQ=255)                                     gctgGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAc                                  >  2:270518/1‑91 (MQ=255)                                       tgGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGAccc                                <  2:852322/91‑1 (MQ=255)                                               gATGCAGATGCGGAGTACGCGGCAGTGATCGACATCGACCTGGCGGATATCAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGAcg                        >  2:94546/1‑91 (MQ=255)                                                        gCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGc               <  2:338451/91‑1 (MQ=255)                                                        gCGGAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGc               <  1:270518/91‑1 (MQ=255)                                                          ggAATACGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTg             <  1:761969/91‑1 (MQ=255)                                                               aCGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGc        <  1:489169/91‑1 (MQ=255)                                                               aCGCGGCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGc        >  2:742971/1‑91 (MQ=255)                                                                    gCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTAc   <  1:581137/91‑1 (MQ=255)                                                                    gCAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTAc   <  2:350670/91‑1 (MQ=255)                                                                     cAGTGATCGACATCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACa  >  1:53040/1‑91 (MQ=255)                                                                                 tCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCg       >  2:389831/1‑74 (MQ=255)                                                                                 tCGACCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCg       <  1:389831/74‑1 (MQ=255)                                                                                     |                                                                           TCAGGGCATGGAAAAATGGCTGGCGAATCCTGAGCTGCTGGAAGCCGATGCAGATGCGGAATACGCGGCAGTGATCGACATCGATCTGGCGGATATTAAAGAGCCAATCCTGTGTGCTCCGAACGACCCGGATGACGCGCGTCCGCTGTCTGCGGTACA  >  NZ_CP009273/130036‑130194

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