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%	185.8 / NA	57	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 (28/29);  total (28/29)

ATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAA  >  NC_000913/4427948‑4428213                                                                                                                                    |                                                                                                                                       atatTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACcgcg                                                                                                                                 <  1:901651/139‑1 (MQ=255)           aTCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATaaa                                                                                                                       >  1:563935/1‑139 (MQ=255)                         gCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                                         <  2:804957/139‑1 (MQ=255)                                   gTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGaa                                                                                               >  1:1179683/1‑139 (MQ=255)                                       aTCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAgg                                                                                           <  2:1188480/139‑1 (MQ=255)                                             tCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                            <  2:623757/132‑1 (MQ=255)                                             tCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                            >  1:623757/1‑132 (MQ=255)                                             tCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGt                                                                                     >  1:1374506/1‑139 (MQ=255)                                              cAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGCCAATCCTTAGCTGGCTTTTGGCTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGtt                                                                                    >  1:476506/1‑139 (MQ=255)                                              cAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATa                                                                                                                         >  2:538152/1‑102 (MQ=255)                                              cAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATa                                                                                                                         <  1:538152/102‑1 (MQ=255)                                                 tATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                         <  2:408019/131‑1 (MQ=255)                                                 tATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                         >  1:81703/1‑131 (MQ=255)                                                 tATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                         >  1:408019/1‑131 (MQ=255)                                                 tATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                                         <  2:81703/131‑1 (MQ=255)                                                        tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGCCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTccc                                                                          <  2:476506/139‑1 (MQ=255)                                                        tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTccc                                                                          >  2:230263/1‑139 (MQ=255)                                                         cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                                              >  2:1085159/1‑102 (MQ=255)                                                         cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                                              <  1:1085159/102‑1 (MQ=255)                                                                ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGa                                                                                          <  1:1167280/115‑1 (MQ=255)                                                                ggCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGa                                                                                          >  2:1167280/1‑115 (MQ=255)                                                                     tGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGtt                                                             <  2:358487/139‑1 (MQ=255)                                                                     tGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGtt                                                             >  2:611273/1‑139 (MQ=255)                                                                             aGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAAc                                                                                                                     >  2:531697/1‑75 (MQ=255)                                                                             aGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAAc                                                                                                                     <  1:531697/75‑1 (MQ=255)                                                                                      gTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcccc                                                                         >  2:474319/1‑110 (MQ=255)                                                                                      gTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcccc                                                                         <  1:474319/110‑1 (MQ=255)                                                                                         gtataaGAGAC‑AGCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                                            <  2:464654/83‑1 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGaa                                                                                                                           <  2:902117/57‑1 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGaa                                                                                                                           >  1:902117/1‑57 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                                             <  2:498630/103‑1 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                                             >  1:498630/1‑103 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCt                                         >  2:729393/1‑139 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTG‑TTATTCCCCGTATCa                                                                   >  1:434926/1‑112 (MQ=255)                                                                                         atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTG‑TTATTCCCCGTATCa                                                                   <  2:434926/112‑1 (MQ=255)                                                                                                agaCAATCCTTAGCTGGCCTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTc                                  >  2:418140/1‑139 (MQ=255)                                                                                                   cAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg                               <  1:611273/139‑1 (MQ=255)                                                                                                     aTCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCg                             <  2:1374506/139‑1 (MQ=255)                                                                                                      tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACg                                                                                                                             <  1:726272/42‑1 (MQ=255)                                                                                                      tCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACg                                                                                                                             >  2:726272/1‑42 (MQ=255)                                                                                                          tAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAg                                                          >  2:841171/1‑105 (MQ=255)                                                                                                          tAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAg                                                          <  1:841171/105‑1 (MQ=255)                                                                                                             cTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCa                     <  1:729393/139‑1 (MQ=255)                                                                                                               ggCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCAcc                   <  2:1179683/139‑1 (MQ=255)                                                                                                                gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                    >  2:169042/1‑57 (MQ=255)                                                                                                                gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCgg                                                                                                    <  1:169042/57‑1 (MQ=255)                                                                                                                     tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCt                                                                                       >  1:964841/1‑65 (MQ=255)                                                                                                                     tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCt                                                                                       <  2:964841/65‑1 (MQ=255)                                                                                                                     tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTAt                                                                     <  2:29843/83‑1 (MQ=255)                                                                                                                     tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTAt                                                                     >  1:29843/1‑83 (MQ=255)                                                                                                                      ggTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGcc            >  1:897164/1‑139 (MQ=255)                                                                                                                          gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                                         <  2:1408239/42‑1 (MQ=255)                                                                                                                          gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                                         >  1:1408239/1‑42 (MQ=255)                                                                                                                          gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGtt        <  1:418140/139‑1 (MQ=255)                                                                                                                           aCCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGt                                                                                     <  2:838914/61‑1 (MQ=255)                                                                                                                           aCCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGt                                                                                     >  1:838914/1‑61 (MQ=255)                                                                                                                                tCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCaa  >  2:569673/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       ATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAA  >  NC_000913/4427948‑4428213

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