Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,209,276	G→A	pseudogene (177/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,209,276	0	G	A	100.0%	194.3 / NA	50	pseudogene (177/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/27);  total (23/27)

GACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCGGCAAT  >  NC_000913/2209132‑2209422                                                                                                                                                 |                                                                                                                                                   gACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAAc                                                                                                                                              <  2:731434/151‑1 (MQ=255)  acacCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACg                                                                                                                                             >  1:151344/1‑151 (MQ=255)    acCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTAACGGTGTCCGTACCCAGCTAGCGAACGGc                                                                                                                                           >  1:62888/1‑151 (MQ=255)      cGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCa                                                                                                                                         <  2:89097/151‑1 (MQ=255)       gCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACg                                                                                                                                             <  2:28604/146‑1 (MQ=255)             ccTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTg                                                                                                                                  <  1:785447/151‑1 (MQ=255)                gACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGaa                                                                                                                                <  2:737770/150‑1 (MQ=255)                gACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGaa                                                                                                                                >  2:87795/1‑150 (MQ=255)                gACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGa                                                                                                                                 <  1:276019/149‑1 (MQ=255)                  cGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGaa                                                                                                                                >  2:509532/1‑148 (MQ=255)                            cAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGc                                                                                                                   <  2:767732/151‑1 (MQ=255)                                   ggTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCg                                                                                                            >  1:276068/1‑151 (MQ=255)                                   ggTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCg                                                                                                            <  1:574626/151‑1 (MQ=255)                                      tttGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCggg                                                                                                          <  1:384028/150‑1 (MQ=255)                                        tGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGcc                                                                                                       >  1:775879/1‑151 (MQ=255)                                           tAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAAc                                                                                                    >  1:575095/1‑151 (MQ=255)                                            aaGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACg                                                                                                   <  1:757622/151‑1 (MQ=255)                                             aGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACgg                                                                                                  >  1:776232/1‑151 (MQ=255)                                                gcAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTAt                                                                                               >  2:322756/1‑151 (MQ=255)                                                 cAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTAtg                                                                                              >  1:839125/1‑151 (MQ=255)                                                    tCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGc                                                                                           >  2:610475/1‑151 (MQ=255)                                                     cTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGcc                                                                                          >  1:48228/1‑151 (MQ=255)                                                        aaGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTc                                                                                       <  1:513259/151‑1 (MQ=255)                                                          ggAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCAt                                                                                     <  2:659990/151‑1 (MQ=255)                                                                ggAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGt                                                                                <  2:518770/150‑1 (MQ=255)                                                                  aaaaTGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGAt                                                                             <  2:438918/151‑1 (MQ=255)                                                                      tGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCt                                                                          <  2:627727/150‑1 (MQ=255)                                                                       ggCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTgg                                                                        >  1:806580/1‑151 (MQ=255)                                                                          aaCGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGgtgt                                                                     >  1:106323/1‑151 (MQ=255)                                                                               cccGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGt                                                                <  2:48138/151‑1 (MQ=255)                                                                                ccGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGtt                                                               >  2:796529/1‑151 (MQ=255)                                                                                   gagaTCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATc                                                            >  1:658840/1‑151 (MQ=255)                                                                                            cATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACttt                                                   >  2:163692/1‑151 (MQ=255)                                                                                            cATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACttt                                                   <  2:504157/151‑1 (MQ=255)                                                                                               gagaACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTtgt                                                 <  1:829195/150‑1 (MQ=255)                                                                                                  aaCGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGcc                                              >  2:383402/1‑150 (MQ=255)                                                                                                     gacgCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGccacca                                          <  2:335579/151‑1 (MQ=255)                                                                                                          tttAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTc                                       >  2:413218/1‑149 (MQ=255)                                                                                                             aaTCGCCACCTTCAGTGACGGTGTCCGTCCCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCCCTCAGTCt                                  <  2:271851/151‑1 (MQ=255)                                                                                                                     ccTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGaaaa                          <  2:362252/151‑1 (MQ=255)                                                                                                                     ccTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGaaaa                          <  1:206402/151‑1 (MQ=255)                                                                                                                      cTTCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGaaaa                          >  1:859222/1‑150 (MQ=255)                                                                                                                        tCAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGaaaaa                         <  2:1041/149‑1 (MQ=255)                                                                                                                         cAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAaagaa                      <  1:574479/151‑1 (MQ=255)                                                                                                                         cAGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAaaga                       >  1:165046/1‑150 (MQ=255)                                                                                                                          aGTGACGGTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAaaga                       >  1:851299/1‑149 (MQ=255)                                                                                                                                ggTGTCCGTACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGt               <  2:711107/151‑1 (MQ=255)                                                                                                                                 gtgtCCGCACCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTg              >  1:863769/1‑151 (MQ=255)                                                                                                                                         aCCCAGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCgg      <  1:644317/151‑1 (MQ=255)                                                                                                                                             aGCTAGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCGGCAAt  <  1:491547/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCGGAGATCAGCCATGAGAACGACGCTTTAATCGCCACCTTCAGTGACGGTGTCCGTACCCAGCTGGCGAACGGCCAGGCACTGAAAGAGGCTCAGTGCAGTTGCGGGGCCAACGGTATGTGCCGTCATCGCGTGATGCTGGTGTTAAGTTATCAACGACTTTGTGCCACCACTCAGTCTACGGAAAAAGAAGAAGAGTGGGATCCGGCAAT  >  NC_000913/2209132‑2209422

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