Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,428,079	C→T	pseudogene (386/402 nt) *213* (TAG→TAA)	ytfA → ytfB ←	pseudogene, related to transcriptional regulators OapA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,428,079	0	C	T	100.0%	207.6 / NA	63	pseudogene (386/402 nt) *213* (TAG→TAA)	ytfA ytfB	pseudogene, related to transcriptional regulators OapA family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/30);  total (33/30)

TCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTC  >  NC_000913/4427959‑4428193                                                                                                                         |                                                                                                                   tCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACcgcg                                                                                                             >  1:665181/1‑128 (MQ=255) tCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACcgcg                                                                                                             <  2:665181/128‑1 (MQ=255) tCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATaaaa                                                                                                  >  2:817543/1‑139 (MQ=255)              gCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                     <  2:516139/139‑1 (MQ=255)                ttGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                          >  1:647386/1‑132 (MQ=255)                ttGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                          <  2:647386/132‑1 (MQ=255)                ttGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                   >  2:514348/1‑139 (MQ=255)                   gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCa                                                                                            <  1:1151296/127‑1 (MQ=255)                   gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCa                                                                                            >  2:1151296/1‑127 (MQ=255)                   gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                   >  1:1245337/1‑136 (MQ=255)                   gCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTgg                                                                                   <  2:1245337/136‑1 (MQ=255)                         tttATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAAc                                                                          <  2:1242089/139‑1 (MQ=255)                             tCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAAc                                                                          <  2:655862/135‑1 (MQ=255)                             tCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAAc                                                                          >  1:655862/1‑135 (MQ=255)                              cacaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                     <  1:817543/139‑1 (MQ=255)                                              cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                     <  1:1419622/123‑1 (MQ=255)                                              cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                     >  2:1419622/1‑123 (MQ=255)                                               aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                            <  1:431459/83‑1 (MQ=255)                                               aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                            >  2:431459/1‑83 (MQ=255)                                               aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                        <  2:1360824/119‑1 (MQ=255)                                               aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                        >  1:1360824/1‑119 (MQ=255)                                               aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCg                                                    >  1:360284/1‑139 (MQ=255)                                                ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                            >  1:1509196/1‑82 (MQ=255)                                                ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                            <  2:1509196/82‑1 (MQ=255)                                                ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                   >  1:287167/1‑139 (MQ=255)                                                ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                   >  2:25635/1‑139 (MQ=255)                                                    aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               <  2:287167/139‑1 (MQ=255)                                                    aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               <  2:360284/139‑1 (MQ=255)                                                    aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               <  1:514348/139‑1 (MQ=255)                                                     ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAgg                                                                       >  1:747749/1‑114 (MQ=255)                                                     ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAgg                                                                       <  2:747749/114‑1 (MQ=255)                                                     ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGa                                                              <  2:1100826/123‑1 (MQ=255)                                                     ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGa                                                              >  1:1100826/1‑123 (MQ=255)                                                     ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCaa                                              >  2:368855/1‑139 (MQ=255)                                                          tGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCgg                                                                                       >  1:527418/1‑93 (MQ=255)                                                          tGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCgg                                                                                       <  2:527418/93‑1 (MQ=255)                                                            gAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCggg                                                                               <  2:436993/99‑1 (MQ=255)                                                            gAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCggg                                                                               >  1:436993/1‑99 (MQ=255)                                                                 cAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGATTAAAACTGCCATCCGGCt                                                                                     <  2:838156/88‑1 (MQ=255)                                                                 cAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                     >  1:838156/1‑88 (MQ=255)                                                                 cAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTCGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAAGGGTTAAGCCCg                                  >  1:986735/1‑139 (MQ=255)                                                                              atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               <  1:1301822/113‑1 (MQ=255)                                                                              atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               >  2:1301822/1‑113 (MQ=255)                                                                                        cAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg           >  2:190439/1‑139 (MQ=255)                                                                                          aTCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCg         >  2:363983/1‑139 (MQ=255)                                                                                           tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               >  2:241103/1‑100 (MQ=255)                                                                                           tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                               <  1:241103/100‑1 (MQ=255)                                                                                           tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTa                                        >  1:1001238/1‑107 (MQ=255)                                                                                           tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTa                                        <  2:1001238/107‑1 (MQ=255)                                                                                             cTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAt      >  1:146582/1‑139 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                     <  1:427418/57‑1 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                     >  2:427418/1‑57 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                >  1:14919/1‑62 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                <  2:14919/62‑1 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGccc                                   >  1:1463156/1‑107 (MQ=255)                                                                                                aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGccc                                   <  2:1463156/107‑1 (MQ=255)                                                                                                 gCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTc  <  1:25635/139‑1 (MQ=255)                                                                                                      cTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                <  2:1099060/56‑1 (MQ=255)                                                                                                      cTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                >  2:1199928/1‑56 (MQ=255)                                                                                                      cTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                <  1:1199928/56‑1 (MQ=255)                                                                                                      cTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                >  1:1099060/1‑56 (MQ=255)                                                                                                      cTTTTGGCTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAg                                      <  2:1452955/98‑1 (MQ=255)                                                                                                      cTTTTGGCTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAg                                      >  1:1452955/1‑98 (MQ=255)                                                                                                                         |                                                                                                                   TCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTC  >  NC_000913/4427959‑4428193

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