Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	658,216	C→A	pseudogene (186/788 nt)	ybeM →	deamidated glutathione amidase, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	658,216	0	C	A	100.0%	209.6 / NA	60	pseudogene (186/788 nt)	ybeM	deamidated glutathione amidase, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/26);  total (34/26)

ACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCG  >  NC_000913/658080‑658349                                                                                                                                         |                                                                                                                                      aCGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGaa                                                                                                                                    <  1:203524/140‑1 (MQ=255)   gCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAAtt                                                                                                                                  >  1:199367/1‑140 (MQ=255)   gCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAAtt                                                                                                                                  >  1:358413/1‑140 (MQ=255)        gaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCg                                                                                                                             <  1:315982/140‑1 (MQ=255)             gtgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACtt                                                                                                                        >  2:89775/1‑140 (MQ=255)               gCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTa                                                                                                                      >  2:154948/1‑140 (MQ=255)                 cTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTa                                                                                                                      >  2:306821/1‑138 (MQ=255)                 cTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTa                                                                                                                      <  1:306821/138‑1 (MQ=255)                 cTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACg                                                                                                                    >  1:905833/1‑140 (MQ=255)                 cTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACg                                                                                                                    <  2:675068/140‑1 (MQ=255)                   cGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGc                                                                                                                  >  2:122270/1‑140 (MQ=255)                   cGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGc                                                                                                                  >  2:762017/1‑140 (MQ=255)                    gTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCg                                                                                                                 <  2:358413/140‑1 (MQ=255)                      tGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCgag                                                                                                               <  1:620586/140‑1 (MQ=255)                           gcgcAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGt                                                                                                          <  2:789250/140‑1 (MQ=255)                               aGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAAc                                                                                                      >  1:26490/1‑140 (MQ=255)                                ggcggcGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACg                                                                                                     <  2:26490/140‑1 (MQ=255)                                ggcggcGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACg                                                                                                     >  1:868159/1‑140 (MQ=255)                                 gcggcgGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGt                                                                                                    <  1:626965/140‑1 (MQ=255)                                       gAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACatg                                                                                              <  1:504818/140‑1 (MQ=255)                                           acgacgCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACAtgatga                                                                                          >  1:285028/1‑140 (MQ=255)                                           acgacgCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACAtgatga                                                                                          >  2:342354/1‑140 (MQ=255)                                              acgCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATgacga                                                                                       <  2:436678/140‑1 (MQ=255)                                                   tCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAAtt                                                                                  <  2:868159/140‑1 (MQ=255)                                                      cTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTg                                                                               <  1:27478/140‑1 (MQ=255)                                                       tGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGa                                                                              <  2:426965/140‑1 (MQ=255)                                                            gCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGAtt                                                                         <  1:389874/140‑1 (MQ=255)                                                                 gCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGt                                                                    >  2:541767/1‑140 (MQ=255)                                                                 gCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGAAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGCCGACAATTCTGACGATTCATGt                                                                    >  2:200973/1‑140 (MQ=255)                                                                  ccGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGtt                                                                   <  1:666428/140‑1 (MQ=255)                                                                          aTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAAc                                                           >  2:781105/1‑140 (MQ=255)                                                                                  cgcgcgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgc                                                   <  2:269026/140‑1 (MQ=255)                                                                                   gcgcgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcg                                                  <  1:342354/140‑1 (MQ=255)                                                                                      cgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTTATGTTCCTTCAACGCCGGGgcgcgc                                               >  2:510992/1‑140 (MQ=255)                                                                                      cgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               <  2:199367/140‑1 (MQ=255)                                                                                      cgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               <  1:89775/140‑1 (MQ=255)                                                                                      cgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               >  2:542522/1‑140 (MQ=255)                                                                                      cgATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               >  2:583141/1‑140 (MQ=255)                                                                                        atgatCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCAt                                             <  2:905833/140‑1 (MQ=255)                                                                                                 atgCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               <  2:156560/129‑1 (MQ=255)                                                                                                 atgCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcgcgc                                               >  1:156560/1‑129 (MQ=255)                                                                                                 atgCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGc                                    >  2:255431/1‑140 (MQ=255)                                                                                                     aGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggt                                >  1:558851/1‑140 (MQ=255)                                                                                                       aTCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGGCAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggtgg                              >  1:27634/1‑140 (MQ=255)                                                                                                       aTCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggtgg                              >  2:60328/1‑140 (MQ=255)                                                                                                       aTCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggtgg                              >  1:656366/1‑140 (MQ=255)                                                                                                       aTCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggtgg                              >  2:409137/1‑140 (MQ=255)                                                                                                       aTCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCtggtgg                              >  2:471402/1‑140 (MQ=255)                                                                                                        tCTATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGc                             >  1:303170/1‑140 (MQ=255)                                                                                                          tATCGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCAc                           >  1:302896/1‑140 (MQ=255)                                                                                                             cGGTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTc                        <  1:72205/140‑1 (MQ=255)                                                                                                               gTTAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcag                      <  1:583141/140‑1 (MQ=255)                                                                                                                 tAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcaggc                    >  1:691975/1‑140 (MQ=255)                                                                                                                 tAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcaggc                    >  1:337300/1‑140 (MQ=255)                                                                                                                 tAAATCAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcaggc                    >  1:308411/1‑140 (MQ=255)                                                                                                                      cAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGaa               >  2:885205/1‑140 (MQ=255)                                                                                                                      cAGCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGaa               >  2:9556/1‑140 (MQ=255)                                                                                                                        gCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcagg                     >  2:391133/1‑132 (MQ=255)                                                                                                                        gCACAGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTcagg                     <  1:391133/132‑1 (MQ=255)                                                                                                                          acaGCTGCTGGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACAt           >  1:800449/1‑140 (MQ=255)                                                                                                                              ctgctgGAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACAtcgtc       >  1:532805/1‑140 (MQ=255)                                                                                                                                   ggAAGGAGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCg  <  1:781105/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      ACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCGCCCG  >  NC_000913/658080‑658349

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