Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,826,277	G→A	P216L (CCC→CTC)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,277	0	G	A	100.0%	78.8 / NA	24	P216L (CCC→CTC)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

AGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGC  >  NC_000913/1826143‑1826389                                                                                                                                       |                                                                                                                 agtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAAATGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAAAGCAAc                                                                                                             >  1:515988/1‑140 (MQ=255)            aTTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCTAAg                                                                                                  <  2:490538/140‑1 (MQ=255)               tGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGcc                                                                                                     <  2:674104/134‑1 (MQ=255)               tGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGcc                                                                                                     >  1:674104/1‑134 (MQ=255)                 cGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCa                                                                                             <  1:35943/140‑1 (MQ=255)                 cGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCa                                                                                             <  2:349038/140‑1 (MQ=255)                      ccaccCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTa                                                                                        <  2:433811/140‑1 (MQ=255)                        accCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACCCTCACCAGAAAACAGCGCAGCAAATTCGCTGGCAGTTACTTCAAACTGGCGAAAATCGTTTTTCCCAAAAAACAAAGCCTTGCCAAATTCCACCGGACa                                                                                      >  2:35943/1‑140 (MQ=255)                                                     gATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGa                                                         <  2:48547/140‑1 (MQ=255)                                                      aTACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGaa                                                        >  1:544830/1‑140 (MQ=255)                                                      aTACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGaa                                                        >  1:593643/1‑140 (MQ=255)                                                      aTACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGaa                                                        >  2:265715/1‑140 (MQ=255)                                                        aCCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAAt                                                      >  1:696423/1‑140 (MQ=255)                                                          cGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGgcg                                                             <  1:208953/131‑1 (MQ=255)                                                          cGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGgcg                                                             >  2:208953/1‑131 (MQ=255)                                                           gACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATggg                                                   >  2:802367/1‑140 (MQ=255)                                                                 ctcACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACg                                             >  2:356277/1‑140 (MQ=255)                                                                        gATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTaccacc                                      <  2:696423/140‑1 (MQ=255)                                                                           agcagcGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAgg                                   >  1:413207/1‑140 (MQ=255)                                                                                cgcagcAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCt                              >  1:370645/1‑140 (MQ=255)                                                                                  cagcaATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGTCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGa                            >  1:486049/1‑140 (MQ=255)                                                                                     caATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAgg                                   >  1:585933/1‑130 (MQ=255)                                                                                     caATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAgg                                   <  2:585933/130‑1 (MQ=255)                                                                                            gcTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCa                  <  2:593643/140‑1 (MQ=255)                                                                                              tGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGAGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAgc                <  1:265715/140‑1 (MQ=255)                                                                                                            aaCTGGCGAAGATCGTTTTCCCAAAAGAGCACCGTTTTGCCAAGTTCCAGCGTACAGGTCAGCGCGCAAAAATGTCTGGCGCTAAAATGGGTAACCGTCCCCCCAGGTTCCTGATGGAACACCAGCCCCTCCAGCCccgc  <  2:486049/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                 AGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGC  >  NC_000913/1826143‑1826389

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