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%	202.7 / NA	62	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 (30/32);  total (30/32)

TATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCAT  >  NC_000913/4427949‑4428198                                                                                                                                   |                                                                                                                        tattatTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCgcgc                                                                                                                 <  1:498350/139‑1 (MQ=255)      tttGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGaa                                                                                                            <  1:195256/139‑1 (MQ=255)          aTCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATaaa                                                                                                        >  1:436857/1‑139 (MQ=255)                     tACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCg                                                                                             >  1:817737/1‑139 (MQ=255)                         cTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                         >  1:60426/1‑139 (MQ=255)                            ggCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCg                                                                                      <  1:330749/139‑1 (MQ=255)                                     tATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAg                                                                             <  1:1570642/139‑1 (MQ=255)                                        cacaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                               <  2:1319158/118‑1 (MQ=255)                                        cacaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                               >  1:1319158/1‑118 (MQ=255)                                        cacaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          <  2:817737/139‑1 (MQ=255)                                        cacaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          <  1:798154/139‑1 (MQ=255)                                             cAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGtt                                                                     >  1:22663/1‑139 (MQ=255)                                                     tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTc                                                             <  1:1036413/139‑1 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                          >  1:878494/1‑108 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                          <  2:878494/108‑1 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          >  1:220751/1‑124 (MQ=255)                                                       tCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          <  2:220751/124‑1 (MQ=255)                                                        cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTcccc                                                          <  2:4708/139‑1 (MQ=255)                                                         aCCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCg                                                         >  1:812578/1‑139 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                              >  1:689576/1‑133 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGattatt                                                              <  2:689576/133‑1 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                        <  2:675529/139‑1 (MQ=255)                                                          ccccAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGt                                                        <  2:22663/139‑1 (MQ=255)                                                              aGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                    <  2:436857/139‑1 (MQ=255)                                                                  aCTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATgg                                                >  1:1726796/1‑139 (MQ=255)                                                                     ggAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTGTTCCCCGTATCAATGGTTa                                             >  2:798419/1‑139 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                               <  2:764887/70‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATc                                                                                               >  1:764887/1‑70 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                          >  1:65452/1‑75 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCt                                                                                          <  2:65452/75‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGtt                                                                     <  2:896536/96‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGtt                                                                     >  1:896536/1‑96 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGt                                                                      <  2:1701491/95‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGt                                                                      >  1:1701491/1‑95 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatt                                                                 >  1:1013850/1‑100 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCa                                                    >  2:1451929/1‑113 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTGCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTAGCa                                                    <  1:1451929/113‑1 (MQ=255)                                                                                        atcaAAGAGACAATCCTTAGCTGGCTTTTGGTAGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGatt                                                                 <  2:1013850/100‑1 (MQ=255)                                                                                         tcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          <  1:441201/90‑1 (MQ=255)                                                                                         tcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          >  2:441201/1‑90 (MQ=255)                                                                                          caAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTgg                        <  1:620498/139‑1 (MQ=255)                                                                                                  cAATCCTTAGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg                >  1:843236/1‑139 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          >  2:1302416/1‑73 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGAc                                                                          <  1:1302416/73‑1 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt        <  2:1726796/139‑1 (MQ=255)                                                                                                          aGCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt        >  2:357115/1‑139 (MQ=255)                                                                                                           gCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGa                                                                                 <  2:839376/65‑1 (MQ=255)                                                                                                           gCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGa                                                                                 >  1:839376/1‑65 (MQ=255)                                                                                                           gCTGGCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTGTTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTc       <  1:798419/139‑1 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg                <  1:1034671/126‑1 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg                >  2:1034671/1‑126 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCa   >  1:1265161/1‑139 (MQ=255)                                                                                                               gCTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCa   >  2:718641/1‑139 (MQ=255)                                                                                                                cTTTTGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCAt  >  2:1325129/1‑139 (MQ=255)                                                                                                                    tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGc               <  2:146315/122‑1 (MQ=255)                                                                                                                    tGGTTGACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGc               >  1:146315/1‑122 (MQ=255)                                                                                                                         gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                         >  1:757102/1‑43 (MQ=255)                                                                                                                         gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTg                                                                                         <  2:757102/43‑1 (MQ=255)                                                                                                                         gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGc                                                                                       <  1:303561/45‑1 (MQ=255)                                                                                                                         gACCCTTCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGc                                                                                       >  2:303561/1‑45 (MQ=255)                                                                                                                              ttCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCggg                                                                                    >  2:1012233/1‑43 (MQ=255)                                                                                                                              ttCATTTACCGCGCACGAATAAAACTGCCATCCGGCTGGCggg                                                                                    <  1:1012233/43‑1 (MQ=255)                                                                                                                                   |                                                                                                                        TATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCAT  >  NC_000913/4427949‑4428198

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