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%	165.6 / NA	46	E126E (GAA→GAG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/20);  total (26/20)

GTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA  >  NC_000913/1826429‑1826673                                                                                                                      |                                                                                                                                gTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAg                                                                                                            >  1:483191/1‑139 (MQ=255)       ttGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTtct                                                                                                      >  1:1068744/1‑139 (MQ=255)             tAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCtt                                                                                                <  2:1068744/139‑1 (MQ=255)                     cGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCaaa                                                                                        >  1:1419446/1‑139 (MQ=255)                       aaCTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATa                                                                                      >  1:762982/1‑139 (MQ=255)                              gATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCa                                                                               >  2:174924/1‑139 (MQ=255)                                 gCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACg                                                                            >  1:993076/1‑139 (MQ=255)                                         ccACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCttt                                                                                               <  2:1192253/112‑1 (MQ=255)                                         ccACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCttt                                                                                               >  1:1192253/1‑112 (MQ=255)                                         ccACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCaa                                                                    <  2:172324/139‑1 (MQ=255)                                                      tATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                                                       <  2:762982/139‑1 (MQ=255)                                                      tATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                                                       <  1:343034/139‑1 (MQ=255)                                                           aGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGc                                                  >  1:900991/1‑139 (MQ=255)                                                                aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGa                                                            <  2:600374/124‑1 (MQ=255)                                                                aaGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGa                                                            >  1:600374/1‑124 (MQ=255)                                                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGc                                                                                                  <  1:981538/84‑1 (MQ=255)                                                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGc                                                                                                  >  2:981538/1‑84 (MQ=255)                                                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                           >  1:980493/1‑139 (MQ=255)                                                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                           >  1:1165106/1‑139 (MQ=255)                                                                  ggTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCa                                           >  2:959963/1‑139 (MQ=255)                                                                    tGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAAt                                         <  2:1018522/139‑1 (MQ=255)                                                                      ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                                                       >  1:1640049/1‑123 (MQ=255)                                                                      ccAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCat                                                       <  2:1640049/123‑1 (MQ=255)                                                                       cAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAAc                                      >  2:312811/1‑139 (MQ=255)                                                                       cAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAAc                                      >  1:456081/1‑139 (MQ=255)                                                                            cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCg                                                             >  1:1282588/1‑111 (MQ=255)                                                                            cacaGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCg                                                             <  2:1282588/111‑1 (MQ=255)                                                                              caGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTccc                               >  2:1220223/1‑139 (MQ=255)                                                                              caGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCT‑CAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCt                              >  2:573189/1‑139 (MQ=255)                                                                               aGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCt                              <  1:959963/139‑1 (MQ=255)                                                                                 aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGAtt                                                          <  2:698705/109‑1 (MQ=255)                                                                                 aTTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGAtt                                                          >  1:698705/1‑109 (MQ=255)                                                                                   tCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCT‑CAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCtt                         <  1:573189/139‑1 (MQ=255)                                                                                    cTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCtt                         <  2:900991/139‑1 (MQ=255)                                                                                       tGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAg                      >  1:997602/1‑139 (MQ=255)                                                                                       tGCCCTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAg                      >  2:566873/1‑139 (MQ=255)                                                                                         cccTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCata                                                      <  2:310197/105‑1 (MQ=255)                                                                                         cccTGGTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCata                                                      >  1:310197/1‑105 (MQ=255)                                                                                              gTCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCagga               >  2:406991/1‑139 (MQ=255)                                                                                               tCATAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCaggag              <  2:456081/139‑1 (MQ=255)                                                                                                 aTAAAAATCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAggagga            <  1:312811/139‑1 (MQ=255)                                                                                                     aaaTCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTcc        <  1:1220223/139‑1 (MQ=255)                                                                                                       aTCTTCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCg      <  2:980493/139‑1 (MQ=255)                                                                                                           tCCAGGCACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa  >  1:1164218/1‑139 (MQ=255)                                                                                                                 cACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAcc                                                                         <  2:619925/62‑1 (MQ=255)                                                                                                                 cACTGCTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGAcc                                                                         >  1:619925/1‑62 (MQ=255)                                                                                                                      |                                                                                                                                GTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATAAAAATCTTCCAGGCACTGTTCCAGTTCGCGCGCCCGACAGGTTTCTCCGCTTTCAGCAAATAGCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGA  >  NC_000913/1826429‑1826673

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