Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	826,119	C→T	*378* (TAG→TAA)	ybhS ←	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	826,119	0	C	T	100.0%	85.6 / NA	28	*378* (TAG→TAA)	ybhS	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/15);  total (13/15)

GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCA  >  NC_000913/825987‑826243                                                                                                                                     |                                                                                                                             gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                        >  2:373987/1‑139 (MQ=255)                      aCCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                   >  2:280299/1‑139 (MQ=255)                       ccTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGcc                                                                                                  >  1:369576/1‑139 (MQ=255)                          gAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACg                                                                                               >  2:393108/1‑139 (MQ=255)                                gCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGATTTTCAGCGACGTCAGGc                                                                                         >  2:394423/1‑139 (MQ=255)                                  aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCCCACGAACGCCTGCAACTCTTTGCGGATTAACGTCCCTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCg                                                                                       <  1:777444/139‑1 (MQ=255)                                     ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATa                                                                                    >  1:386207/1‑139 (MQ=255)                                         aaaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                   >  1:193497/1‑120 (MQ=255)                                         aaaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGc                                                                                                   <  2:193497/120‑1 (MQ=255)                                          aaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                               <  2:369576/139‑1 (MQ=255)                                          aaTCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGACCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                               <  1:240748/139‑1 (MQ=255)                                                aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCAccgc                                                                         >  1:347297/1‑139 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                               <  1:639377/128‑1 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                               >  2:639377/1‑128 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                    >  1:456079/1‑139 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                    <  1:280299/139‑1 (MQ=255)                                                     cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTATCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                    <  2:464898/139‑1 (MQ=255)                                                        gggTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGa                                                                 <  2:708029/139‑1 (MQ=255)                                                                ggTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAAc                                                         <  2:379411/139‑1 (MQ=255)                                                                   tCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGc                                                      <  2:372160/139‑1 (MQ=255)                                                                          gcaACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTa                                               <  1:612799/139‑1 (MQ=255)                                                                            aaCGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTacc                                             <  2:456079/139‑1 (MQ=255)                                                                               gACTGCTACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTaccacc                                          >  2:773599/1‑139 (MQ=255)                                                                                 cTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAg                                        <  1:300490/139‑1 (MQ=255)                                                                                   gCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                      <  1:598202/139‑1 (MQ=255)                                                                                                    aaCGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCg                     >  2:202968/1‑139 (MQ=255)                                                                                                           aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGa              >  2:117775/1‑139 (MQ=255)                                                                                                                       acttGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGACAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCa  <  1:373987/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                             GCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCA  >  NC_000913/825987‑826243

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