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%	126.2 / NA	40	*326* (TAG→TAA)	ybeQ	Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/21);  total (19/21)

CGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGG  >  NC_000913/675445‑675705                                                                                                                              |                                                                                                                                        cGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCtt                                                                                                                            <  1:320767/139‑1 (MQ=255)                  ggTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTcc                                                                                                           >  2:763510/1‑139 (MQ=255)                          ccAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACa                                                                                                  <  2:703833/139‑1 (MQ=255)                                        gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                    >  1:229973/1‑139 (MQ=255)                                             tttGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACt                                                                               >  2:833112/1‑139 (MQ=255)                                                      tAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTAc                                                                      >  1:744293/1‑139 (MQ=255)                                                        aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                      <  2:45945/121‑1 (MQ=255)                                                        aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                      >  1:364907/1‑121 (MQ=255)                                                        aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                      <  2:364907/121‑1 (MQ=255)                                                        aTAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATATTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCg                                                                                      >  1:45945/1‑121 (MQ=255)                                                          aGTCTGATTTCGAACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATa                                                                  <  1:110013/139‑1 (MQ=255)                                                              tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTcc‑gactga                                                                                                    >  2:502557/1‑101 (MQ=255)                                                              tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTcc‑gactga                                                                                                    <  1:502557/101‑1 (MQ=255)                                                              tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            <  1:690382/125‑1 (MQ=255)                                                              tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            >  2:690382/1‑125 (MQ=255)                                                              tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAgg                                                              >  2:408233/1‑139 (MQ=255)                                                                    cGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            >  1:479549/1‑119 (MQ=255)                                                                    cGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTCGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCGGGTTCCGGATGACTGaa                                                                            <  2:479549/119‑1 (MQ=255)                                                                      cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            >  1:485580/1‑86 (MQ=255)                                                                      cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            <  2:485580/86‑1 (MQ=255)                                                                                   tgaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                    <  2:545393/95‑1 (MQ=255)                                                                                   agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGa                                                                                    >  1:545393/1‑96 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCtt                                                                                                                            <  2:20742/55‑1 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCtt                                                                                                                            >  1:20742/1‑55 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            <  2:540282/71‑1 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTc                                                                                                            >  1:540282/1‑71 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       <  2:229973/139‑1 (MQ=255)                                                                                     aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                       >  2:426179/1‑139 (MQ=255)                                                                                       tatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAg                                     <  2:744293/139‑1 (MQ=255)                                                                                             ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                               >  2:157245/1‑139 (MQ=255)                                                                                              cGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcc                              >  2:455706/1‑139 (MQ=255)                                                                                               gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            >  1:610421/1‑92 (MQ=255)                                                                                               gAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGaa                                                                            <  2:610421/92‑1 (MQ=255)                                                                                                aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTa                                                                                                                                  <  2:427889/38‑1 (MQ=255)                                                                                                aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTa                                                                                                                                  >  1:427889/1‑38 (MQ=255)                                                                                                       ttCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTa                     <  1:220254/139‑1 (MQ=255)                                                                                                                      ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc      <  1:763510/139‑1 (MQ=255)                                                                                                                      ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc      <  1:408233/139‑1 (MQ=255)                                                                                                                          aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg  <  1:426179/139‑1 (MQ=255)                                                                                                                          aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCC‑GACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg  <  1:157245/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                        CGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGG  >  NC_000913/675445‑675705

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