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%	96.2 / NA	29	G282R (GGG→AGG)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/12);  total (17/12)

GGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCAATTTTTACCAGTACGT  >  NC_000913/1336425‑1336917                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    ggTGAATGGATTGATTATCCGGATACACTCGTTGGTACTGTCTCGCACACAACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGACGCCG‑CAATGTTAGGCCAGCCGGTTTCAAGGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAACATTCCGTGAAGTAATTCACACCACAGCCCTGGTTAGCACTGTTTCCCAAATGCTGCTCAAACAGGGCGTGGTggg                                                                                                                                                                                                                                                    >  2:140574/1‑251 (MQ=255)                 aTCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAAttt                                                                                                                                                                                                                                    >  2:183667/1‑251 (MQ=255)                 aTCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAAttt                                                                                                                                                                                                                                    >  2:106166/1‑251 (MQ=255)                               ttttgttctgactcgcccccccacaTGAACAACGGCCTTGGGGTGCTGGGGGGGGGGGTTGGGTGGATCGAAGCAAAAGCCGACAA‑GTTAGGCCCGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACCTACCGGAAAATTACGGGAAGGTTTTTCCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTc                                                                                                                                                                                                                       <  1:256828/227‑1 (MQ=255)                                                      catgacaaaCGGCCTTGGCGTGCGGGGGTGGGGCGTTGTTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGTTTTCCATGCTTACCCCGGATGTAGTGGGCTTCAAACTTACCGAAAAATTACGTGAAGTTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATc                                                                                                                                                                                               <  1:84551/244‑1 (MQ=255)                                                       aTGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGGTcg                                                                                                                                                                                              >  1:150782/1‑251 (MQ=255)                                                                  ccTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGTCTTAAAACTTACCGAAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTg                                                                                                                                                                                   <  1:318704/251‑1 (MQ=255)                                                                                          gggtggGATGAAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATTCTTATCCCGGATGTAGTGGGCTTCAAACTTCCCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTg                                                                                                                                                           <  1:198113/250‑1 (MQ=255)                                                                                                     aaGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGc                                                                                                                                                >  2:261943/1‑251 (MQ=255)                                                                                                                    aTGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCg                                                                                                                                   >  1:185491/1‑250 (MQ=255)                                                                                                                      gTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAAttt                                                                                                                                                                                                                                    >  1:203732/1‑151 (MQ=255)                                                                                                                        tAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCt                                                                                                                              >  1:256189/1‑251 (MQ=255)                                                                                                                             cAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAc                                                                                                                         >  1:286592/1‑251 (MQ=255)                                                                                                                                                  ccGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGCAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTa                                                                                                      >  2:245638/1‑249 (MQ=255)                                                                                                                                                       tGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAACCTCGATTACATGCGTTTAAGCGGgc                                                                                               >  1:292915/1‑251 (MQ=255)                                                                                                                                                                 ttCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGa                                                                                        >  1:15072/1‑248 (MQ=255)                                                                                                                                                                     aaCTTACCGAAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAg                                                                                 <  1:183667/251‑1 (MQ=255)                                                                                                                                                                         tACCGAAAAATTTCGGGAAGTTTTTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc                                                                              <  2:15072/250‑1 (MQ=255)                                                                                                                                                                            cGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGt                                                                          >  1:175529/1‑251 (MQ=255)                                                                                                                                                                             ggAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGtt                                                                         >  1:118544/1‑251 (MQ=255)                                                                                                                                                                                       gTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaaa                                                               <  1:2478/251‑1 (MQ=255)                                                                                                                                                                                           aGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCCATCGATGCTGTAACCCTCGCTTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATg                                                           >  1:34004/1‑251 (MQ=255)                                                                                                                                                                                                  aCCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAg                                                     <  2:286592/250‑1 (MQ=255)                                                                                                                                                                                                   acGCCACAGACCTGGTTCTCACTGATACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGCTGGTCTGGATTCAATACCGTTGGCGGATCGCGCCACCATTGCCAATATGTAGCCAGAATATGGGGCCACCTGTGGCTTCTTCCCACTCGATGCTGTAACCCTCGATTACATACGTTTAAGCGGGCGAAGCGGAGATCAGTGCGAGTTGGTCGAAAAATATGCCATAgcg                                                   >  2:345495/2‑251 (MQ=255)                                                                                                                                                                                                         caGACCTGGTTCTCACTTGTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGc                                             <  2:185491/251‑1 (MQ=255)                                                                                                                                                                                                             ccTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtgt                                         >  1:208510/1‑251 (MQ=255)                                                                                                                                                                                                                ggTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGAt                                                                                                                  <  2:212093/175‑1 (MQ=255)                                                                                                                                                                                                                              cccAAATGCTGCGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGAt                        <  1:245638/251‑1 (MQ=255)                                                                                                                                                                                                                                     gctgcGCAAACATGGCGTGGTGAGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGCTTACATGCGTTTAAGAGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCCAAGCGACGGGCATGTGGCGTAACCCGGGCGATGAACCa                  >  2:26818/1‑250 (MQ=255)                                                                                                                                                                                                                                                    cGTGGTGAGGAAAGTCTTTGAATGTTATGGGGATGGGCTGGATTCACTACCGTCGGTGGATCGCGCCACCATTGCCAATATGCCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCCGTAACCCACGAGTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGGCGAAAAATATGGCAAAGCGCATGGCATGGGGCGGAACGTGGGCGAGGAACCAATTTATACCAGTACGt  <  2:305307/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    GGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCG‑CAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGGAAAATTACGTGAAGGTATTACCGCCACAGACCTGGTTCTCACTGTTACCCAAATGCTGCGCAAACATGGCGTGGTGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAACCAATTTTTACCAGTACGT  >  NC_000913/1336425‑1336917

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: ‑