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%	249.7 / NA	77	pseudogene (381/381 nt)	yddK	pseudogene, leucine‑rich protein, putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/36);  total (41/36)

TGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTT  >  NC_000913/1544276‑1544506                                                                                                             |                                                                                                                           tGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  1:479801/1‑139 (MQ=255) tGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  1:924383/1‑139 (MQ=255)        aTGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCa                                                                                                                 <  2:161720/113‑1 (MQ=255)        aTGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCa                                                                                                                 >  1:161720/1‑113 (MQ=255)        aTGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCa                                                                                       >  2:193116/1‑139 (MQ=255)          gTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                     <  2:924383/139‑1 (MQ=255)                  tCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              <  2:975427/122‑1 (MQ=255)                  tCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  1:975427/1‑122 (MQ=255)                  tCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                             <  1:12093/139‑1 (MQ=255)                     aCCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                          >  1:225038/1‑139 (MQ=255)                                 gTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTa                                                               >  2:161357/1‑138 (MQ=255)                                 gTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTa                                                               <  1:161357/138‑1 (MQ=255)                                 gTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTAt                                                              >  1:47235/1‑139 (MQ=255)                                   gAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGAGGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTg                                                            >  1:792099/1‑139 (MQ=255)                                     aTCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGcc                                                          >  1:641313/1‑139 (MQ=255)                                        tAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  1:295697/1‑100 (MQ=255)                                        tAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              <  2:295697/100‑1 (MQ=255)                                         aaaGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAt                                                      <  2:122007/139‑1 (MQ=255)                                          aaGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAtt                                                     >  1:482032/1‑139 (MQ=255)                                               ttACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  2:1092480/1‑93 (MQ=255)                                               ttACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              <  1:1092480/93‑1 (MQ=255)                                                tACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               <  2:1182876/139‑1 (MQ=255)                                                tACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               >  1:160644/1‑139 (MQ=255)                                                tACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               >  2:279908/1‑139 (MQ=255)                                                 aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                     >  1:424588/1‑100 (MQ=255)                                                 aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                     <  2:424588/100‑1 (MQ=255)                                                 aCATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaataa                                              >  2:422349/1‑139 (MQ=255)                                                  cATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                             <  2:792099/139‑1 (MQ=255)                                                  cATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTaataat                                             <  1:193116/139‑1 (MQ=255)                                                       tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTc                                                                       <  2:569597/108‑1 (MQ=255)                                                       tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTc                                                                       >  1:569597/1‑108 (MQ=255)                                                       tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGtta                                        <  1:351226/139‑1 (MQ=255)                                                       tAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGtta                                        <  2:225038/139‑1 (MQ=255)                                                        aaCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                       >  2:335962/1‑139 (MQ=255)                                                        aaCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                       <  1:279908/139‑1 (MQ=255)                                                         aCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                      >  1:854412/1‑139 (MQ=255)                                                            atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGt                                                                                 >  1:1110507/1‑93 (MQ=255)                                                            atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGt                                                                                 <  2:1110507/93‑1 (MQ=255)                                                            atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACa                                   <  2:854412/139‑1 (MQ=255)                                                             tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACag                                  >  1:887688/1‑139 (MQ=255)                                                                   caacaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAGCTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                            >  2:977976/1‑139 (MQ=255)                                                                   caacaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                            <  2:641313/139‑1 (MQ=255)                                                                   caacaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                            <  2:479801/139‑1 (MQ=255)                                                                   caacaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                            >  1:1174637/1‑139 (MQ=255)                                                                     acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              <  1:564809/71‑1 (MQ=255)                                                                     acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                              >  2:564809/1‑71 (MQ=255)                                                                       aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAg                                                                                      <  1:363694/77‑1 (MQ=255)                                                                       aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAg                                                                                      >  2:363694/1‑77 (MQ=255)                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTg                                                        <  1:171867/105‑1 (MQ=255)                                                                         gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTg                                                        >  2:171867/1‑105 (MQ=255)                                                                           acAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAt                    <  1:422349/138‑1 (MQ=255)                                                                            cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAGCTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAtt                   <  1:977976/139‑1 (MQ=255)                                                                                   ttATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               <  1:448622/104‑1 (MQ=255)                                                                                   ttATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               >  2:448622/1‑104 (MQ=255)                                                                                    tataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTaa           >  2:504846/1‑139 (MQ=255)                                                                                     ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                      <  2:629064/111‑1 (MQ=255)                                                                                     ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttatt                                      >  1:629064/1‑111 (MQ=255)                                                                                        gAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATc                                                                                             <  2:176115/53‑1 (MQ=255)                                                                                        gAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATc                                                                                             >  1:176115/1‑53 (MQ=255)                                                                                           ttAGTGATGATTTATTCATTAAAATGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTc                                                                       >  1:826439/1‑72 (MQ=255)                                                                                           ttAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTCTCCAGAAGGTCTATCGGTTc                                                                       <  2:826439/72‑1 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAGGGTCTAt                                                                             >  1:805552/1‑64 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                             <  2:805552/64‑1 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                           <  2:188083/82‑1 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                           >  1:188083/1‑82 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               <  2:684484/94‑1 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                               >  1:684484/1‑94 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt  >  1:306525/1‑139 (MQ=255)                                                                                             aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt  >  2:1070306/1‑139 (MQ=255)                                                                                                   gatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                     >  1:1086938/1‑66 (MQ=255)                                                                                                   gatTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                     <  2:1086938/66‑1 (MQ=255)                                                                                                     tttATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                           <  2:1195603/74‑1 (MQ=255)                                                                                                     tttATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                           >  1:1195603/1‑74 (MQ=255)                                                                                                         ttCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                     <  2:1218759/44‑1 (MQ=255)                                                                                                         ttCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                     >  1:1218759/1‑44 (MQ=255)                                                                                                         ttCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAtt                   >  1:315837/1‑110 (MQ=255)                                                                                                         ttCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAtt                   <  2:315837/110‑1 (MQ=255)                                                                                                             |                                                                                                                           TGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTT  >  NC_000913/1544276‑1544506

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