Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	675,570	C→T	*326* (TAG→TAA)	ybeQ ←	Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	675,570	0	C	T	100.0%	124.3 / NA	39	*326* (TAG→TAA)	ybeQ	Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/19);  total (20/19)

GAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGG  >  NC_000913/675435‑675705                                                                                                                                        |                                                                                                                                        gAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTAc                                                                                                                                      >  2:668425/1‑139 (MQ=255)      gCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTaa                                                                                                                                 <  1:543872/139‑1 (MQ=255)               aTGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGt                                                                                                                        >  1:34710/1‑139 (MQ=255)                           tgGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            <  2:81068/139‑1 (MQ=255)                           tgGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            <  2:424064/139‑1 (MQ=255)                                          ccAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTg                                                                                            >  1:195819/1‑139 (MQ=255)                                                  gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                    >  2:70213/1‑139 (MQ=255)                                                  gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                    >  2:432645/1‑139 (MQ=255)                                                                        tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAgg                                                              >  2:1034695/1‑139 (MQ=255)                                                                          aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACa                                                                <  1:98695/135‑1 (MQ=255)                                                                          aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACa                                                                >  2:98695/1‑135 (MQ=255)                                                                          aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTc                                                            >  1:846531/1‑139 (MQ=255)                                                                                         aaaGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGgag                                             <  1:842906/139‑1 (MQ=255)                                                                                             agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTccgact                                                                                                       >  1:745495/1‑74 (MQ=255)                                                                                             agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTccgact                                                                                                       <  2:745495/78‑5 (MQ=255)                                                                                               tatatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  1:552820/2‑139 (MQ=255)                                                                                               aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGGGATCGAc                                       >  2:560942/1‑139 (MQ=255)                                                                                               aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  1:699544/1‑139 (MQ=255)                                                                                               aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  1:100577/1‑139 (MQ=255)                                                                                               aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  2:1099326/1‑139 (MQ=255)                                                                                               aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  2:1106132/1‑139 (MQ=255)                                                                                                 tatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            <  1:356496/100‑1 (MQ=255)                                                                                                 tatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            >  2:356496/1‑100 (MQ=255)                                                                                                       ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTct                                 <  1:62917/137‑1 (MQ=255)                                                                                                       ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTct                                 >  2:62917/1‑137 (MQ=255)                                                                                                         gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGa                                        >  2:833567/1‑128 (MQ=255)                                                                                                         gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGa                                        <  1:833567/128‑1 (MQ=255)                                                                                                         gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcct                             >  1:157278/1‑139 (MQ=255)                                                                                                         gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                               <  2:872896/137‑1 (MQ=255)                                                                                                         gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                               >  1:872896/1‑137 (MQ=255)                                                                                                             gAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGGCAGGTCTTCCTtctc                         <  2:34710/139‑1 (MQ=255)                                                                                                               aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                       <  1:432645/139‑1 (MQ=255)                                                                                                                 ttCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTa                     <  1:4795/139‑1 (MQ=255)                                                                                                                  tCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAg                    <  2:157278/139‑1 (MQ=255)                                                                                                                    cccGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTg                  <  1:1099326/139‑1 (MQ=255)                                                                                                                                ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc      <  2:699544/139‑1 (MQ=255)                                                                                                                                ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc      <  2:100577/139‑1 (MQ=255)                                                                                                                                ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc      <  1:1034695/139‑1 (MQ=255)                                                                                                                                    aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg  <  1:1106132/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGG  >  NC_000913/675435‑675705

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