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%	141.7 / NA	40	I167I (ATT→ATC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/18);  total (22/18)

CTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACTG  >  NC_000913/4367749‑4367904                                                                               |                                                                              cTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGAt                                                                             >  1:830422/1‑81 (MQ=255)    cAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATc                                                                            <  1:830370/79‑1 (MQ=255)      gCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCAt                                                                        >  2:830422/1‑81 (MQ=255)         gggTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTa                                                                     >  1:830424/1‑81 (MQ=255)          ggTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTAc                                                                    <  1:830372/81‑1 (MQ=255)          ggTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTAc                                                                    <  2:830370/81‑1 (MQ=255)              cGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAg                                                                 <  1:830373/80‑1 (MQ=255)               gACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCt                                                               >  1:830419/1‑81 (MQ=255)                aCCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCt                                                               >  1:830414/1‑80 (MQ=255)                   cATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAAt                                                           >  1:830438/1‑81 (MQ=255)                      tttCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAg                                                        >  2:830414/1‑81 (MQ=255)                          aGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGa                                                     <  2:830373/80‑1 (MQ=255)                           ggATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAggg                                                      >  1:830429/1‑78 (MQ=255)                           ggATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAag                                                   <  2:830374/81‑1 (MQ=255)                             aTGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAagag                                                 >  2:830440/1‑81 (MQ=255)                                tCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTa                                               <  1:830376/80‑1 (MQ=255)                                tCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTa                                               <  2:830376/80‑1 (MQ=255)                                      aaTTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGc                                         >  2:830420/1‑80 (MQ=255)                                         tCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATa                                     >  1:830450/1‑81 (MQ=255)                                           gACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACg                                   >  1:830431/1‑81 (MQ=255)                                                cTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAAc                              >  1:830449/1‑81 (MQ=255)                                                cTTTACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAAc                              <  1:830378/81‑1 (MQ=255)                                                  ttACGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCg                            >  2:830425/1‑81 (MQ=255)                                                     cGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGt                          <  2:830378/80‑1 (MQ=255)                                                     cGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTc                         >  2:830441/1‑81 (MQ=255)                                                     cGTTCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTc                         >  2:830438/1‑81 (MQ=255)                                                       ttCAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCgg                       >  1:830427/1‑81 (MQ=255)                                                         cAAAGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCggg                      >  2:830435/1‑80 (MQ=255)                                                          aaaGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGt                     <  1:830380/80‑1 (MQ=255)                                                            aGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAgg                  >  1:830446/1‑81 (MQ=255)                                                            aGCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAgg                  >  1:830425/1‑81 (MQ=255)                                                             gCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGc                 <  1:830381/81‑1 (MQ=255)                                                             gCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGc                 <  1:830382/81‑1 (MQ=255)                                                             gCCTTCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGc                 <  1:830383/81‑1 (MQ=255)                                                                 tCACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgtcg             >  1:830448/1‑81 (MQ=255)                                                                  cACGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGcgtc              <  2:830383/79‑1 (MQ=255)                                                                   aCGCAGTTGGTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGtt           <  1:830385/81‑1 (MQ=255)                                                                            gTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACt   >  2:830433/1‑80 (MQ=255)                                                                            gTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACt   <  1:830386/80‑1 (MQ=255)                                                                            gTTGATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACTg  <  1:830387/81‑1 (MQ=255)                                                                               |                                                                              CTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGGGAAGAGTAAACTGCGATACGGAAACCGGTCGGGTAGGCGTCGTTAGTGGACTG  >  NC_000913/4367749‑4367904

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