Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,367,827	A→G	I167I (ATT→ATC)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,367,827	0	A	G	100.0%	151.2 / NA	43	I167I (ATT→ATC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/21);  total (22/21)

CTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACT  >  NC_000913/4367749‑4367903                                                                               |                                                                             cTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGAt                                                                            <  1:898924/81‑1 (MQ=255)    cAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATc                                                                           <  1:898925/79‑1 (MQ=255)      gCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCAt                                                                       >  1:898981/1‑81 (MQ=255)          ggTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCt                                                                     <  1:898926/79‑1 (MQ=255)          ggTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCt                                                                     <  1:898928/79‑1 (MQ=255)          ggTACGACCCATTTTCAGGATGTCCTGGAATACGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTAc                                                                   <  2:898927/81‑1 (MQ=255)           gTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAATGCCTTCACGCAGTTGGTTGATCGCATCTAc                                                                   >  2:899010/1‑80 (MQ=255)           gTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTAcc                                                                  <  2:898929/81‑1 (MQ=255)           gTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTAc                                                                   >  1:898961/1‑80 (MQ=255)             aCGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAg                                                                >  2:898974/1‑81 (MQ=255)              cGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGc                                                               <  2:898930/81‑1 (MQ=255)               gACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCt                                                              >  2:898995/1‑81 (MQ=255)                 cccATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTa                                                            <  1:898931/81‑1 (MQ=255)                       ttCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAg                                                       <  2:898932/80‑1 (MQ=255)                       ttCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAg                                                       <  2:898933/80‑1 (MQ=255)                        tCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAgg                                                      >  2:898962/1‑80 (MQ=255)                         cAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAggg                                                     <  2:898934/80‑1 (MQ=255)                              tGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGt                                               <  1:898936/81‑1 (MQ=255)                              tGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGt                                               <  1:898935/81‑1 (MQ=255)                                 ccTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTaaa                                            >  2:898988/1‑81 (MQ=255)                                     gAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGc                                        <  1:898937/81‑1 (MQ=255)                                          cGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATAc                                   >  2:898978/1‑81 (MQ=255)                                          cGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATAc                                   >  1:898965/1‑81 (MQ=255)                                            aCAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACg                                  >  1:898982/1‑80 (MQ=255)                                            aCAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACg                                  >  1:898979/1‑80 (MQ=255)                                                  ttACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCg                           >  2:898989/1‑81 (MQ=255)                                                   tACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAAcc                            <  2:898938/79‑1 (MQ=255)                                                    aCGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGt                         >  2:898973/1‑81 (MQ=255)                                                     cGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTc                        <  1:898939/81‑1 (MQ=255)                                                       ttCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCgg                      >  2:899001/1‑81 (MQ=255)                                                        tCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCggg                     >  2:898966/1‑81 (MQ=255)                                                         cAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCggg                     <  2:898940/80‑1 (MQ=255)                                                          aaaGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCggg                     >  1:898987/1‑79 (MQ=255)                                                          aaaGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGt                    <  1:898941/80‑1 (MQ=255)                                                               cTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgt              <  1:898942/81‑1 (MQ=255)                                                                ttCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgtc             >  1:898967/1‑81 (MQ=255)                                                                 tCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgtcg            >  1:899007/1‑81 (MQ=255)                                                                  cACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgtcgt           >  1:898985/1‑81 (MQ=255)                                                                   aCGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGtt          <  2:898943/81‑1 (MQ=255)                                                                     gCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAg        >  2:898997/1‑81 (MQ=255)                                                                       aGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTg      <  2:898944/81‑1 (MQ=255)                                                                           ggTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACt  >  1:898972/1‑81 (MQ=255)                                                                           ggTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACt  >  1:898980/1‑81 (MQ=255)                                                                               |                                                                             CTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACT  >  NC_000913/4367749‑4367903

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