Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,405,459	G→T	R258L (CGT→CTT)	purA →	adenylosuccinate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,405,459	0	G	T	100.0%	139.2 / NA	40	R258L (CGT→CTT)	purA	adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/19);  total (21/19)

ATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAC  >  NC_000913/4405336‑4405591                                                                                                                            |                                                                                                                                     aTTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGtt                                                                                                                      <  1:659798/140‑1 (MQ=255)             ttGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTc                                                                                                          <  1:75798/140‑1 (MQ=255)                aaGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCaaa                                                                                                       >  1:729532/1‑140 (MQ=255)                   gTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCt                                                                                                    <  2:691800/140‑1 (MQ=255)                         aGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTcc                                                                                              <  2:473783/140‑1 (MQ=255)                            gTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACt                                                                                           <  2:715536/140‑1 (MQ=255)                                  tgctgGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCgtgt                                                                                      >  2:831588/1‑139 (MQ=255)                                  tgctgGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCgtgt                                                                                      <  1:831588/139‑1 (MQ=255)                                       gATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGc                                                                                <  1:705532/140‑1 (MQ=255)                                                cACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGtt                                                                       <  2:729532/140‑1 (MQ=255)                                                 aCGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTg                                                                                 <  1:395621/129‑1 (MQ=255)                                                 aCGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTg                                                                                 >  2:395621/1‑129 (MQ=255)                                                  cGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTcc                                                                     <  2:52017/140‑1 (MQ=255)                                                    gTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCccg                                                                   >  1:60235/1‑140 (MQ=255)                                                          aTCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGCGTAGGTGCAGGTCCGTTCCCGACCGaa                                                             >  2:365485/1‑140 (MQ=255)                                                              gtacgtaACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGt                                                         <  2:303096/140‑1 (MQ=255)                                                               tacgtaACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGtt                                                        >  2:694533/1‑140 (MQ=255)                                                               tacgtaACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGtt                                                        >  1:121640/1‑140 (MQ=255)                                                                    aaCTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTtgatg                                                   >  1:221324/1‑140 (MQ=255)                                                                     aCTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTtgatga                                                  >  1:649752/1‑140 (MQ=255)                                                                         cttcCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACt                                              >  1:408779/1‑140 (MQ=255)                                                                           tcCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTgg                                            >  2:206696/1‑140 (MQ=255)                                                                             cAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCg                                          <  1:694533/140‑1 (MQ=255)                                                                               acacCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAg                                        <  2:649752/140‑1 (MQ=255)                                                                                   cacTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTcc                                    >  2:692720/1‑140 (MQ=255)                                                                                        ctgGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGc                               <  2:408779/140‑1 (MQ=255)                                                                                           gtggcgtggcGACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAg                            >  1:336050/1‑140 (MQ=255)                                                                                             ggcgtggcgACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCa                          <  1:353425/140‑1 (MQ=255)                                                                                               cgtggcgACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAgg                        <  1:692720/140‑1 (MQ=255)                                                                                                gtggcgACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAggg                       >  1:609442/1‑140 (MQ=255)                                                                                                  ggcgACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTa                     <  2:221324/140‑1 (MQ=255)                                                                                                  ggcgACCGGTTCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTa                     >  2:434336/1‑140 (MQ=255)                                                                                                            tCCGGCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCgg           <  1:206696/140‑1 (MQ=255)                                                                                                                gCCTGGGCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCa       >  1:918138/1‑140 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTTCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAAc     >  2:343736/1‑137 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTTCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAAc     <  1:343736/137‑1 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAc  >  2:532090/1‑140 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAc  >  2:481044/1‑140 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAc  >  2:404271/1‑140 (MQ=255)                                                                                                                     ggCCCGCTTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAc  >  1:182225/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                     ATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTAC  >  NC_000913/4405336‑4405591

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