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

ACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCA  >  NC_000913/4310169‑4310424                                                                                                                           |                                                                                                                                      aCGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGc                                                                                                                       >  2:1353967/1‑139 (MQ=255)   gTTCCGCAGCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTCTGCCGCTGCaa                                                                                                                     >  2:1046764/1‑139 (MQ=255)   gTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCaa                                                                                                                     >  1:272386/1‑139 (MQ=255)     tCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAAc                                                                                                                   >  1:1055709/1‑139 (MQ=255)        gCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACgct                                                                                                                >  2:567155/1‑139 (MQ=255)         cACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACgctg                                                                                                               >  2:382874/1‑139 (MQ=255)             tCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                           <  1:283123/139‑1 (MQ=255)              cAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGcc                                                                                                          >  2:350987/1‑139 (MQ=255)                      gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                  >  1:1445543/1‑139 (MQ=255)                      gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                  <  2:1055709/139‑1 (MQ=255)                      gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                  >  1:595625/1‑139 (MQ=255)                      gTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                  <  2:146534/139‑1 (MQ=255)                       tttCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                 <  2:208107/139‑1 (MQ=255)                       tttCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                 <  1:820612/139‑1 (MQ=255)                       tttCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                 <  2:1303448/139‑1 (MQ=255)                               gcAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGtat                                                                                         >  1:122197/1‑139 (MQ=255)                                   ggTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCg                                                                                     <  2:417013/139‑1 (MQ=255)                                        cTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGc                                                                                <  2:122197/139‑1 (MQ=255)                                         tGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCCGCTGCCGTCTTTCATCACCGTATAGCGGTCGCa                                                                               >  2:346452/1‑139 (MQ=255)                                             aTTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCaaa                                                                                                                    >  2:701591/1‑98 (MQ=255)                                             aTTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCaaa                                                                                                                    <  1:701591/98‑1 (MQ=255)                                              ttCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCttt                                                                                                    >  1:434912/1‑113 (MQ=255)                                              ttCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCttt                                                                                                    <  2:434912/113‑1 (MQ=255)                                                     tCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAAt                                                                   <  1:567155/139‑1 (MQ=255)                                                         cgcgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                               >  1:323506/1‑139 (MQ=255)                                                           cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                             <  2:595625/139‑1 (MQ=255)                                                           cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                             <  2:1081307/139‑1 (MQ=255)                                                           cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                             <  1:382874/139‑1 (MQ=255)                                                                     gTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                                   <  2:424589/139‑1 (MQ=255)                                                                     gTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACCCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                                   <  1:1046764/139‑1 (MQ=255)                                                                      ttttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                  <  2:319824/139‑1 (MQ=255)                                                                           cAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGata                                             <  1:1308164/139‑1 (MQ=255)                                                                                  cgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                      <  2:819534/139‑1 (MQ=255)                                                                                  cgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                      <  1:1353967/139‑1 (MQ=255)                                                                                  cgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                      <  2:272386/139‑1 (MQ=255)                                                                                           aTCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTAc                             <  1:346452/139‑1 (MQ=255)                                                                                               gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                                   <  1:208599/113‑1 (MQ=255)                                                                                               gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                                   >  2:208599/1‑113 (MQ=255)                                                                                                   gacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCttt                                                                                                    <  1:980359/60‑1 (MQ=255)                                                                                                   gacgATATCGTCATATGACACATCACTTACTATGCCGCTGAAAACGCTGCTGCCGTCttt                                                                                                    >  2:980359/1‑60 (MQ=255)                                                                                                     cgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACg                   <  1:47701/139‑1 (MQ=255)                                                                                                      gATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                  >  1:637506/1‑107 (MQ=255)                                                                                                      gATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                  <  2:637506/107‑1 (MQ=255)                                                                                                             tCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGTGAGATATAGACGATGGCCGTACCCTCTTTACGCaa                >  2:1500937/1‑134 (MQ=255)                                                                                                             tCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGTGAGATATAGACGATGGCCGTACCCTCTTTACGCaa                <  1:1500937/134‑1 (MQ=255)                                                                                                             tCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGAt           <  2:1255657/139‑1 (MQ=255)                                                                                                              cATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACt              <  2:1116180/135‑1 (MQ=255)                                                                                                              cATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACt              >  1:1116180/1‑135 (MQ=255)                                                                                                                tttGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCa        >  2:576709/1‑139 (MQ=255)                                                                                                                  tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATg      >  2:1254975/1‑139 (MQ=255)                                                                                                                   gACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                 >  1:446739/1‑47 (MQ=255)                                                                                                                   gACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                 <  2:446739/47‑1 (MQ=255)                                                                                                                    acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt    <  1:804455/139‑1 (MQ=255)                                                                                                                    acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt    >  2:804455/1‑139 (MQ=255)                                                                                                                    acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt    >  1:1208879/1‑139 (MQ=255)                                                                                                                      acaTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCa  <  1:350987/139‑1 (MQ=255)                                                                                                                        aTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                               >  1:225016/1‑76 (MQ=255)                                                                                                                        aTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                               <  2:225016/76‑1 (MQ=255)                                                                                                                           |                                                                                                                                      ACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCA  >  NC_000913/4310169‑4310424

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