Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,254	A→G	pseudogene (252/1029 nt)	yjhR →	pseudogene, helicase family; putative frameshift suppressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,535,254	0	A	G	100.0%	196.2 / NA	55	pseudogene (252/1029 nt)	yjhR	pseudogene, helicase family; putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/30);  total (25/30)

TGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAG  >  NC_000913/4535106‑4535396                                                                                                                                                     |                                                                                                                                               tggtggCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCt                                                                                                                                              >  2:136161/1‑151 (MQ=255) tggtggCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCt                                                                                                                                              >  1:42635/1‑151 (MQ=255)        acaTAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTg                                                                                                                                       <  1:313705/151‑1 (MQ=255)            aaGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTgggg                                                                                                                                    >  2:321938/1‑150 (MQ=255)            aaGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTgggg                                                                                                                                    >  2:657490/1‑150 (MQ=255)                  gatATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACtgtg                                                                                                                              <  1:562440/150‑1 (MQ=255)                   atatAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCa                                                                                                                            <  2:681944/151‑1 (MQ=255)                   atatAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCa                                                                                                                            <  2:339199/151‑1 (MQ=255)                         agCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCtt                                                                                                                      >  2:296555/1‑151 (MQ=255)                            gcCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCa                                                                                                                   <  1:804275/151‑1 (MQ=255)                             ccATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTg                                                                                                                     >  2:909516/1‑148 (MQ=255)                              cATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAgg                                                                                                                 >  2:715278/1‑151 (MQ=255)                              cATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAgg                                                                                                                 <  1:330219/151‑1 (MQ=255)                                        aaCCACTGTACAAATTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGaa                                                                                                         >  1:36612/1‑149 (MQ=255)                                                    aaGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTtct                                                                                           <  1:160118/151‑1 (MQ=255)                                                     aGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTtctc                                                                                          >  1:360798/1‑151 (MQ=255)                                                         gttgGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCtttt                                                                                      <  1:489354/151‑1 (MQ=255)                                                         gttgGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCttt                                                                                       <  1:122251/150‑1 (MQ=255)                                                            ggTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTtctc                                                                                   >  1:793460/1‑151 (MQ=255)                                                            ggTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTtct                                                                                    <  2:226069/150‑1 (MQ=255)                                                               gttgtGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTg                                                                                <  2:383259/151‑1 (MQ=255)                                                                 tgtgACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGt                                                                               <  1:196207/150‑1 (MQ=255)                                                                    gACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTAc                                                                           >  2:416753/1‑151 (MQ=255)                                                                       gCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCa                                                                        <  1:442319/151‑1 (MQ=255)                                                                       gCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGAGGGGACTGTGCAATCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTc                                                                         >  2:262526/1‑150 (MQ=255)                                                                         cATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCaa                                                                       >  1:510407/1‑150 (MQ=255)                                                                         cATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCaa                                                                       >  1:780711/1‑150 (MQ=255)                                                                          aTTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCaaaa                                                                     <  2:65573/151‑1 (MQ=255)                                                                          aTTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCaaaa                                                                     <  1:259263/151‑1 (MQ=255)                                                                              tCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACAt                                                                  <  2:250763/150‑1 (MQ=255)                                                                                ggCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGaa                                                               <  2:715218/151‑1 (MQ=255)                                                                                         tAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAg                                                       <  1:699688/150‑1 (MQ=255)                                                                                         tAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGAAGGCAg                                                       <  1:112990/150‑1 (MQ=255)                                                                                             gCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGttt                                                   <  1:744083/150‑1 (MQ=255)                                                                                             gCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGttt                                                   <  1:650405/150‑1 (MQ=255)                                                                                                  cAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTtct                                                                                           <  2:336051/105‑1 (MQ=255)                                                                                                  cAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTtct                                                                                           >  1:336077/1‑105 (MQ=255)                                                                                                   aaaaTGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGAt                                             >  1:876331/1‑150 (MQ=255)                                                                                                        gTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAAc                                       >  2:336182/1‑151 (MQ=255)                                                                                                                  aaGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATcc                             <  1:757962/151‑1 (MQ=255)                                                                                                                  aaGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATcc                             <  1:651680/151‑1 (MQ=255)                                                                                                                   aGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATcc                             >  2:99829/1‑150 (MQ=255)                                                                                                                   aGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCt                            >  1:310215/1‑151 (MQ=255)                                                                                                                   aGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCt                            <  1:132243/151‑1 (MQ=255)                                                                                                                    gCTGGAAATTAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTc                           <  2:444879/151‑1 (MQ=255)                                                                                                                            ttAATGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGc                   >  1:678998/1‑151 (MQ=255)                                                                                                                              aaTGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTg                 >  2:871643/1‑151 (MQ=255)                                                                                                                               aTGGCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGt                >  1:540675/1‑151 (MQ=255)                                                                                                                                 ggCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGtct              >  1:378338/1‑151 (MQ=255)                                                                                                                                 ggCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGtc               <  1:2785/150‑1 (MQ=255)                                                                                                                                  gCAAAGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGTAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGtc               >  2:861425/1‑149 (MQ=255)                                                                                                                                    aaaGACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAc           <  1:420346/151‑1 (MQ=255)                                                                                                                                       gACGAGCAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCAcgcg        <  1:507533/151‑1 (MQ=255)                                                                                                                                             cAGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAg  <  1:253186/151‑1 (MQ=255)                                                                                                                                              aGGGTTTGCTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTaa    >  2:583313/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                               TGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAG  >  NC_000913/4535106‑4535396

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