Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,500,060	C→T	pseudogene (475/942 nt)	yjgX ←	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,500,060	0	C	T	100.0%	189.9 / NA	59	pseudogene (475/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/30);  total (29/30)

AATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTG  >  NC_000913/4499925‑4500193                                                                                                                                        |                                                                                                                                      aaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                                    <  1:329920/139‑1 (MQ=255)       gTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAtt                                                                                                                              <  2:490242/139‑1 (MQ=255)          aGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTgg                                                                                                                           >  2:1542951/1‑139 (MQ=255)             aCCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                        <  2:472702/139‑1 (MQ=255)             aCCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                        <  2:513160/139‑1 (MQ=255)                    aaGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAccc                                                                                                                 >  2:475884/1‑139 (MQ=255)                                  gCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGaa                                                                                                   >  1:1804658/1‑139 (MQ=255)                                      gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAtt                                                                                               <  1:1542951/139‑1 (MQ=255)                                               gTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGga                                                                                      <  2:1754646/139‑1 (MQ=255)                                                ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                                    <  1:870166/92‑1 (MQ=255)                                                ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                                    >  2:870166/1‑92 (MQ=255)                                                ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGgag                                                                                     >  2:1357343/1‑139 (MQ=255)                                                         ttACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTAt                                                                            >  2:1703010/1‑139 (MQ=255)                                                           aCTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                        >  2:1737129/1‑93 (MQ=255)                                                           aCTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                        <  1:1737129/93‑1 (MQ=255)                                                               aTGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCa                                                                      >  2:1408211/1‑139 (MQ=255)                                                                   tGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  <  2:592170/139‑1 (MQ=255)                                                                   tGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  <  1:1357343/139‑1 (MQ=255)                                                                    ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGAt                                                                                                                         <  2:1098289/83‑1 (MQ=255)                                                                    ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  >  1:1379001/1‑138 (MQ=255)                                                                    ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  <  2:1379001/138‑1 (MQ=255)                                                                    ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCa                                                                 >  2:47381/1‑139 (MQ=255)                                                                    ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGATACACGGATTTGGAt                                                                                                                         >  1:1098289/1‑83 (MQ=255)                                                                     gCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                                    <  2:15217/71‑1 (MQ=255)                                                                     gCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                                    >  1:15217/1‑71 (MQ=255)                                                                          gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                   >  1:236980/1‑83 (MQ=255)                                                                          gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                   <  2:236980/83‑1 (MQ=255)                                                                           tAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                  >  1:662287/1‑115 (MQ=255)                                                                           tAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                  <  2:662287/115‑1 (MQ=255)                                                                              cAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGttt                                                       <  2:147400/139‑1 (MQ=255)                                                                                ggCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGGTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                                                     <  1:1703010/139‑1 (MQ=255)                                                                                   aaTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                   >  1:648324/1‑74 (MQ=255)                                                                                   aaTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                   <  2:648324/74‑1 (MQ=255)                                                                                     tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                  <  1:484025/105‑1 (MQ=255)                                                                                     tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                  >  2:484025/1‑105 (MQ=255)                                                                                     tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                                             <  2:333890/126‑1 (MQ=255)                                                                                     tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                                             >  1:333890/1‑126 (MQ=255)                                                                                     tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaa                                                >  2:675400/1‑139 (MQ=255)                                                                                                tACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  >  1:860265/1‑110 (MQ=255)                                                                                                tACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                  <  2:860265/110‑1 (MQ=255)                                                                                                  cTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCa                                   <  2:951923/139‑1 (MQ=255)                                                                                                    gCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTCCTGTTTCCTGTAACAACAGTTTGTCATc                                 <  2:1162365/139‑1 (MQ=255)                                                                                                     cTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCt                                >  1:1008612/1‑139 (MQ=255)                                                                                                          ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                                       <  2:658160/95‑1 (MQ=255)                                                                                                          ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                                       >  1:658160/1‑95 (MQ=255)                                                                                                          ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGAGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTGTGTCATCTGtatt                           >  1:1242430/1‑139 (MQ=255)                                                                                                               ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCg                                                                                        <  2:1572229/73‑1 (MQ=255)                                                                                                               ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCg                                                                                        >  1:1572229/1‑73 (MQ=255)                                                                                                               ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACggg                      >  1:1060654/1‑139 (MQ=255)                                                                                                                  ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaag                   <  1:675400/139‑1 (MQ=255)                                                                                                                    cggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGga                                                                                      <  1:551552/66‑1 (MQ=255)                                                                                                                    cggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                                                     <  1:1195659/99‑1 (MQ=255)                                                                                                                    cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGga                                                                                      >  2:551552/1‑70 (MQ=255)                                                                                                                    cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                                                     >  2:1195659/1‑103 (MQ=255)                                                                                                                    cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcg                 >  1:1709876/1‑139 (MQ=255)                                                                                                                     ggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcga                >  2:1065076/4‑139 (MQ=255)                                                                                                                        tAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcgaagc             <  1:475884/139‑1 (MQ=255)                                                                                                                        tAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaagcgaagc             <  2:1008612/139‑1 (MQ=255)                                                                                                                             gTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCt        >  1:445971/1‑139 (MQ=255)                                                                                                                             gTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCt        <  1:1065076/139‑1 (MQ=255)                                                                                                                                aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATc                                                                           >  1:1208461/1‑69 (MQ=255)                                                                                                                                aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATc                                                                           <  2:1208461/69‑1 (MQ=255)                                                                                                                                 aGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCt    >  2:183453/1‑139 (MQ=255)                                                                                                                                   gCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTg  <  2:445971/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      AATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTG  >  NC_000913/4499925‑4500193

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