Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,415	G→T	E228* (GAA→TAA)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,415	0	G	T	100.0%	143.6 / NA	41	E228* (GAA→TAA)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/20);  total (21/20)

CGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCG  >  NC_000913/4282‑4549                                                                                                                                      |                                                                                                                                       cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTg                                                                                                                                  >  1:456036/1‑140 (MQ=255)     gATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTgc                                                                                                                              <  1:844540/140‑1 (MQ=255)      aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTgcg                                                                                                                             <  2:784838/140‑1 (MQ=255)      aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTgcg                                                                                                                             <  1:951415/140‑1 (MQ=255)         tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAg                                                                                                                          <  2:172292/140‑1 (MQ=255)         tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAg                                                                                                                          <  1:609491/140‑1 (MQ=255)                       gcgcTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGAc                                                                                                            <  1:595140/140‑1 (MQ=255)                             gTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCaa                                                                                                      >  1:316254/1‑140 (MQ=255)                              tGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAAc                                                                                                     <  2:456036/140‑1 (MQ=255)                                     ggCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTgg                                                                                              >  1:342534/1‑140 (MQ=255)                                       cGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGtt                                                                                            >  1:591841/1‑140 (MQ=255)                                       cGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGtt                                                                                            <  2:720143/140‑1 (MQ=255)                                           tgatgatgaAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGtct                                                                                        <  2:904110/140‑1 (MQ=255)                                             atgatgaAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCg                                                                                      <  1:33259/140‑1 (MQ=255)                                                 tgaAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGc                                                                                  <  1:46440/140‑1 (MQ=255)                                                 tgaAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCAGTGc                                                                                  <  2:332368/140‑1 (MQ=255)                                                   aagaaCTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCa                                                                                >  2:49467/1‑140 (MQ=255)                                                                ggCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCg                                                                   <  1:49467/140‑1 (MQ=255)                                                                ggCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCg                                                                   >  2:211210/1‑140 (MQ=255)                                                                  cgcTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGc                                                                 >  2:350017/1‑140 (MQ=255)                                                                                 cGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggg                                                  >  1:292730/1‑139 (MQ=255)                                                                                 cGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggg                                                  >  1:286918/1‑139 (MQ=255)                                                                                   gCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                >  1:853135/1‑140 (MQ=255)                                                                                   gCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                >  2:893474/1‑140 (MQ=255)                                                                                   gCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                >  1:184535/1‑140 (MQ=255)                                                                                          cGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCg                                         <  2:606668/140‑1 (MQ=255)                                                                                            aTTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGaa                                       >  1:317690/1‑140 (MQ=255)                                                                                                  aTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACt                                 <  1:893474/140‑1 (MQ=255)                                                                                                      gCCGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGt                             <  2:342534/140‑1 (MQ=255)                                                                                                       ccGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCg                                         >  2:838749/1‑127 (MQ=255)                                                                                                       ccGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCg                                         <  1:838749/127‑1 (MQ=255)                                                                                                       ccGTTTGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTc                            >  1:99288/1‑140 (MQ=255)                                                                                                          tttGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGc                         >  1:885017/1‑140 (MQ=255)                                                                                                           ttGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGGAGTCACTCGGGCTGcc                        >  2:286979/1‑140 (MQ=255)                                                                                                            tGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGTCGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCg                       >  2:360508/1‑140 (MQ=255)                                                                                                            tGCTGGCGCAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACCCGGGCTGCCg                       >  2:92194/1‑140 (MQ=255)                                                                                                                  cgcAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGaaa                 <  1:674178/140‑1 (MQ=255)                                                                                                                    cAGATTTGCTACTACTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACg               <  2:853135/140‑1 (MQ=255)                                                                                                                            ctactaCTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTAttt       >  1:456429/1‑139 (MQ=255)                                                                                                                             tactacTTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATtgc      <  1:211210/140‑1 (MQ=255)                                                                                                                                 actTTTAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGGCTGCCGGGGAAAAGGTTTATTGCTGCg  >  2:314705/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       CGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCG  >  NC_000913/4282‑4549

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: ‑