Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,825,986	C→T	R313H (CGC→CAC)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,825,986	0	C	T	100.0%	176.2 / NA	51	R313H (CGC→CAC)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/24);  total (27/24)

TACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACT  >  NC_000913/1825852‑1826116                                                                                                                                       |                                                                                                                                   taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAAt                                                                                                                               >  2:147594/1‑140 (MQ=255) taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAAt                                                                                                                               >  2:918312/1‑140 (MQ=255) taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAAt                                                                                                                               >  2:226221/1‑140 (MQ=255)      gATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCaa                                                                                                                          >  1:873672/1‑140 (MQ=255)        tGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACg                                                                                                                        <  1:147594/140‑1 (MQ=255)         gACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGc                                                                                                                       >  1:236192/1‑140 (MQ=255)                      taaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTa                                                                                                          <  2:260304/140‑1 (MQ=255)                        aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAgg                                                                                                        >  1:873936/1‑140 (MQ=255)                         aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTAc                                                                                                                    <  2:797735/127‑1 (MQ=255)                         aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTAc                                                                                                                    >  1:797735/1‑127 (MQ=255)                         aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAggg                                                                                                       >  1:765418/1‑140 (MQ=255)                         aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAggg                                                                                                       >  1:711076/1‑140 (MQ=255)                          gCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGa                                                                                                      <  1:226221/140‑1 (MQ=255)                          gCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGa                                                                                                      <  2:641238/140‑1 (MQ=255)                            aaaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                                                    <  2:873672/140‑1 (MQ=255)                            aaaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAACCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                                                    <  2:862305/140‑1 (MQ=255)                              aTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                                                  >  2:49283/1‑140 (MQ=255)                                aaCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTGTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                                                    >  2:632794/1‑136 (MQ=255)                                aaCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                                                    <  1:632794/136‑1 (MQ=255)                                  ccAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              <  2:873936/140‑1 (MQ=255)                                  ccAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              <  2:731683/140‑1 (MQ=255)                                    aGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                                            <  2:765418/140‑1 (MQ=255)                                    aGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                                            <  1:49283/140‑1 (MQ=255)                                        ttAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActct                                                                                        <  1:913847/140‑1 (MQ=255)                                           aTCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCATGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              <  1:862275/131‑1 (MQ=255)                                           aTCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCATGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              >  2:862275/1‑131 (MQ=255)                                            tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcat                                                                                    >  2:870805/1‑140 (MQ=255)                                              tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatca                                                                                  <  2:711076/140‑1 (MQ=255)                                                tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              <  1:537548/126‑1 (MQ=255)                                                tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                                              >  2:537548/1‑126 (MQ=255)                                                tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATg                                                                                >  2:532076/1‑140 (MQ=255)                                                      tttATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTAc                                                                          >  1:8673/1‑140 (MQ=255)                                                          taAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTa                                                                      <  1:532076/140‑1 (MQ=255)                                                                        aaTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTcc                                                        >  1:341472/1‑140 (MQ=255)                                                                               gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                                                 >  1:568223/1‑140 (MQ=255)                                                                               gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                                                 >  2:858507/1‑140 (MQ=255)                                                                               gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                                                 >  1:633498/1‑140 (MQ=255)                                                                               gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAAT‑CCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcg                                                >  2:694587/1‑140 (MQ=255)                                                                                cAGAATGATTTCTACGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcg                                                >  1:859737/1‑140 (MQ=255)                                                                                   aaTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAAT‑CCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAg                                            <  1:694587/140‑1 (MQ=255)                                                                                      gATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa                                          <  2:568223/140‑1 (MQ=255)                                                                                        tttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAAt                                        >  1:380002/1‑140 (MQ=255)                                                                                         ttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg                                       <  2:633498/140‑1 (MQ=255)                                                                                                     gTATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCaa                           <  2:348350/140‑1 (MQ=255)                                                                                                      tATTAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCaaa                          <  1:918312/140‑1 (MQ=255)                                                                                                         tAGCTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACgg                       >  2:569551/1‑140 (MQ=255)                                                                                                            cTTATTTTTTCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCAt                    >  1:279820/1‑140 (MQ=255)                                                                                                                     tCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATaaa                 <  2:456726/134‑1 (MQ=255)                                                                                                                     tCGAGCATTAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATaaa                 >  1:456726/1‑134 (MQ=255)                                                                                                                             tAATCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCAc   >  1:703793/1‑140 (MQ=255)                                                                                                                              aaTCCCGCGTGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACt  <  2:8673/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   TACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACT  >  NC_000913/1825852‑1826116

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