Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,944	G→A	pseudogene (1845/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,944	0	G	A	100.0%	234.2 / NA	72	pseudogene (1845/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/33);  total (39/33)

ACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCC  >  NC_000913/2210815‑2211077                                                                                                                                  |                                                                                                                                      aCCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCg                                                                                                                              <  2:1440881/139‑1 (MQ=255)    ttCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                           >  1:215258/1‑139 (MQ=255)    ttCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAAt                                                                                                                           <  2:1120108/139‑1 (MQ=255)            gCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCAc                                                                                                                   >  2:391502/1‑139 (MQ=255)                tttGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACcgcg                                                                                                               <  1:655918/139‑1 (MQ=255)                 ttGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                              <  2:150471/139‑1 (MQ=255)                            gTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                   <  2:215258/139‑1 (MQ=255)                                       atcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                        >  1:324563/1‑139 (MQ=255)                                                    aCGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                           <  1:1155855/139‑1 (MQ=255)                                                    aCGCGGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                           <  1:313660/139‑1 (MQ=255)                                                          ccAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGt                                                                     <  1:1357270/139‑1 (MQ=255)                                                          ccAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGt                                                                     >  1:12475/1‑139 (MQ=255)                                                                gTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACa                                                               >  2:1409645/1‑139 (MQ=255)                                                                  tACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGa                                                             >  2:80415/1‑139 (MQ=255)                                                                            aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   >  1:1514467/1‑139 (MQ=255)                                                                            aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   >  1:1260582/1‑139 (MQ=255)                                                                            aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   <  2:12475/139‑1 (MQ=255)                                                                            aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   <  1:391502/139‑1 (MQ=255)                                                                              tGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTCGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGGTGGGACAATTTCGTAATACAGAAAGCGAGGCACg                                                 >  2:1139892/1‑139 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGc                                                                                                                                   >  2:1085962/1‑56 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGc                                                                                                                                   <  1:1085962/56‑1 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                              >  1:314693/1‑109 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                              <  2:314693/109‑1 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   <  1:725615/136‑1 (MQ=255)                                                                               gCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCa                                                   >  2:725615/1‑136 (MQ=255)                                                                                   cATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAg                                            <  2:854214/139‑1 (MQ=255)                                                                                   cATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAg                                            <  2:1514467/139‑1 (MQ=255)                                                                                     tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                   >  2:232761/1‑82 (MQ=255)                                                                                     tAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTAt                                                                                                   <  1:232761/82‑1 (MQ=255)                                                                                          gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               >  2:877942/1‑81 (MQ=255)                                                                                          gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               <  1:877942/81‑1 (MQ=255)                                                                                          gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                     >  2:837171/1‑139 (MQ=255)                                                                                           cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTgg                                                                                                                                <  2:1002009/47‑1 (MQ=255)                                                                                           cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTgg                                                                                                                                >  1:1002009/1‑47 (MQ=255)                                                                                           cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATg                                    >  1:877142/1‑139 (MQ=255)                                                                                           cTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATg                                    <  2:324563/139‑1 (MQ=255)                                                                                             ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                        <  1:1542131/85‑1 (MQ=255)                                                                                             ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                        >  2:1542131/1‑85 (MQ=255)                                                                                             ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                  >  2:694150/1‑139 (MQ=255)                                                                                                gcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATaa                               >  2:710122/1‑139 (MQ=255)                                                                                                gcTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATaa                               <  1:837171/139‑1 (MQ=255)                                                                                                  tATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACg                             <  1:1619096/139‑1 (MQ=255)                                                                                                   aTCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACgg                            >  2:1424222/1‑139 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCa                                                                                                                    <  2:539157/40‑1 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCa                                                                                                                    >  1:539157/1‑40 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                           >  1:270484/1‑81 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                           <  2:270484/81‑1 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAAtt                                                                         <  2:1155546/83‑1 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAAtt                                                                         >  1:1155546/1‑83 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                 >  2:850936/1‑139 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                 >  1:1556592/1‑139 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                 >  1:721130/1‑139 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                 >  2:378247/1‑139 (MQ=255)                                                                                                              agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                 >  2:1571014/1‑139 (MQ=255)                                                                                                                tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               <  2:1648609/59‑1 (MQ=255)                                                                                                                tGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               >  1:1648609/1‑59 (MQ=255)                                                                                                                  cTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGa             >  2:46506/1‑139 (MQ=255)                                                                                                                     ttGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                >  1:269021/1‑53 (MQ=255)                                                                                                                     ttGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                <  2:269021/53‑1 (MQ=255)                                                                                                                     ttGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa            <  2:758972/137‑1 (MQ=255)                                                                                                                     ttGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa            >  1:758972/1‑137 (MQ=255)                                                                                                                      tGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACaa         <  1:1571014/139‑1 (MQ=255)                                                                                                                        gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               >  2:889953/1‑51 (MQ=255)                                                                                                                        gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGc                                                                                               <  1:889953/51‑1 (MQ=255)                                                                                                                            ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGtttt                      >  2:404062/1‑120 (MQ=255)                                                                                                                            ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGtttt                      <  1:404062/120‑1 (MQ=255)                                                                                                                            ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTc   >  2:1163938/1‑139 (MQ=255)                                                                                                                             gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                        <  1:460160/53‑1 (MQ=255)                                                                                                                             gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                        >  2:460160/1‑53 (MQ=255)                                                                                                                             gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTcc  >  1:1358327/1‑139 (MQ=255)                                                                                                                              cATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTa                                                                                                    <  1:1632475/40‑1 (MQ=255)                                                                                                                              cATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTa                                                                                                    >  2:1632475/1‑40 (MQ=255)                                                                                                                                  |                                                                                                                                      ACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCC  >  NC_000913/2210815‑2211077

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: ‑