Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,985,520	(C)5→6	coding (609/990 nt)	araF ←	arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,985,515	1	.	C	100.0%	100.3 / NA	29	A205G (GCA→GGA)	araF	arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/13);  total (16/13)

GCCGATGATATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATG‑CCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGCC  >  NC_000913/1985376‑1985646                                                                                                                                             |                                                                                                                                    gccGATGATATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATg                                                                                                                                      <  2:1081780/140‑1 (MQ=255)              gCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGt                                                                                                                         >  1:17503/1‑140 (MQ=255)                           ccTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTggt                                                                                                            <  2:623619/140‑1 (MQ=255)                           ccTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTggt                                                                                                            <  2:252868/140‑1 (MQ=255)                             tGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTggtag                                                                                                          >  1:683086/1‑140 (MQ=255)                               aCCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTggtaggt                                                                                                        >  1:466927/1‑140 (MQ=255)                                 cTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTAc                                                                                                      <  2:4925/140‑1 (MQ=255)                                       tCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                                                                >  1:171111/1‑140 (MQ=255)                                          cgcgTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                                                             >  1:517820/1‑140 (MQ=255)                                           gcgTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAtt                                                                                            >  1:571115/1‑140 (MQ=255)                                                               tGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                                                                        >  2:1100494/1‑140 (MQ=255)                                                                        tCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                                                               <  2:517820/140‑1 (MQ=255)                                                                         cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                                                              >  2:251453/1‑140 (MQ=255)                                                                          aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                             <  2:171111/140‑1 (MQ=255)                                                                          aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                             <  2:592859/140‑1 (MQ=255)                                                                          aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                                                             <  1:240843/140‑1 (MQ=255)                                                                            aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCa                                                           >  1:848283/1‑140 (MQ=255)                                                                                    tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa                                                   <  2:731255/140‑1 (MQ=255)                                                                                            tGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgtag                                           >  1:442797/1‑140 (MQ=255)                                                                                             gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgtagt                                          >  1:5057/1‑140 (MQ=255)                                                                                                   cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCAt                                                      >  2:43238/1‑122 (MQ=255)                                                                                                   cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCAt                                                      <  1:43238/122‑1 (MQ=255)                                                                                                       cGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTAcggcggcgg                                >  1:1020019/1‑140 (MQ=255)                                                                                                        ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCggcggg                               >  1:964664/1‑140 (MQ=255)                                                                                                        ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCggcggg                               <  2:442797/140‑1 (MQ=255)                                                                                                        ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCggcggg                               <  2:571115/140‑1 (MQ=255)                                                                                                          aTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCggcgggcg                             <  2:964664/140‑1 (MQ=255)                                                                                                                                     gTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTa        >  2:868821/1‑134 (MQ=255)                                                                                                                                     gTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTa        <  1:868821/134‑1 (MQ=255)                                                                                                                                     gTCAAATGCCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGcc  >  2:362088/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                    GCCGATGATATCGGCCGCTTTAAAGCCCTGACCTTCCGTCGCGCGTACGCCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATG‑CCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCGGCGGGCGGTATCCAGTTCGTTGGCGGTAATCGCC  >  NC_000913/1985376‑1985646

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