Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,527,562	A→G	L140P (CTG→CCG)	sgcA ←	putative phosphotransferase enzyme IIA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,527,562	0	A	G	100.0%	182.7 / NA	50	L140P (CTG→CCG)	sgcA	putative phosphotransferase enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

TGCAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTGGCGCTGACGCAGAGAAGCAG  >  NC_000913/4527446‑4527695                                                                                                                     |                                                                                                                                      tGCAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCtgt                                                                                                                 >  1:632073/1‑139 (MQ=255)   cAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGc                                                                                                               >  1:726993/1‑139 (MQ=255)     ggCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTg                                                                                                             >  1:30269/1‑139 (MQ=255)               aTCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCa                                                                                                   <  2:30269/139‑1 (MQ=255)                  cAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTg                                                                                                <  1:572761/139‑1 (MQ=255)                         tCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAAc                                                                                                                               <  1:56301/101‑1 (MQ=255)                         tCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAAc                                                                                                                               >  2:56301/1‑101 (MQ=255)                                ccGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTg                                                                                                >  2:454504/1‑125 (MQ=255)                                ccGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTg                                                                                                <  1:454504/125‑1 (MQ=255)                                  gTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTg                                                                                                <  1:609418/123‑1 (MQ=255)                                  gTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTg                                                                                                >  2:609418/1‑123 (MQ=255)                                   tCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCaaa                                                                                                                                >  2:1046064/1‑90 (MQ=255)                                   tCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCaaa                                                                                                                                <  1:1046064/90‑1 (MQ=255)                                   tCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAg                                                                                                                             >  1:1404959/1‑93 (MQ=255)                                   tCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAg                                                                                                                             <  2:1404959/93‑1 (MQ=255)                                   tCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGc                                                                               >  2:1357891/1‑139 (MQ=255)                                    cAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCg                                                                              <  2:726993/139‑1 (MQ=255)                                    cAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCg                                                                              <  2:400756/139‑1 (MQ=255)                                      gCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGaa                                                                            <  1:1192468/139‑1 (MQ=255)                                                        gagGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATaaa                                                                    >  2:788860/1‑129 (MQ=255)                                                        gagGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATaaa                                                                    <  1:788860/129‑1 (MQ=255)                                                                   cTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGt                                               <  1:1441146/139‑1 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCtgtg                                                                                                                >  2:928046/1‑61 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCtgtg                                                                                                                <  1:928046/61‑1 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCg                                                                                      <  2:63742/87‑1 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCg                                                                                      >  1:63742/1‑87 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATaa                                                                     <  1:1098733/104‑1 (MQ=255)                                                                                ggTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATaa                                                                     >  2:1098733/1‑104 (MQ=255)                                                                                   gCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATgcg                                                        >  1:254275/1‑114 (MQ=255)                                                                                   gCAGGATCATAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATgcg                                                        <  2:254275/114‑1 (MQ=255)                                                                                            tAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGc                                                                                                                          >  2:125160/1‑39 (MQ=255)                                                                                            tAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGc                                                                                                                          <  1:125160/39‑1 (MQ=255)                                                                                            tAAGCGCCTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTg                      >  2:794428/1‑139 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCg                                                                                      <  2:351015/69‑1 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCg                                                                                      >  1:351015/1‑69 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   <  2:200221/104‑1 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   >  1:200221/1‑104 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGTATCCGTGgcg                   >  2:1380682/1‑136 (MQ=255)                                                                                                  ccTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGTATCCGTGgcg                   <  1:1380682/136‑1 (MQ=255)                                                                                                   cTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTg                                                             >  2:807445/1‑93 (MQ=255)                                                                                                   cTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   <  1:554446/103‑1 (MQ=255)                                                                                                   cTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   >  2:554446/1‑103 (MQ=255)                                                                                                   cTCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCGGGCAGTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTg                                                             <  1:807445/93‑1 (MQ=255)                                                                                                    tCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   <  1:657107/102‑1 (MQ=255)                                                                                                    tCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   >  2:657107/1‑102 (MQ=255)                                                                                                    tCCGTTAACCTGATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTGGCGCTGAc              >  1:705154/1‑139 (MQ=255)                                                                                                               gATACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTGGCGCTGACGCAGAGAagca   >  2:1019340/1‑139 (MQ=255)                                                                                                                aTACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   >  1:1142913/1‑90 (MQ=255)                                                                                                                aTACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGaa                                                   <  2:1142913/90‑1 (MQ=255)                                                                                                                aTACCGGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTGGCGCTGACGCAGAGAagcag  <  2:632073/139‑1 (MQ=255)                                                                                                                     |                                                                                                                                      TGCAGGCTCCCAAAATCCAGATTGTCCAATGCCGTCAGCTCACGACCATAGTGCAGAGGATTTGCGCTGGCGAGAGAAGGGTGCAGGATCATAAGCGCCTCCGTTAACCTGATACCAGCGCAAACAGTGCATCGTCTGTGCTGGCATTTCCCACTGCCGTTAAGCGTTGTGGCGAATCAATAAACTGGCTGATGCGCTGAATGGTCAAAATATGTGCGTTGGCATCCGTGGCGCTGACGCAGAGAAGCAG  >  NC_000913/4527446‑4527695

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