Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,788,238	G→A	pseudogene (1842/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,788,238	0	G	A	100.0%	215.4 / NA	66	pseudogene (1842/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/34);  total (32/34)

AGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTAGTCGTGG  >  NC_000913/2788112‑2788369                                                                                                                               |                                                                                                                                    aGCAATGGCATATATAGCCGTTCTGTCGTTGAGCAGGCCTTTAGTTTATTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTCTGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTAAATAAAGGt                                                                                                                         <  2:418409/139‑1 (MQ=255)      tGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    <  2:309585/139‑1 (MQ=255)      tGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    <  1:495939/139‑1 (MQ=255)         cATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTAc                                                                                                                 <  1:660074/139‑1 (MQ=255)                aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                          >  2:494272/1‑139 (MQ=255)                aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                          >  1:657597/1‑139 (MQ=255)                aGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                          <  2:657597/139‑1 (MQ=255)                  ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                                       >  1:61005/1‑124 (MQ=255)                  ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTg                                                                                                                       <  2:61005/124‑1 (MQ=255)                  ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagca                                                                                                        >  1:24284/1‑139 (MQ=255)                  ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagca                                                                                                        >  2:394223/1‑139 (MQ=255)                  ccGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagca                                                                                                        >  2:789673/1‑139 (MQ=255)                   cGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    <  2:619200/126‑1 (MQ=255)                   cGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    >  1:619200/1‑126 (MQ=255)                   cGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcagcag                                                                                                       <  2:504640/139‑1 (MQ=255)                     ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAg                                                                                                                           >  1:142681/1‑117 (MQ=255)                     ctctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAg                                                                                                                           <  2:142681/117‑1 (MQ=255)                       ctGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGctct                                                                                                   >  1:215640/1‑139 (MQ=255)                            gTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAAt                                                                                              >  2:715892/1‑139 (MQ=255)                                  cagGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTc                                                                                        >  1:550920/1‑139 (MQ=255)                                   agGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                                       <  2:215640/139‑1 (MQ=255)                                   agGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTcc                                                                                       <  1:494272/139‑1 (MQ=255)                                    ggCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                      <  2:24284/139‑1 (MQ=255)                                     gCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     <  1:715892/139‑1 (MQ=255)                                       cTTTAGTTTACTGGGATGGCCAAGTATTAAAATCGACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                                   <  1:538172/139‑1 (MQ=255)                                        tttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                                  <  1:789673/139‑1 (MQ=255)                                        tttAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                                  <  2:550920/139‑1 (MQ=255)                                                 cTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                          <  2:788402/106‑1 (MQ=255)                                                 cTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATcag                                                                                                          >  1:788402/1‑106 (MQ=255)                                                    ggATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTGCAATGAATTCCACGCATTATATATAGAt                                                                      >  1:503271/1‑139 (MQ=255)                                                      aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                                  >  1:795525/1‑125 (MQ=255)                                                      aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCa                                                                                  <  2:795525/125‑1 (MQ=255)                                                      aTGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATCATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGAtca                                                                    <  2:503271/139‑3 (MQ=255)                                                         gCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCAGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    <  2:58166/88‑1 (MQ=255)                                                         gCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCAGCGGTGAGGCGTAACTTAATAAAGGTTGctc                                                                                                                    >  1:58166/1‑88 (MQ=255)                                                                         aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAg                                                                                                                           <  2:66753/65‑1 (MQ=255)                                                                         aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAg                                                                                                                           >  1:66753/1‑65 (MQ=255)                                                                         aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     >  1:426611/1‑103 (MQ=255)                                                                         aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     <  2:426611/103‑1 (MQ=255)                                                                         aaCGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCa                                                 >  2:86685/1‑139 (MQ=255)                                                                                 aaaTTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAggg                                         >  2:17225/1‑139 (MQ=255)                                                                                                     tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     >  2:220778/1‑75 (MQ=255)                                                                                                     tCCACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     <  1:220778/75‑1 (MQ=255)                                                                                                       cACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                                  <  2:305786/92‑1 (MQ=255)                                                                                                       cACGCGTGTCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAg                                                                  >  1:305786/1‑92 (MQ=255)                                                                                                              gtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     >  2:96295/1‑66 (MQ=255)                                                                                                              gtCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCAc                                                                                     <  1:96295/66‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                      >  2:551305/1‑64 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCa                                                                                      <  1:551305/64‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                                   >  1:215567/1‑67 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGc                                                                                   <  2:215567/67‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATtatat                                                                            <  1:764544/74‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATtatat                                                                            >  2:764544/1‑74 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGa                                                                       >  2:310674/1‑79 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGa                                                                       <  1:310674/79‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATaa                                                                    <  1:44094/82‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATaa                                                                    >  2:44094/1‑82 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCGGGAATATGGAGCGCATTAGAGGGTGTATCTGCAACTAACAACGTTCTGGCCtt           >  2:228520/1‑139 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGgt                                        <  1:683138/110‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGgt                                        >  2:683138/1‑110 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTc                  <  2:788334/132‑1 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTc                  >  1:788334/1‑132 (MQ=255)                                                                                                               tCGGCGGTGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCtt           >  1:448994/1‑139 (MQ=255)                                                                                                                      tGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCAttt                                              <  1:325455/97‑1 (MQ=255)                                                                                                                      tGAGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCAttt                                              >  2:325455/1‑97 (MQ=255)                                                                                                                        aGGCGTAACTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTAGTCGTgg  <  1:86685/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                    AGCAATGGCATATATAGCCGCTCTGCCGTTCAGCAGGCCTTTAGTTTACTGGGATGGCCAAGTATTAAAATCAACGACTGAAATTAAAAATTATGATGGGTCCACGCGTGTCGGCGGTGAGGCGTAGCTTAATAAAGGTTGCTCTACCTATCAGCAGCTCTACAATGAATTCCACGCATTATATATAGATAAAGCAGGAATATGGAGCGCATTTGAGGGTGTATCTGCAACTAAAAACGTTCTGGCCTTTAGTCGTGG  >  NC_000913/2788112‑2788369

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