Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,336,674	G→A	G282R (GGG→AGG)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,336,674	0	G	A	100.0%	104.9 / NA	31	G282R (GGG→AGG)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/13);  total (18/13)

AATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCAATTTTTACCAGT  >  NC_000913/1336429‑1336913                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aaTGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGaa                                                                                                                                                                                                                                             >  1:97949/1‑250 (MQ=255)    gggtttcttacgcgtatccactctgtGGTGCTGACTCGCACCCCACCATGATCACCGTCCTTGGCG‑GTTGGGGGTGGGGCGTGGGTGGGATGGAAGCAGAGTCCGCAATGTTAGGCCGGCCGGTTTCCATTCTTCCCCCGGATGTAGTGGGCTTCAAACTTCCCGGAAAATGACGTGAGGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATtc                                                                                                                                                                                                                                         <  2:261799/226‑1 (MQ=255)                 ggATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAAGCAGCAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGCTGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATg                                                                                                                                                                                                                            >  1:273119/1‑251 (MQ=255)                 ggATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATg                                                                                                                                                                                                                            >  1:47697/1‑251 (MQ=255)                   aTACCCTCGTTGTTACGGACTCGCACACCACCATGATCACCGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAGGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGTAGTTAGTTGGCTTCAAACTTACCGGAAAATTACGTGAGGGTATTACCGCCACAGACCTGGTTCTCGCTGTTACCCAAATGCTTCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGt                                                                                                                                                                                                                          <  1:187178/251‑1 (MQ=255)                                         gcccaccacCAGGATCAACGGCCTGGGCGTGCTGGGG‑TGGGGCGTGGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGt                                                                                                                                                                                                    <  1:114095/248‑1 (MQ=255)                                               caccagtaacAACGGCCTTGGGGTGTTGGGG‑TGGGGGGTTGGTGGGATCGAAGCGAAACCCGAAAGGTTAGTCCAGCCGTTTTCCATGCTTATCCCGGATGTAGTGGGTTTAAAACTTACAGGAAAATTACGGTAAGTTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCgg                                                                                                                                                                                              <  2:242130/242‑1 (MQ=255)                                                   aTGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGACGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGACTTTTATGGTGATGGTCTGGATTCACTACCGTTGGc                                                                                                                                                                                                >  2:265273/1‑245 (MQ=255)                                                                                ggggCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCCTGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTAGcccga                                                                                                                                                               >  2:408371/1‑247 (MQ=255)                                                                                                 gAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTgg                                                                                                                                             >  1:107871/1‑251 (MQ=255)                                                                                                        aaGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTAAAACTTTCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATTGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCttcttc                                                                                                                                      <  1:456564/251‑1 (MQ=255)                                                                                                              catttttAGGCCAGCCGGTTTCCTTGCTTACCCCTGATGTAGTGGGCTTCAAACTTCCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATc                                                                                                                                <  2:149281/246‑1 (MQ=255)                                                                                                                              ggTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCg                                                                                                                >  1:14098/1‑251 (MQ=255)                                                                                                                                   ccATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTAc                                                                                                           >  1:76947/1‑251 (MQ=255)                                                                                                                                              ccGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCt                                                                                                                                                                                                                  >  1:418621/1‑137 (MQ=255)                                                                                                                                                 gATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCgg                                                                                             >  1:155247/1‑251 (MQ=255)                                                                                                                                                           gCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGa                                                                                    >  1:321328/1‑250 (MQ=255)                                                                                                                                                                     tACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc                                                                          <  1:4089/250‑1 (MQ=255)                                                                                                                                                                     tACCGAAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc                                                                          <  2:97949/250‑1 (MQ=255)                                                                                                                                                                         ggAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGtt                                                                     >  1:280736/1‑251 (MQ=255)                                                                                                                                                                                  cgtgatggaatatccGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCTCTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATTGTGCCACCTGTGGCTTCTTCCCAAGCGATGCTGTATCCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaa                                                            <  1:17563/238‑1 (MQ=255)                                                                                                                                                                                   gTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaa                                                            >  1:303529/1‑250 (MQ=255)                                                                                                                                                                                              aCCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGACACCTGTGGCTTCTTCCCAATCGCTGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAgc                                                >  1:282316/1‑251 (MQ=255)                                                                                                                                                                                               ccGCCACAGACCTGGTTCTCACTGTTTCCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAAGTTTATGGTGATGGTCTGGATTCTCTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTTCCACCTGTGGCTTCTTCCCAATCGGTTGTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGGCGAAAAATAGGCCTAAgcg                                               <  2:160159/251‑1 (MQ=255)                                                                                                                                                                                                      aGACCTGGTTCTCACTGTAACCCAAATGCTGCGCAAACATGGAGTGGAGAGGAAATTCGTCGAATTTTATGGTGCTGGTCTGGATTCACTACCGTTGGAGGACAGCGCCACCATTGCACATATGTCGCCAGAATATGGGGCAACCTGAGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACTTGCGTTTAAGCGGGCGCAGCGAAGATCGGGTCGGGTTGGTCGAAAAATATGCCAAAGCGCAggga                                         >  2:280515/1‑249 (MQ=255)                                                                                                                                                                                                       gACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAt                                       >  1:148780/1‑251 (MQ=255)                                                                                                                                                                                                            ggTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtgtg                                    <  2:321328/249‑1 (MQ=255)                                                                                                                                                                                                                        tACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAAGTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGc                       <  2:307604/250‑1 (MQ=255)                                                                                                                                                                                                                              aaTGCTGCGCAAACATGTCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAc                <  1:408371/251‑1 (MQ=255)                                                                                                                                                                                                                                  ctgcGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTATGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCAAt            >  1:225102/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCCCCACCATTGCCAATATGTCGCCAGACTATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCACTTTTTACCAGt  >  1:416800/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                AATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGG‑TGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCAATTTTTACCAGT  >  NC_000913/1336429‑1336913

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