Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,310,291	G→A	S230S (AGC→AGT)	alsA ←	D‑allose ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,310,291	0	G	A	100.0%	104.4 / NA	33	S230S (AGC→AGT)	alsA	D‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

TCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAG  >  NC_000913/4310157‑4310425                                                                                                                                       |                                                                                                                                       tCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTa                                                                                                                                    >  2:766586/1‑139 (MQ=255) tCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTa                                                                                                                                    >  2:650171/1‑139 (MQ=255)   aCGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACt                                                                                                                                  >  2:532967/1‑139 (MQ=255)         ggTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCg                                                                                                                            >  1:657435/1‑139 (MQ=255)           tGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCt                                                                                                                          >  1:309832/1‑139 (MQ=255)           tGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCt                                                                                                                          >  1:659250/1‑139 (MQ=255)                      aCCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACgctgc                                                                                                               <  1:403626/139‑1 (MQ=255)                         tCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                            <  1:82863/139‑1 (MQ=255)                           aaaaaCCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTGGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCg                                                                                                          >  2:579315/1‑139 (MQ=255)                                ccGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCttt                                                                                                     <  1:67068/139‑1 (MQ=255)                                  gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATg                                                                                                                               >  1:56149/1‑111 (MQ=255)                                  gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATg                                                                                                                               <  2:56149/111‑1 (MQ=255)                                  gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                   <  1:450214/139‑1 (MQ=255)                                     tCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcatca                                                                                                <  2:122577/139‑1 (MQ=255)                                     tCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcatca                                                                                                >  1:369046/1‑139 (MQ=255)                                                           tctcCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                  <  2:692353/115‑1 (MQ=255)                                                           tctcCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                  >  1:692353/1‑115 (MQ=255)                                                              ccTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCg                                                                       >  2:196253/1‑139 (MQ=255)                                                                       cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                              >  2:763183/1‑139 (MQ=255)                                                                         ttAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCa                                                            >  1:683107/1‑139 (MQ=255)                                                                                ggTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                     <  2:120874/139‑1 (MQ=255)                                                                                  ttttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                   <  1:382390/139‑1 (MQ=255)                                                                                  ttttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                   <  1:650171/139‑1 (MQ=255)                                                                                  ttttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                   <  1:161223/139‑1 (MQ=255)                                                                                    ttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCgag                                                 <  1:178342/139‑1 (MQ=255)                                                                                              cgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                       <  1:651249/139‑1 (MQ=255)                                                                                                      cATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTa                               <  1:532967/139‑1 (MQ=255)                                                                                                               gacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTa                      >  1:11832/1‑139 (MQ=255)                                                                                                                 cgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACg                    <  1:766586/139‑1 (MQ=255)                                                                                                                              tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTa                      <  2:20695/124‑1 (MQ=255)                                                                                                                              tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTa                      >  1:20695/1‑124 (MQ=255)                                                                                                                              tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATg       >  2:571149/1‑139 (MQ=255)                                                                                                                                   caTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAg  >  1:368384/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       TCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAG  >  NC_000913/4310157‑4310425

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