Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	29,918	A→G	S90G (AGC→GGC)	carA →	carbamoyl phosphate synthetase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,918	0	A	G	100.0%	109.4 / NA	32	S90G (AGC→GGC)	carA	carbamoyl phosphate synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/15);  total (17/15)

TTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATG  >  NC_000913/29794‑30038                                                                                                                             |                                                                                                                         ttCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTc                                                                                                                  >  2:71849/1‑133 (MQ=255)     cGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTa                                                                                                              <  2:304292/133‑1 (MQ=255)             cGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGaa                                                                                                      <  1:71849/133‑1 (MQ=255)                tACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGAc                                                                                                   <  2:45552/133‑1 (MQ=255)                  ctctTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAACACct                                                                                                 <  2:470919/133‑1 (MQ=255)                          taTCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTcttctt                                                                                         >  1:421161/1‑133 (MQ=255)                             ccccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTAcc                                                                                      >  1:12026/1‑133 (MQ=255)                             ccccATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAAGACCGAATACCTCTCTTCTTAcc                                                                                      >  1:249004/1‑133 (MQ=255)                             ccccATATTGGCAATGTCGGCACCAATGACGCCAATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTAcc                                                                                      >  1:322060/1‑133 (MQ=255)                                cATATTGGCATTGTCGGCACCAATTACGCCGAGGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGa                                                                                   <  2:205559/133‑1 (MQ=255)                                 atatTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGaa                                                                                  >  1:339592/1‑133 (MQ=255)                                 atatTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGaa                                                                                  >  2:350433/1‑133 (MQ=255)                                   atTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAAc                                                                                <  1:2554/133‑1 (MQ=255)                                   atTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAAc                                                                                >  2:83428/1‑133 (MQ=255)                                       gCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAAATTCCGTAAGACCGAAGACCTCTCTTGTTACTTGAAACGCCa                                                                            >  2:140876/1‑133 (MQ=255)                                         aaTGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATa                                                                          <  2:355471/133‑1 (MQ=255)                                             tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACAt                                                                      >  1:346876/1‑133 (MQ=255)                                                                  aggaaTCTTCGCAGGTACATGCAAGATGTCTGTTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGa                                                 <  2:249004/131‑1 (MQ=255)                                                                              aGGTCCATTCACAAGTTCCGGTTATTCGCCCCCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCt                                     <  1:140876/133‑1 (MQ=255)                                                                              aGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCt                                     <  2:455082/133‑1 (MQ=255)                                                                               ggTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTAGCGATTGACGATATAGATACCCGTAAGCTg                                    >  1:209272/1‑133 (MQ=255)                                                                                   cATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCAATACCCGTAGGCTGAcgc                                >  1:193720/1‑133 (MQ=255)                                                                                          aaGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACt                         >  1:65215/1‑133 (MQ=255)                                                                                            ggTCTGGTGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTgc                       >  1:421496/1‑133 (MQ=255)                                                                                                   tGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGaa                <  2:12026/133‑1 (MQ=255)                                                                                                   tGATTCGCGACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGaa                <  2:421496/133‑1 (MQ=255)                                                                                                        cgcgACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGgcg           <  2:339592/133‑1 (MQ=255)                                                                                                        cgcgACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGgcg           >  1:126763/1‑133 (MQ=255)                                                                                                        cgcgACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGgcg           >  1:104707/1‑133 (MQ=255)                                                                                                          cgACCTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGca         <  1:350433/133‑1 (MQ=255)                                                                                                              cTGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGa     <  1:83428/133‑1 (MQ=255)                                                                                                               tGCCGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGaa    >  2:114070/1‑133 (MQ=255)                                                                                                                 ccGCTGATTGCCGGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATg  >  1:124684/1‑133 (MQ=255)                                                                                                                             |                                                                                                                         TTCTCGTCAAATCGTTACTCTTACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATG  >  NC_000913/29794‑30038

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