Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,535,441	C→T	pseudogene (439/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,441	0	C	T	100.0%	123.1 / NA	39	pseudogene (439/1029 nt)	yjhR	pseudogene, helicase family, putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/19);  total (20/19)

TGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTG  >  NC_000913/4535317‑4535570                                                                                                                             |                                                                                                                                  tGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgccg                                                                                                                     <  1:335336/139‑1 (MQ=255)     tACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAg                                                                                                                 >  1:212169/1‑139 (MQ=255)       cTCGAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGGATTTTCGCCGCGAggg                                                                                                               >  2:546642/1‑139 (MQ=255)       cTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAggg                                                                                                               >  1:201117/1‑139 (MQ=255)            aaCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAtt                                                                                                          >  2:245801/1‑139 (MQ=255)                 gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaa                                                                                                     >  1:117006/1‑139 (MQ=255)                      tGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                                <  1:612135/139‑1 (MQ=255)                       ggCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATAtc                                                                                               >  1:69195/1‑139 (MQ=255)                        gCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATAtct                                                                                              >  2:90721/1‑139 (MQ=255)                             tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                         <  1:546642/139‑1 (MQ=255)                               tCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                       <  2:201117/139‑1 (MQ=255)                                  tGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCa                                                                                    <  2:85881/139‑1 (MQ=255)                                    aTAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGa                                                                                  <  1:245801/139‑1 (MQ=255)                                           cAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATaataat                                                                           >  2:231505/1‑139 (MQ=255)                                               aCTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGt                                                                       >  1:266724/1‑139 (MQ=255)                                                             tGCTGTTTCACGCGCTATAGATAGTGTTCTGGGATTTGGCGAGATGGAGCTTATCGAGGGGGAGGCCGCATTTTCGCCGTGAGGGTGATTGGCAAAATAGCTGTTTTCTTCAGATAATAATGCGTAGAAATTTGAGttt                                                         >  1:29817/1‑139 (MQ=255)                                                                      aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                     <  2:59419/102‑1 (MQ=255)                                                                      aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                     >  1:59419/1‑102 (MQ=255)                                                                      aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                <  2:266724/139‑1 (MQ=255)                                                                       cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                               <  2:117006/139‑1 (MQ=255)                                                                            tAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGAt                                          >  1:217677/1‑139 (MQ=255)                                                                                  tAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAAtt                                    >  2:329206/1‑139 (MQ=255)                                                                                      tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                          >  1:265197/1‑49 (MQ=255)                                                                                      tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCAtttt                                                                                                                          <  2:265197/49‑1 (MQ=255)                                                                                      tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                <  2:109061/139‑1 (MQ=255)                                                                                      tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                <  2:212169/139‑1 (MQ=255)                                                                                      tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                <  2:224448/139‑1 (MQ=255)                                                                                        tCTGGTATTTGGCGATATGGGCCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAACTATCTCTTttctt                                                                                      >  1:471503/1‑83 (MQ=255)                                                                                        tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAATTATCTCTTttctt                                                                                      <  2:471503/83‑1 (MQ=255)                                                                                        tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCa                                                                                    <  1:153229/85‑1 (MQ=255)                                                                                        tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCa                                                                                    >  2:153229/1‑85 (MQ=255)                                                                                              aTTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATaca                        <  1:231505/139‑1 (MQ=255)                                                                                               tttGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAGCGACAGGATTTAATTTCCGCACATacac                       >  1:582920/1‑139 (MQ=255)                                                                                                 tGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACaa                     >  1:21397/1‑139 (MQ=255)                                                                                                      atatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTc                >  1:325006/1‑139 (MQ=255)                                                                                                               ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa       <  1:90721/139‑1 (MQ=255)                                                                                                               ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa       <  2:69195/139‑1 (MQ=255)                                                                                                                 aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATg     >  1:406613/1‑139 (MQ=255)                                                                                                                 aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATg     >  1:428841/1‑139 (MQ=255)                                                                                                                    gagaTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGgtg  <  1:329206/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                  TGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTG  >  NC_000913/4535317‑4535570

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