Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,156,086	G→A	S25N (AGC→AAC)	argB →	acetylglutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,156,086	0	G	A	100.0%	199.3 / NA	54	S25N (AGC→AAC)	argB	acetylglutamate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/29);  total (25/29)

GCAATGATGAATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAGCGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGATGAG  >  NC_000913/4156010‑4156162                                                                             |                                                                             gCAATGATGAATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCa                                                                          >  2:823431/1‑81 (MQ=255)       atgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGt                                                                     >  2:823430/1‑80 (MQ=255)       atgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGt                                                                     <  2:823380/80‑1 (MQ=255)       atgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGt                                                                     <  2:823379/80‑1 (MQ=255)       atgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTg                                                                    <  2:823381/81‑1 (MQ=255)        tgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGa                                                                   <  1:823383/81‑1 (MQ=255)        tgaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGa                                                                   <  1:823382/81‑1 (MQ=255)         gaATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGaa                                                                  >  1:823432/1‑81 (MQ=255)           aTCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAAt                                                                 >  1:823420/1‑80 (MQ=255)            tCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTa                                                               <  1:823385/81‑1 (MQ=255)            tCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTa                                                               <  2:823384/81‑1 (MQ=255)                ttaattaTCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGt                                                           >  1:823451/1‑81 (MQ=255)                  aattaTCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTg                                                          <  2:823387/80‑1 (MQ=255)                  aattaTCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTg                                                          <  1:823386/80‑1 (MQ=255)                   attaTCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAg                                                        <  2:823388/81‑1 (MQ=255)                      aTCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTca                                                     >  1:823389/1‑81 (MQ=255)                        cAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTcaca                                                   >  2:823439/1‑81 (MQ=255)                         aaaCTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACAt                                                  >  1:823449/1‑81 (MQ=255)                          aaCTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACAt                                                  >  2:823436/1‑80 (MQ=255)                               ggcggcGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGc                                              <  1:823390/79‑1 (MQ=255)                               ggcggcGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGc                                              <  2:823391/79‑1 (MQ=255)                                gcggcgTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTc                                           <  1:823392/81‑1 (MQ=255)                                 cggcgTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTcc                                          >  1:823442/1‑81 (MQ=255)                                  ggcgTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCg                                         >  2:823441/1‑81 (MQ=255)                                  ggcgTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCg                                         <  1:823393/81‑1 (MQ=255)                                    cgTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCt                                       >  1:823446/1‑81 (MQ=255)                                       aCTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGt                                    >  1:823444/1‑81 (MQ=255)                                       aCTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGt                                    >  2:823445/1‑81 (MQ=255)                                        ctgctgGATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGt                                    >  2:823435/1‑80 (MQ=255)                                              gATAGTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATtgtg                             <  2:823394/81‑1 (MQ=255)                                                  gTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcg                         <  2:823395/81‑1 (MQ=255)                                                  gTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcg                         <  1:823396/81‑1 (MQ=255)                                                  gTGAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAc                          <  1:823397/80‑1 (MQ=255)                                                    gAAGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggc                       >  1:823433/1‑81 (MQ=255)                                                     aaGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggc                       <  2:823398/80‑1 (MQ=255)                                                     aaGAGGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggc                       <  2:823399/80‑1 (MQ=255)                                                       gagGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggcggc                    >  2:823452/1‑81 (MQ=255)                                                       gagGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggcggc                    <  2:823400/81‑1 (MQ=255)                                                       gagGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggcggc                    >  2:823443/1‑81 (MQ=255)                                                        agGCGCTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCAcggcggcg                   >  2:823464/1‑81 (MQ=255)                                                            gcTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTg               >  2:823466/1‑81 (MQ=255)                                                            gcTGGAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTg               >  1:823448/1‑81 (MQ=255)                                                                gAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCgtg           <  1:823402/81‑1 (MQ=255)                                                                gAACGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCg             <  2:823401/79‑1 (MQ=255)                                                                   cGTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCgtggtg        <  1:823403/81‑1 (MQ=255)                                                                    gTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCgtggtg        >  1:823453/1‑80 (MQ=255)                                                                    gTCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGtggtgg       <  2:823404/81‑1 (MQ=255)                                                                     tCTGTTTAACGCACTGGTGAATTATCGTGGGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGa      >  1:823440/1‑81 (MQ=255)                                                                     tCTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGa      <  1:823405/81‑1 (MQ=255)                                                                      cTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGAt     >  2:823447/1‑81 (MQ=255)                                                                      cTGTTTAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGAt     >  1:823456/1‑81 (MQ=255)                                                                         tttAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGAt     <  2:823407/78‑1 (MQ=255)                                                                         tttAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGAt     <  1:823408/78‑1 (MQ=255)                                                                         tttAACGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGATGAg  <  2:823406/81‑1 (MQ=255)                                                                             |                                                                             GCAATGATGAATCCATTAATTATCAAACTGGGCGGCGTACTGCTGGATAGTGAAGAGGCGCTGGAACGTCTGTTTAGCGCACTGGTGAATTATCGTGAGTCACATCAGCGTCCGCTGGTGATTGTGCACGGCGGCGGTTGCGTGGTGGATGAG  >  NC_000913/4156010‑4156162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑