Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,544,384	C→T	pseudogene (381/381 nt)	yddK ←	pseudogene, leucine‑rich protein, putative glycoportein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,544,384	0	C	T	100.0%	133.2 / NA	42	pseudogene (381/381 nt)	yddK	pseudogene, leucine‑rich protein, putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/22);  total (20/22)

AAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTT  >  NC_000913/1544260‑1544514                                                                                                                             |                                                                                                                                   aaaCGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaaa                                                                                                                               >  2:450365/1‑130 (MQ=255) aaaCGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaaa                                                                                                                               <  1:450365/130‑1 (MQ=255)       ggTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGaa                                                                                                                     <  2:417963/134‑1 (MQ=255)       ggTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGaa                                                                                                                     >  1:417963/1‑134 (MQ=255)         ttgaTCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTggagg                                                                                                              >  2:487036/4‑139 (MQ=255)                  gAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATc                                                                                                     <  1:752212/139‑1 (MQ=255)                                tttCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATAAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCt                                                                                       >  2:601242/1‑139 (MQ=255)                                         ttCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTcc                                                                              >  1:275156/1‑139 (MQ=255)                                                                  aatatTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                     <  1:250875/138‑1 (MQ=255)                                                                    tatTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCg                                                   <  1:252333/139‑1 (MQ=255)                                                                         aCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                              >  1:91219/1‑139 (MQ=255)                                                                            atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACa                                           <  1:291722/139‑1 (MQ=255)                                                                             tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAtt                                                             <  2:249827/120‑1 (MQ=255)                                                                             tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAtt                                                             >  1:249827/1‑120 (MQ=255)                                                                              attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                         <  1:125717/139‑1 (MQ=255)                                                                              attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                         <  2:33294/139‑1 (MQ=255)                                                                              attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                         <  1:601242/139‑1 (MQ=255)                                                                               ttatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagta                                        >  1:364137/1‑139 (MQ=255)                                                                                tatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaa                                       <  2:364137/139‑1 (MQ=255)                                                                                     acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                               >  1:816758/1‑126 (MQ=255)                                                                                     acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                               <  2:816758/126‑1 (MQ=255)                                                                                       aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGt                                >  2:423620/1‑139 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAGGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             >  2:194397/1‑92 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAGGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             <  1:194397/92‑1 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAg                                                                                              <  1:518473/75‑1 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAg                                                                                              >  2:518473/1‑75 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                 <  1:744958/88‑1 (MQ=255)                                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                 >  2:744958/1‑88 (MQ=255)                                                                                            cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             >  1:448343/1‑89 (MQ=255)                                                                                            cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             <  2:448343/89‑1 (MQ=255)                                                                                              gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                       >  2:260958/1‑109 (MQ=255)                                                                                              gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                       <  1:260958/109‑1 (MQ=255)                                                                                              gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCg                                                   <  2:273663/113‑1 (MQ=255)                                                                                              gACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCg                                                   >  1:273663/1‑113 (MQ=255)                                                                                                     ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGATCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAg                  >  2:538838/1‑139 (MQ=255)                                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             <  1:162703/72‑1 (MQ=255)                                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                             >  2:162703/1‑72 (MQ=255)                                                                                                               tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                     >  2:424187/1‑94 (MQ=255)                                                                                                               tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                                     <  1:424187/94‑1 (MQ=255)                                                                                                               tgatgaTTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtttt        >  2:292515/1‑139 (MQ=255)                                                                                                                    atTTA‑TCATTAAAAAGGCCCTTAACTGGAGGGAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACtt  <  1:186342/139‑1 (MQ=255)                                                                                                                     tttATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACtt  <  2:91219/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                   AAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTT  >  NC_000913/1544260‑1544514

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