Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,826,546	T→C	E126E (GAA→GAG)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,546	0	T	C	100.0%	105.6 / NA	30	E126E (GAA→GAG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/14);  total (16/14)

GCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG  >  NC_000913/1826461‑1826671                                                                                      |                                                                                                                              gCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACg                                                                          >  1:753181/1‑139 (MQ=255)      gAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGc                                                                     <  1:232770/139‑1 (MQ=255)                 tGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGa                                                          <  2:1272972/139‑1 (MQ=255)                        tGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTcgcg                                                                                                                    <  1:417076/74‑1 (MQ=255)                        tGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTcgcg                                                                                                                    >  2:417076/1‑74 (MQ=255)                            gATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTTCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGTTTCATATCCCt                                               <  2:762437/139‑1 (MQ=255)                             aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                              <  1:508986/139‑1 (MQ=255)                             aTCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                              >  2:508986/1‑139 (MQ=255)                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                              >  2:760488/1‑134 (MQ=255)                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTa                                              <  1:760488/134‑1 (MQ=255)                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                         >  2:23085/1‑139 (MQ=255)                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                         >  2:860558/1‑139 (MQ=255)                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                         >  2:480589/1‑139 (MQ=255)                                    tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAAt                                       >  1:522065/1‑139 (MQ=255)                                    tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAAt                                       >  1:1084390/1‑139 (MQ=255)                                     gCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAc                                                                        >  2:1004502/1‑105 (MQ=255)                                     gCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAc                                                                        <  1:1004502/105‑1 (MQ=255)                                       cAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAAc                                    >  2:383771/1‑139 (MQ=255)                                          cgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAg                                                                                   >  1:1041967/1‑89 (MQ=255)                                          cgcACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAg                                                                                   <  2:1041967/89‑1 (MQ=255)                                                         cccTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGGtt                       <  2:762090/134‑1 (MQ=255)                                                         cccTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCtt                       >  1:762090/1‑134 (MQ=255)                                                            tGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCag               <  2:501567/139‑1 (MQ=255)                                                             ggTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGa                                                          <  2:133298/95‑1 (MQ=255)                                                             ggTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGa                                                          >  1:133298/1‑95 (MQ=255)                                                              gTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCagga             >  2:687087/1‑139 (MQ=255)                                                                   aaaaaTCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGc                                                                                         <  2:120268/58‑1 (MQ=255)                                                                   aaaaaTCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGc                                                                                         >  1:120268/1‑58 (MQ=255)                                                                         cTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg  <  2:753181/139‑1 (MQ=255)                                                                         cTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAg  <  1:480589/139‑1 (MQ=255)                                                                                      |                                                                                                                              GCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAG  >  NC_000913/1826461‑1826671

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