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%	161.5 / NA	50	*378* (TAG→TAA)	ybhS	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/23);  total (27/23)

TGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGC  >  NC_000913/825985‑826249                                                                                                                                       |                                                                                                                                   tGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAAt                                                                                                                                <  1:852564/139‑1 (MQ=255)   gcggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATcc                                                                                                                              >  2:996288/1‑139 (MQ=255)     ggcgAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAg                                                                                                                            <  1:363910/139‑1 (MQ=255)                        aCCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCg                                                                                                                         <  2:1062748/123‑1 (MQ=255)                        aCCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCg                                                                                                                         >  1:1062748/1‑123 (MQ=255)                         ccTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGcc                                                                                                        <  1:633387/139‑1 (MQ=255)                             aaTTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGt                                                                                                    <  1:996288/139‑1 (MQ=255)                               ttAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCa                                                                                                  <  2:1058553/139‑1 (MQ=255)                                    aCGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCg                                                                                             >  2:43268/1‑139 (MQ=255)                                       ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                                                >  1:266858/1‑133 (MQ=255)                                       ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAgg                                                                                                <  2:266858/133‑1 (MQ=255)                                       ggTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATa                                                                                          >  1:994679/1‑139 (MQ=255)                                        gTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaa                                                                                         >  1:186694/1‑139 (MQ=255)                                        gTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATaa                                                                                         >  2:474336/1‑139 (MQ=255)                                                aGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                     >  1:993611/1‑135 (MQ=255)                                                aGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                     <  2:993611/135‑1 (MQ=255)                                                  aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGc                                                                                               >  2:626774/1‑123 (MQ=255)                                                  aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGc                                                                                               <  1:626774/123‑1 (MQ=255)                                                  aaTCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCAccgc                                                                               >  1:1163780/1‑139 (MQ=255)                                                       cgcgGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                          <  2:1163780/139‑1 (MQ=255)                                                                      cgcgCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                                           >  2:950519/1‑139 (MQ=255)                                                                      cgcgCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCa                                                           <  1:43268/139‑1 (MQ=255)                                                                                cGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTaccac                                                 <  2:128303/139‑1 (MQ=255)                                                                                      cAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAg                                                                                                 <  2:1203606/85‑1 (MQ=255)                                                                                      cAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAg                                                                                                 >  1:1203606/1‑85 (MQ=255)                                                                                         ctctTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcaccac                                        <  2:186694/139‑1 (MQ=255)                                                                                           ctTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcac                                           <  2:638141/134‑1 (MQ=255)                                                                                           ctTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGcac                                           >  1:638141/1‑134 (MQ=255)                                                                                              tGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                            >  1:226110/1‑130 (MQ=255)                                                                                              tGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGca                                            <  2:226110/130‑1 (MQ=255)                                                                                                cGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAATAACAGCACGTTTa                                                     <  1:819414/119‑1 (MQ=255)                                                                                                cGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTa                                                     >  2:819414/1‑119 (MQ=255)                                                                                                   aTTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTc                              >  1:96052/1‑139 (MQ=255)                                                                                                   aTTAACGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTc                              >  1:999248/1‑139 (MQ=255)                                                                                                      aaCGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCg                           >  2:560608/1‑139 (MQ=255)                                                                                                        cGTCCATAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGc                         <  2:96052/139‑1 (MQ=255)                                                                                                             aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                     <  2:559464/74‑1 (MQ=255)                                                                                                             aTAAGCGATGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACAt                                                                                     >  1:559464/1‑74 (MQ=255)                                                                                                                    aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGtt                                                                                                                    >  2:824276/1‑36 (MQ=255)                                                                                                                    aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                          <  2:861301/78‑1 (MQ=255)                                                                                                                    aTGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAg                                                                          >  1:861301/1‑78 (MQ=255)                                                                                                                    aTGAAACATGCACTTCTCCTTAATCCAGCCGACGtt                                                                                                                    <  1:824276/36‑1 (MQ=255)                                                                                                                     tGAAACATGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCtt            >  1:945419/1‑139 (MQ=255)                                                                                                                            tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAggg     >  2:272514/1‑139 (MQ=255)                                                                                                                            tGCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAggg     >  2:887289/1‑139 (MQ=255)                                                                                                                             gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAg                >  2:982834/1‑127 (MQ=255)                                                                                                                             gCTCTTCTCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAg                <  1:982834/127‑1 (MQ=255)                                                                                                                               tcttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGc  >  2:28546/1‑139 (MQ=255)                                                                                                                                 ttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTcc                             <  2:472445/110‑1 (MQ=255)                                                                                                                                 ttctCCTTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTcc                             >  1:472445/1‑110 (MQ=255)                                                                                                                                       |                                                                                                                                   TGGCGGCGAACGGGAACAGGATCACCTGAATTAGCACGGGTAAAATCAGAATCGCGCGGGTTTGCGGTTCGCGCAGCAACGACTGCAACTCTTTGCGGATTAACGTCCATAAGCGATGAAACATGCTCTTCTCCCTAATCCAGCCGACGTTTGGTTTTCAGCCACGTCAGGCCGATAAACATCACCGCCGAAGCGATCAAAAACAGCACGTTTACCACCAGCACCACTGGAATATTCCCGGCGAGGAACAGGCTTTGCAGGGTGC  >  NC_000913/825985‑826249

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