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

TGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGA  >  NC_000913/4310167‑4310426                                                                                                                             |                                                                                                                                        tGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCt                                                                                                                           <  2:1322709/139‑1 (MQ=255)    cGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCa                                                                                                                        <  1:1496753/139‑1 (MQ=255)           cACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACgctg                                                                                                                 >  1:286183/1‑139 (MQ=255)               tCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                             >  1:1070665/1‑139 (MQ=255)                 aaaaaCCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCaaa                                                                                                                      <  2:1701/128‑1 (MQ=255)                 aaaaaCCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCaaa                                                                                                                      >  1:1701/1‑128 (MQ=255)                 aaaaaCCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCg                                                                                                           <  2:1070665/139‑1 (MQ=255)                         tttCGTGCGCAAGGTTTCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                   <  1:1126704/139‑1 (MQ=255)                                          cTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGc                                                                                  >  2:993961/1‑139 (MQ=255)                                          cTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGc                                                                                  >  2:1123333/1‑139 (MQ=255)                                                           cgcgTTAAAACGGTTTTGCAGTTCGCGGTCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcg                                                                                                    <  1:1019005/104‑2 (MQ=255)                                                           cgcgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtca                                                                                                    >  2:1019005/1‑104 (MQ=255)                                                           cgcgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                                 >  1:1600441/1‑139 (MQ=255)                                                             cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                   <  1:345223/103‑1 (MQ=255)                                                             cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTtcat                                                                                                   >  2:345223/1‑103 (MQ=255)                                                             cgTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                               <  1:1344616/139‑1 (MQ=255)                                                                      ggTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:280456/139‑1 (MQ=255)                                                                      ggTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:859177/139‑1 (MQ=255)                                                                        ttttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCg                                                    <  2:1600441/139‑1 (MQ=255)                                                                         tttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:1237372/136‑1 (MQ=255)                                                                         tttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      >  2:1237372/1‑136 (MQ=255)                                                                         tttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACACtgctgc                                                                                                             <  1:591372/81‑1 (MQ=255)                                                                         tttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACACtgctgc                                                                                                             >  2:591372/1‑81 (MQ=255)                                                                          ttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCg                                                                <  2:1458012/125‑1 (MQ=255)                                                                          ttGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCg                                                                >  1:1458012/1‑125 (MQ=255)                                                                                 tCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGtat                                                                                           >  1:734882/1‑91 (MQ=255)                                                                                 tCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGtat                                                                                           <  2:734882/91‑1 (MQ=255)                                                                                  cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                                          >  1:1563453/1‑139 (MQ=255)                                                                                  cgcgGCCTACCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATACCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                                          >  1:1144521/1‑139 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  2:665697/119‑1 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      >  1:665697/1‑119 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                        <  1:479290/133‑1 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                                        >  2:479290/1‑133 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCg                                  <  2:1144521/139‑1 (MQ=255)                                                                                          aCCATCAGACGGACGATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCg                                  <  2:1563453/139‑1 (MQ=255)                                                                                                 gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGTCGTACCctc                           <  2:286183/139‑1 (MQ=255)                                                                                                 gacggacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCctc                           >  2:607447/1‑139 (MQ=255)                                                                                                     gacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                             <  1:1451049/53‑1 (MQ=255)                                                                                                     gacgATATCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCtgctgc                                                                                                             >  2:1451049/1‑53 (MQ=255)                                                                                                            tCGTCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACt                <  1:607447/139‑1 (MQ=255)                                                                                                               tCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:863411/98‑1 (MQ=255)                                                                                                               tCATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      >  2:863411/1‑98 (MQ=255)                                                                                                                cATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      >  2:383349/1‑97 (MQ=255)                                                                                                                cATTTGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:383349/97‑1 (MQ=255)                                                                                                                    tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      >  2:1258028/1‑93 (MQ=255)                                                                                                                    tGACACATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                                      <  1:1258028/93‑1 (MQ=255)                                                                                                                      acacATCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGAt      <  1:993961/139‑1 (MQ=255)                                                                                                                          aTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACg                     <  2:1250958/120‑1 (MQ=255)                                                                                                                          aTCACTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGa  <  1:1123333/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                        TGACGTTCCGCACCTCAAAAACCGTTTCGTGCGCAAGGTTGCTGACATTCTCCTTCATCGCGTTAAAACGGTTTTGCAGTTCGCGGCCTACCATCAGACGGACGATATCGTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGA  >  NC_000913/4310167‑4310426

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