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%	107.3 / NA	35	S230S (AGC→AGT)	alsA	D‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/15);  total (20/15)

TCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAT  >  NC_000913/4310157‑4310422                                                                                                                                       |                                                                                                                                    tCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTa                                                                                                                                 >  1:625470/1‑139 (MQ=255)  cACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTAc                                                                                                                                >  1:158658/1‑139 (MQ=255)                        cTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGctgctg                                                                                                          >  1:127033/1‑139 (MQ=255)                         tCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                         <  2:28251/139‑1 (MQ=255)                                  gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                >  2:692416/1‑139 (MQ=255)                                  gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                >  1:637595/1‑139 (MQ=255)                                     tCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcatca                                                                                             <  1:527897/139‑1 (MQ=255)                                               ggTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                               >  2:523236/1‑127 (MQ=255)                                               ggTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                               <  1:523236/127‑1 (MQ=255)                                                   gCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCg                                                                               >  1:776339/1‑139 (MQ=255)                                                         aTTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATa                                                                         >  1:157047/1‑139 (MQ=255)                                                              ccTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACg                                                                                                                <  2:335820/95‑1 (MQ=255)                                                              ccTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACg                                                                                                                >  1:335820/1‑95 (MQ=255)                                                                     cgcgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                             >  1:490811/1‑139 (MQ=255)                                                                       cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                           <  1:210277/139‑1 (MQ=255)                                                                       cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                           >  1:358808/1‑139 (MQ=255)                                                                       cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                           >  2:209628/1‑139 (MQ=255)                                                                       cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATCTCCGc                                                           >  1:538974/1‑139 (MQ=255)                                                                                            cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                               >  1:257166/1‑82 (MQ=255)                                                                                            cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                               <  2:257166/82‑1 (MQ=255)                                                                                            cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                                      <  2:127033/139‑1 (MQ=255)                                                                                            cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                                      >  2:558834/1‑139 (MQ=255)                                                                                                    aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                    >  1:718100/1‑133 (MQ=255)                                                                                                    aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                    <  2:718100/133‑1 (MQ=255)                                                                                                       aTCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTAc                           <  1:85447/139‑1 (MQ=255)                                                                                                       aTCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTAc                           <  2:538974/139‑1 (MQ=255)                                                                                                           gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                    >  1:620813/1‑126 (MQ=255)                                                                                                           gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                    <  2:620813/126‑1 (MQ=255)                                                                                                           gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCctc                       >  1:473506/1‑139 (MQ=255)                                                                                                               gacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTa                   >  2:289569/1‑139 (MQ=255)                                                                                                                 cgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACg                 <  1:209628/139‑1 (MQ=255)                                                                                                                                acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt  <  2:158658/139‑1 (MQ=255)                                                                                                                                acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt  <  2:625470/139‑1 (MQ=255)                                                                                                                                acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt  >  2:703547/1‑139 (MQ=255)                                                                                                                                acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt  <  2:776339/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    TCACGACTGGTGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAT  >  NC_000913/4310157‑4310422

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